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ABSTRACT 


Tables and graphs of the number densities of hydrogen atoms (H} t protons (H + ), free 
electrons (e - ), hydrogen molecules {Hg), negative hydrogen ions {H")> hydrogen diatomic 
molecular ions and hydrogen triatomic molecular ions (H^) in a hydrogen gas or 

plasma in thermodynamic equilibrium were computed for temperatures from 300 to 
40 000 K and pressures from 1. 01325X10 5 to 1. 01325xl0 B N/m 2 (1 to 1000 atm) except 
for a small, high -density region. For temperatures above 2000 K, the Debye- Hiickel 
approximation was employed,. Morse potentials are given for all singly excited attractive 
states of with principal quantum numbers to 17. A FORTRAN IV computer program 
for computing number densities and other properties of a hydrogen plasma is included. 
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COMPONENTS OF A HYDROGEN PLASMA INCLUDING MINOR SPECIES 

by R. W. Patch 
Lewis Research Center 


SUMMARY 

Tables and graphs of the number densities of hydrogen atoms (H), protons (H + ), 
free electrons (e~), hydrogen molecules (Hg), negative hydrogen ions (H~), hydrogen 
diatomic molecular ions (Hg), and hydrogen triatomic molecular ions {Hg) in a hydrogen 
gas or plasma in thermodynamic equilibrium were computed for temperatures from 300 
to 40 000 K and pressures from 1. 01325x10^ to 1. 01325x10^ newtons per square meter 
{1 to 1000 atm) except for a small, high -density region. For temperatures above 2000 K, 
the Bebye-Huckel approximation was employed. The ion Hg had not been included in 
previous calculations and was found to cause large increases in e" and H" number den- 
sities at small degrees of ionization. The results presented in this report are believed 
to be the most accurate available. They are particularly applicable to calculating optical 
absorption coefficients. Also given are Morse potentials for all singly excited attractive 
states of Hg with principal quantum numbers to 17. 

A computer program in FORTRAN IV for computing number densities and other 
properties of a hydrogen plasma with or without coulomb interactions between free par- 
ticles Is presented. Sample problems and a flow diagram are also included. 


INTRODUCTION 

In high -temperature propulsion devices such as gaseous - core nuclear rockets, the 
dominant mechanism of heat transfer is radiant energy exchange between volumes of 
plasma and between the plasma and the wall. To calculate such heat transfer, it is 
necessary to know the opacity of the plasma. For a hydrogen plasma under many condi- 
tions, minor species make major contributions to the opacity, so their number densities 
must be known. 
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In past work (refs, 1 to 6) a wide variety of assumptions has been made, some of 
which were later shown to be invalid, and in no case were the hydrogen triatomic molec- 
ular ion (Hg) or excited electronic states of the hydrogen molecule (Hg) (designated Hg in 
this report) included. This investigation was conducted to provide more accurate esti- 
mates of number densities for hydrogen for temperatures between 300 and 40 000 K and 
pressures between 1. 01325x10^ and 1. 01325x10^ newtons per square meter (1 to 1000 atm). 
The species included above 2000 K were the hydrogen atom (H), the proton (H + ), the 
free electron (e"), the negative hydrogen ion (H~), the hydrogen diatomic molecular ion 
(H g), Hg (including Hg), and Hg, It was not feasible to estimate virial coefficients for 
all the species, so all inter molecular forces involving neutral species were neglected. 

O 

This necessitated excluding a high -density region near 300 K and 10 newtons per square 
meter from the investigation. 

One of the objects of the investigation was to find the effect of including Hg on calcu- 
lated number densities of e“ and H". Once this is known, the effect of Hg on calculated 
electrical conductivity and opacity can be obtained. 


ANALYSIS 

The calculation of number densities over the wide range of 40 000 to 300 K neces- 
sitated three analyses to prevent major underflows in the digital computer program. The 
range from 40 000 to 2000 K was designated "high” temperature, 2000 to 1300 K was 
"medium" temperature, and 1300 to 300 K was "low" temperature. In all three analy- 
ses, values of the fundamental constants were obtained from reference 7. 


High Temperature (2000 to 40 000 K) 

The equilibrium components of the plasma were calculated in this range by mini- 
mizing the Gibbs free energy G of a closed, neutral system at constant temperature T 
and pressure p(ref. 8). (Symbols are given in appendix A. ) Hence, equations for Gibbs 
free energy, pressure, charge conservation, and proton conservation were required and 
are given in this section. The partition functions needed to calculate the Gibbs free 
energy are also given. 

Assumptio ns a nd restri ctions. - The assumptions and restrictions utilized in the 
high -temperature analysis are as follows; 

(1) The species H(H^ and H*), H + , e", Hg(H| and Hg), H", Hg, and Hg were included, 
where and h| are the ground electronic states of H and Hg, respectively, and H* is 
all excited electronic states of H collectively. 
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(2) The plasma was treated as an ideal gas except for coulomb interactions between 
free charged particles, which were treated by the Debye -Huckel approximation. 

(3) No shifts of energy levels of excited electronic states relative to ground electronic 
states were considered. {These shifts are observed to be small (refs. 9 and 10 ), and no 
one has correlated the observed shifts for H (ref. 9) with a theory whose assumptions 
would be valid for the conditions of the experiment. ) 

(4) Only singly excited electronic states of Hg were considered. (Most doubly excited 
electronic states of Hg have negligible number densities. ) 

(5) The free electrons were nondegenerate (i. e. , they had a Boltzmann distribution of 
velocities). 

(6) The contribution of photons to the pressure was neglected. 

(7) Nuclear spin was taken into account in a consistent manner (i. e. , it was not 
included in calculating degeneracies of atoms and atomic ions and was accounted for in 
molecules and molecular ions by using a symmetry number (ref. 11)). An additional 
refinement consistent with this method was included for h| and Hg. 

Gibbs free energy . - To obtain the Debye -Huckel approximation for the Gibbs free 
energy of a plasma, it is first necessary to know the reciprocal Debye length, the excess 
Helmholtz free energy, and the excess pressure. For a plasma the reciprocal Debye 
length is 



where z^ is the net number of elementary charges on the particle (1, 0, or -1) and the 
summation is over all free particles. The excess Helmholtz free energy A ex is the 
addition to the ideal -gas Helmholtz free energy resulting from coulomb interactions and, 
according to the Debye -Huckel theory (refs. 11 and 12), is 


A 


ex 


kTVK 3 

12jt 


( 2 ) 


Consequently, the excess pressure p is 


p 


ex 



(3) 


The excess Gibbs free energy is then 


3 


G ex A ex + p ex^ 


kTV/f 

8tt 
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(4) 


From reference 11 and equation (4) the total Gibbs free energy of a hydrogen plasma is 

7 


G = 

i=l 


kTVx 

Qir 


(5) 


where id is the chemical potential for species i in an ideal gas mixture, hi such 
an ideal mixture, the /i. id of one species is independent of the presence of all other 
species (ref. 11); thus, equation (5) becomes 


Vq i kTVx 3 


G = -kT V N. In 

U1 1 A i N i 


( 6 ) 


where 


A i s 


h" 


3/2 


^2?rra.kTy 


i = l,2,. . .,7 


(7) 


Equation (6) is the desired expression for the Gibbs free energy of a plasma. 

Pressure. - The pressure is found from the ideal -gas law and equation (3) and is 


P - y N . - 

xr 1 o/l„ 


i=l 


24jt 


( 8 ) 


Charge conservation. - For a neutral plasma 


N 3 + N 5 = N 2 +N 6 + N 7 


(9) 


Proton conservation . - In a closed system the number of protons rj’ must be con- 
served, where 


7j* = Nj + N 2 + 2N 4 + Ng + 2N 6 + 3N 7 


(10) 
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Equilibrium . - For thermodynamic equilibrium, the eight variables N. and V were 
simultaneously varied to minimize G (eq. (6)) subject to the constraints given by equa- 
tions (8) to (10). The method for this, given in appendix B, gives the same result as 
using equations (8) to (10) in conjunction with a dissociative equilibrium equation for each 
of the reactions 


H 2 ~2H 

Hg — H| + H 

H" - H + e" 

Hg -*■ H + H + 

and a Saha equation with lowering of the ionization potential (ref. 12) for the reaction 

H - H + + e" 

Inte rnal partitio n functions . - In carrying out the minimization of G, values for q^ 
are required in equation (6). These values must all be referred to the energies of two 
reference states. First, however, the partition functions q: were calculated using the 
ground-state energy of each species as the reference. 

For H atoms 




— 2 -f* 


17 

n=2 


n^n 


(ID 


where 


, -hcT /kT 

(p n = 2n z e e (l: 

and F n is a cutoff factor given in appendix C, The summation is carried out through a 
principal quantum number n of 17 because F n is 0 for larger n for all conditions in 
this report. For H + , ql, is 1. For e~, q^ is 2. The qj for Hg, H”, Hg, and Hg are 
given in subsequent paragraphs. 

For Hg there are contributions made by h| and Hg. The contribution of h| to the 
partition function of Hg neglecting nuclear spin is obtained from equation (6) of refer- 
ence 13 by substituting J for K and is 
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(13) 


Ik* = - E £ < 2J + l)e' hC[V ' <V ’ j3/kT 


Here g e is 1, 0 is 2, and the J summation includes metastable levels. For H| the 
rotational levels are spaced widely enough so the differences in nuclear spin degeneracies 
should be taken into account, at least near 300 K, Equation (13) then obviously becomes 


1h| - E £ v 

d V 


(14) 


where 


^v *" 


g e -hcG v /kT 


j£ (2J +1 )e- hcf<v ’ J )/ kT + i V (2J + l)e- hCf(v ' J)AT 


oddJ 


even J 


(15) 


Values of G v were obtained from reference 14. Values of f(v,J) were calculated by the 
Wentzel-Kramers-Brillouin approximation using the computer program in reference 13. 
The most reliable potential energy was used at each internuclear distance, namely, 

(1) from 0. 37x10" 10 to 0. 40x10" meter the Kolos and Wolniewicz potential (ref. 15) was 
used, (2) from 0. 41x10"^® to 3. 00x10”^ meter the Rydberg -Klein -Rees potential of 
Spindler (ref. 14) was used, (3) from 3. 01x10" 10 to 5, 29x10" 10 meter the Kolos and 
Wolniewicz potential was used, and (4) from 5. 55x10"^ to 25x10"^ meter the 
Hirschf eider and Lowdin potential (ref. 16) was used. Potential energies at internuclear 
distances not included in (1) to (4) were obtained by interpolation. 

The partition function for Hg was obtained by adding the contributions of Hg to 


q 4 " + £ F n y h 


(16) 


* 

where the summation is over all singly excited electronic states Hg with n of 17 or less, 
where n is the principal quantum number of the excited electron in the small- 
inter nuclear -distance designation (ref. 17) and where 


80 1/2 g e 




j^£PkTj 


3//2 -hcf/kT 


rJW o (0) 2r e W x (0) W n (0) 
— — + — 


«h 


(17) 
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and 
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ue 


-u/2 


u = 


1 - e 


-u 


heap 

kT 


9 = 


kT 

hcD 


( 18 ) 


Equation (17) is a classical partition function with a vibrational quantum correction for a 
Morse vibrating rotator. It was obtained from equation (5) in reference 18 by multiplying 
by g e and shifting the reference energy to that of the rotationless ground vibrational 
state of h|. The constants ¥, r e , 0, D, and were obtained as discussed in appen- 
dix D. Where available, experimental energy levels for n of 2 to 4 in the small- 
inter nuclear -distance designation were used in obtaining these constants (see appendix E 
and table I). Values of g fi , T, r e , and D are given in table II . 

The ground state of H" is a 1 So state (ref. 19). Electronically excited states were 
neglected, so had the value 1. 

The ground state of Hg is a state, whereas the ground state of Hg is a 

state. Because they both have the same nuclear spin statistics, qj, is obtainable (ref. 20) 

from equations (14) and (15) by substituting q£ for and K for J. To avoid con- 
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fusion later, £ v is replaced by 12 , so equations (14) and (15) become 


«6 = 

v 


(19) 


where 


12 = 
v 


g e -hcG/kTi 


S (2K+l)e 
oddK 


-hcf (v, K)/kT 


a £ 

2 


(2K+ l)e 


-hcf (v, K)/kT 


evenK 


( 20 ) 


The summations include m eta stable rotational levels. Excited electronic states of Hg 
were neglected. Values for G y were obtained from reference 21, and values for f 
were obtained from reference 13. 

The partition function for Hg was obtained from reference 13, where electron- 
ically excited states were neglected. 

The atom and the free electron e" were selected as the two reference states for 
the q- . The shifts in reference energies from the values used for the qj result in the 
following relations: 
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I 



n 


1 


i 


-hclf/kT 

q l = q l> q 2 = q 2 e ■ q 3 =q 3> 


hcwJ/kT 

q 4 = q 4 e 


hd 5 /kT -hc(lf-W fi )/kT -hc(lf-W 7 )/kT 

q 5 =q 5 e , q 6 * , d 7 = 


r (2D 


Note that the lowering of the ionization potential due to coulomb interactions was omitted 
from equation (21) as this effect is already included in the last term of equation (6). 
Values for l|, Ig, w|, Wg, and are given in table HI. 


Medium Temperature (1300 to 2000 K) 


In this range the high-temperature analysis was inapplicable because of major under- 
flows in the computer. The degree of ionization was so small that the coulomb interaction 
between free charged particles was negligible and consequently could be neglected. The 
number densities of H* and Eg were each less than 2x10'^ times the number density of 
Hg, so the first two could be neglected. The equilibrium components of the plasma were 
calculated by means of equilibrium constants, using partition functions from the previous 
section. 

Assumptions and restrictions . - The assumptions and restrictions utilized in the 
medium-temperature analysis are as follows: 

(1) Only the species H(H^ and H*), e~, Hg(H;!; and Hg), H", and Hg were included. 

(2) The plasma was treated as an ideal gas. 

(3) No shifts of energy levels of excited electronic states relative to ground elec- 
tronic states were considered. 

(4) Only singly excited electronic states of Hg were considered. 

(5) Nuclear spin was taken into account in the same way as in the high -temperature 
analysis. 

Solution. - Let p^ be an arbitrary reference pressure and define 


s’ 




i = l, 2, . . ., 7 


( 22 ) 


For the reaction 2H — Hg the equilibrium constant is (ref. 11) 


8 


(23) 


s’. 


3*2 


-Si- 


kTqWkTqA - 2 

ViVW 


hcw|/kT 
e ’ 


Define 


«i s Vl T 


3/2 



i = l f 2, 


, 7 


(24) 


Combining equations (23) and (24) gives 


Kpl = 



hcwJ/kT 
e “* 


(25) 


For medium temperatures and pressures between 1. 01 3 25x1 0 5 and 1. 01 32 5xl 0 8 newtons 
per square meter (1 to 1000 atm), s'j and s^ are much larger than Sg, s’g, and s^. 
Hence, H(H^ and H ) and Hg(H| and Hg) are treated as major species and e“, H”, and 
Hg as minor species. Thus, dropping the last term in equation (8) gives 


S’ « S’j + S^ 


(26) 


which is Dalton’s law for ideal gases. Combining equations (23) and (26) gives 



1 


2K pl 


(27) 


The value of may then be calculated from equation (26). 

To find slj, Sg, and s^, the values of s'j and from equations (27) and (26) are 
used together with appropriate equilibrium constants and the condition for electrical 
neutrality. For the reaction (3/2 )Hg — Hg + e“ 


S 7 S 3 


= K_„ 


^RT 5 ^ 2 q^ ^lT 5 / 2 q 3 ^ ^RT 5 / 2 q^ ^ 


-hc[(3w|/2) + l{-wJ/kT 


(28) 




Low Temperature {300 to 1300 K) 

In this range the medium -temperature analysis was inapplicable because it caused 

major underflows on the computer. All ionization and electronic excitation were negli- 

*+■ 

gible. and the number density Nt was much less than Nt for pressures from 

c O ^ ^ 

1. 01325x10° to 1, 01325x10 newtons per square meter (1 to 1000 atm). The equilibrium 
components of the gas were calculated from the equilibrium constant K pi given by 
equation (25). 

Assu mptions and r est ricti ons. - The assumptions and restrictions utilized in the 
low -temperature analysis are as follows: 

(1) Only the species and h| were included. 

(2) The gas was treated as an ideal gas. 

(3) Nuclear spin was taken into account in the same way as in the high - temperature 
analysis. 

Solution. - The molecule Hs was treated as the only major species, and the atom H* 

— — — — " A + * 

was treated as the minor species. Hence, N| « N, and from equation (23) and the 
perfect-gas law 
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( 32 ) 



Computer Program 

The high-, medium-, and low -temperature analyses were programmed in FORTRAN 
IV in a single program given in appendix B. 


RESULTS AND DISCUSSION 

In this section the rotational term values for h| are given. Compositions obtained 
using these term values are then presented. The effects on the composition of including 
or omitting and coulomb interactions between free charged particles are given. The 
accuracy of the compositions is discussed, and the compositions are compared with those 
obtained by other investigators. 


Rotational Term Values for 

The calculated rotational term values f for are given in table IV (missing values 

z - 1 

indicate nonexistent states). These values agree to within 700 meters with rotational 
term values calculated from wavelengths of lines measured by Herzberg and Howe 
(ref. 22). The values of f in table IV predict wave numbers of lines within 170 meters"^ 
of the values measured by Stoicheff (ref. 23), Fink, Wiggins, and Rank (ref. 24), and 
Herzberg (ref. 25). Table IV is also in excellent agreement with the calculations from 
first principles of Waech and Bernstein (ref. 26), 


Components of Hydrogen Plasma 

The number densities of species in a hydrogen plasma are given in tables V(a) to (d). 
The range of temperatures and pressures included in these tables is indicated by the area 
without crosshatching in figure 1. Typical cases in table V were checked to make sure 
they satisfied four different dissociative equilibrium equations and a Saha equation with 
lowering of the ionization potential (ref. 12). Number densities from tables V(a) and (d) 
are plotted in figures 2(a) and (b), respectively. By comparing figures 2(a) and (b), it is 
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evident that higher pressure tends to inhibit dissociation and ionization, as expected. 

Effect of H 3 on Composition, Electrical Conductivity, and Opacity 

Of the partition functions of the seven species considered at high temperatures, the 
partition function of Hg is least accurately known {ref. 13). Consequently, the effect of 
including or neglecting Hg on the number densities of e~ and H” was investigated. For 
small degrees of ionization, including Hi substantially increased the number densities of 
e' and H”, as shown in figure 3 for a pressure of 1. 01325x10 newtons per square meter 
(1000 atm). However, number densities of e* and H~ are less dependent on the inclusion 
or neglect of Hg at lower pressures. Including Hg will thus increase the calculated elec- 
trical conductivity (through e - ) and opacity (through bremsstrahlung and H adsorption). 


Effect of Coulomb Interactions Between Free Particles 

Coulomb interactions between free particles reduce the pressure as shown by equa- 
tion (8). The effect can be evaluated by calculating the coulomb compressibility 


Z S B (33) 

X - 

kT £ N* 
i=l 1 

which is 1 for an ideal gas and less than 1 if coulomb interactions between free particles 
are considered. For the conditions in this report Z c was never less than 0.92. 

A much larger effect of coulomb interactions occurs with the number densities of 
certain species. The ratio of calculated number densities of H, e , and H" with and 
without coulomb interactions between free particles is shown in figure 4 for a pressure of 
1. 01325X10 8 newtons per square meter (1000 atm). It is apparent that at this pressure 
coulomb interactions between free charged particles must be included in composition 
calculation if acceptable accuracy is to be achieved. However, at lower pressure the 
effect is smaller. 


Accuracy of Composition 

Based on the assumptions made and the constants used in the three analyses, the 
nonzero number densities in tables V(a) to (d) are accurate to four significant figures 
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with two possible exceptions. At temperatures below about 2500 K the number densities 

* * 

cf and H may be too small because of a feature in the computer program designed to 
prevent major underflows. However 9 at these low temperatures these excited species 
are negligible anyhow. In any case, this possible discrepancy is less important than 
possible inaccuracies in the assumptions made and the constants used, which are dis- 
cussed in the following paragraphs. 

For consistency it was necessary to neglect all intermolecular forces involving 
neutral species* This caused serious inaccuracies in number densities at the highest 
densities, which occurred at low temperatures where only H and Hg were significant. To 
evaluate these inaccuracies, number densities of H and Hg were calculated from the 
relations and second and third virial coefficients of Fisher (ref. 27). In all cases the 

A A 

relative effect of the virial coefficients was larger on Nj than on N^. The inaccuracy 
in N. when virial coefficients were neglected is shown in figure 1. As a result, tem- 
peratures and pressures where the inaccuracy in exceeded 20 percent were excluded 
from this report (see fig. 1). 

The only assumptions made that would be likely to cause inaccuracies exceeding 

20 percent in the unshaded area in figure 1 pertained to electronically excited states of 

H 2 ? H, Hgj and H". The number density of H 2 is possibly only accurate to within 

50 percent because (1) m eta stable rotational states were neglected, (2) shifts in energy 

levels were neglected, (3) doubly excited states were neglected, and (4) the Bet he cutoff 

* 

is somewhat arbitrary. However, it is doubtful if contributes appreciably to the 
absorption coefficient of a hydrogen plasma, and Hg is certainly unimportant thermo- 
dynamically. The number density of H is possibly only accurate to within 50 percent at 
40 000 K because (1) shifts in energy levels were neglected and (2) the Bethe cutoff is 
somewhat arbitrary* The number density of H should be much more accurate at appre- 
ciably lower temperatures. At 40 000 K there are probably appreciable inaccuracies in 
Hg ? and H" number densities because of the neglect of electronically excited states. 
The number densities of these species are negligible at this temperature, so the inac- 
curacies are academic. At appreciably lower temperatures the excited electronic states 
are unimportant. 

There is a limit to the charged-particle density at which the Debye- Hu ckel theory is 
accurate. For a hydrogen plasma, the theory is believed valid provided (refs. 28 and 12) 


r 


c 



K 


2 


N 5 + 


N c 



< 1 


(34) 


However, as r approaches 1, the accuracy can be expected to decrease (ref. 29). For 
the conditions in this report, r Q never exceeded 0, 6. 
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The largest possible sources of error in this report were the constants and potential for Hg, 
which were used to calculate the partition function {ref. 13). These quantities were based on 
Conroy's calculations from first principles (ref. 30), which are believed to be reliable. However, 
there are no experimental vibrational frequencies or dissociation energy for Hg. 


Comparison of Composition with Other Investigators 

Calculated number densities for Hg and e" obtained in this report and by two other investi- 
gators are shown in figure 5. The number density of Hg obtained from Krascella and Patch's 
results (ref. 1) was low at high temperatures because K pl was extrapolated from lower temper- 
atures. The number density of Hg obtained from Rosenbaum and Levitt's results (ref. 2) was low 
at high temperature because metastable rotational states were neglected. The number density of 
e" obtained from Krascella and Patch's results (ref, 1) was high because the lowering of the ion- 
isation potential was calculated from Ecker and Weizel’s analysis (ref. 31), which was incorrect 
(ref. 32). The number density Ng of e" obtained from Rosenbaum and Levitt's results was low 

because coulomb interactions between free charged particles were neglected. At the pressure for 
. * * 
figure 5, Hg does not affect N g appreciably for the range of Ng shown therein. 

CONCLUDING REMARKS 

The number densities of the components of a hydrogen plasma including minor species were 
calculated. Coulomb interactions between free charged particles were considered for temper- 
atures above 2000 K, Below this temperature they are insignificant. The hydrogen triatomic 
molecular ion (Hg), which had been neglected in previous calculations , was included and was 
found to increase the number densities of e" and R~ substantially for small degrees of ionization. 
The largest possible source of error in this report was the constants for Hg, which have not been 
experimentally verified. * Because of the inclusion of minor species, the results of this report 
will be particularly useful for calculating optical absorption coefficients. 

Lewis Research Center, 

National Aeronautics and Space Administration, 

Cleveland, Ohio, July 9, 1968, 

122-28-02-17-22. 


1 Since this report was written, Leventhal and Friedman (J. Chem. Phys, , vol. 49, no. 4, 
Aug. 15, 1968, pp. 1974-1975) have experimentally determined the dissociation energy of Dg. 
Allowing for differences in vibrational zero point energies of Hg and Dg and also Hg and Dg, their 
results agree within the experimental error with the dissociation energy of Hg used in refer- 
ence 13 and in this report. 
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APPENDIX A 


SYMBOLS 


A 


a n’ a n-l ,a n+l 

a 

o 

C 


c 

D 

9 


9 i (i=01, 02, ... ,05) 
d 


e 

F 

n 

f ? f(v, J) t f(v, K) 
G 

G v 

g 

H 

h 

I n' *n- V J n+1 


AI 

J 

K 


^1' ^p5' K p7 
k 


Helmholtz free energy, J 
radius of electron orbit, m 
first Bohr radius, m 

excess of negative charges over positive charges, moles 

speed of light, m/sec 

Morse well depth, (see fig. 8) 

correction column vector 

elements of 9 

* i 

correction to Dieke's {ref. 33) term values for triplet Hg, m 
charge of electron, C 
cutoff factor 

rotational term values, m -1 (see fig. 8) 

Gibbs free energy, J 

vibrational term value, (see fig. 8) 

reduced Gibbs free energy, moles 

electronic degeneracy 

hump height, m~* (see figs. 8 and 9) 

Plancks constant, J sec 

clamped -nuclei ideal -gas ionization potential, (see fig. 8) 

ideal -gas ionization potential (for rotationless ground vibrational 
states in the case of molecules), m (see fig. 8) 

absolute value of change in ionization potential, m -1 

total angular momentum quantum number 

quantum number for total angular momentum apart from electron 
spin 

equilibrium constants 
Boltzmann constant, j/K 
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m. (i=l,2, ... ,7) 

a 

N 

A 

N i (i=l, 2, ... ,7) 

N ln 

N 4 (i=l,2 7) 

N o 

n 

n + 

n % (1=1,2,.. .,7) 

Alij (1=1,2, ... ,7) 

A In nj (1=1,2 , ... ,1) 
P 

Po’ P o 

P t (1=1)2, — ,7) 
q: (1=1,2,..., 7) 



q i (1=1,2,. . 

R 


• ! 7) 


r 


r 


c 


o 

s’ 


B[ (1=1,2, 


7) 


mass of species i, kg 
total number density, 1/m 

number density of species i, l/m^ 

number density of H with principal quantum number n, l/m"* 

number of particles of species i 

Avogadr o' s numb er , 1/ mole 

principal quantum number 

effective quantum number 

moles of species i, moles 

change in n., moles 

change in natural logarithm of m 

pressure, N/m 

O 

reference pressures, N /m* 

o 

partial pressure of species i, N/m 

internal partition function of species i relative to ground elec- 
tronic state (rotationless ground vibrational state in case of 
molecules and molecular ions) 

contribution of ground electronic state to q’^ 

ideal -gas internal partition functions of species i relative to 
internal energy of e' and 

universal gas constant, J/(mole)(K) 

vibrational quantum correction 

Rydberg constant, m -1 

internuclear distance, m 

ratio of equivalent concentration to critical equivalent concen- 
tration 

equilibrium internuclear distance, m 

half the average distance between neighboring molecules, m 
nondimensionalized pressure, p/p^ 
nondimensionalized partial pressure, p./p^ 
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s 

T 

T 

T + 


T 


e 


U(r) 


u 

V 

AV 

A In V 

v 

w 


w o (e),w I (0),w n (e) 

Wj. (i = 1, 2, 5,6,7) 
X 

x- 1 

Y 

n 

Z 


z. (i=l,2, ... ,7) 
Ofj (i=l, 2, . . . , 7) 
j3 


5 


nondimensionalized pressure, p/p o 
temperature, K 

* 

energy difference between bottom of potential well of state of H« 

+ i & 

and rotationless ground vibrational state of Hg, m (see fig. 8) 

energy difference between rotationless ground vibrational states 
of Hg and H|, m" 1 (see fig. 8) 

electronic term value for H* relative to m - * 

potential energy measured from rotationless ground vibrational 
state of h|, m"* 

dimensionless fundamental vibration frequency 

<a 

volume of system, m 
3 

change in V, m 

change in natural logarithm of V 
vibrational quantum number 

energy difference between top of potential well and rotationless 
ground vibrational state, m"* (see fig. 8 and table HI) 

integrals called W{0), W(l), and W(2), respectively, in ref. 18 

9 

multipliers for An. in eq. (B14), 1/mole 
matrix for iterative corrections, 1/mole 
inverse of matrix for iterative corrections, moles 
initial column vector, 1/mole 
first element in Y, 1/mole 

number of elementary charges ion (or atom) possesses after 
ionization 

coulomb compressibility 

net number of elementary charges on species i (1, 0, or -1) 

3 3/2 

mass parameter for species i, (m )(K ' )/mole 
Morse parameter, m -1 

contribution of a state of Hg to (no cutoff) 

term value difference between dissociation asymptotes, m~* 

(see fig. 9) 
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^v 

V 

V’ 


K 

Aj (1=1,2,... ,7) 
M 

(i=l,2, . . . , 7) 

P 

a 


(j) 

U>' 

Subscripts: 

DH 

d 

ex 

id 

0 

1 

2 

3 

4 

5 

6 
7 


electric permittivity of free space, C^/N-m^ 
contribution of v states of h| to q^ 
number of gram atoms of protons, moles 
number of protons 
reduced temperature 
reciprocal Debye length, m 

o 

characteristic volume for translation for species 1, m 
reduced mass of two atoms composing a diatomic molecule, Is® 
chemical potential of species i, J 
density, kg/m a 
symmetry number 

contribution of H with principal quantum number n to before 
employing cutoff 

contribution of states with vibrational quantum number v to 
Morse oscillator vibrational constant, m~* 
fundamental vibrational frequency , m“* (see fig. 8) 

with coulomb interactions between free charged particles 
desired value 

excess due to coulomb interactions between free charged particles 
ideal gas 

previous iteration 
hydrogen atom, H 
proton, H + 
free electron, e~ 
hydrogen molecule, Hg 
negative hydrogen ion, H" 

hydrogen diatomic molecular ion, Hg 
hydrogen triatomic molecular ion, Hg 
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Superscripts: 

* electronically excited 

t electronic ground state 

Unless otherwise indicated, partial derivatives are taken with respect to one of the 
eight variables n^, n 2 , n^, n^, ng, n^, or v, and the other seven variables and T 
are understood to be held constant. The moles used in this report are gram moles. 
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Illlii 


APPENDIX B 

COMPUTER PROGRAM FOR COMPONENTS OF HYDROGEN PLASMA 
INCLUDING MINOR SPECIES 

The computer program described in this appendix computes number densities and 
other properties of a hydrogen plasma either with or without coulomb interactions between 
free charged particles. The assumptions and restrictions have already been given in the 
main text. 


Analyses 

Three analyses were used. The low -temperature analysis for temperatures from 
300 to 1300 K has already been given, as has the medium -temperature analysis for tem- 
peratures from 1300 to 2000 K. Neither of these analyses required iteration. 

The first part of the high -temperature analysis for temperatures from 2000 to 
40 000 K has also been given (eqs. (1) to (21)). In this section these equations are put in 
a convenient form for iteration, and the iteration procedure is given. 

Substituting equation (1) into equation (6), dividing by RT, and converting from num- 
bers of particles to moles give 



Substituting equation (1) into equation (8), dividing by p o , and converting from numbers 
of particles to moles give 



The desired value of s is p^/ P Q , which is designated s^. Converting equation (10) 
from numbers of particles to moles yields 
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j] sr = Dj + n 2 + 2n^ + n 5 +■ 2n g + 3n^ 


<BS) 


The desired value of tj can be any constant and is designated rj d> From equation (9) 


n 3 + n g s n 2 + n g + (B4) 


During the iteration, this equality may not be met exactly, so C is defined as 


C = n 3 + n 5 - n 2 - n 6 - n ? 


(B5) 


The desired value of C is, of course, 0, 

To find a minimum in g (eq. (Bl)) subject to three constraints (eqs. (B2), (B3), 
and (B5)), all four equations are expanded in Taylor series. 


s=s ° + Zfe) 0 ini+ (#)o 4 vt iZ($) <ani,2+ j©) 0 

i=l 0 i=l X l7 0 0 


i&vy 


6 7 


Z&) 0 iniiv+ ZZfe) 0 in ‘^ 

i=l 0 1-1 j=i+l 0 


In taking partial derivatives of g, q^ and were assumed to be constants. 

7 


s d =s 


o+ ZW 0 ^M„ 

i=l u 


AV 


(B6) 


(B7) 


= ^0 + 


Zfe) 0 Ani 

i=l u 


(B8) 
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(B9) 


o = c 0 + 


Zfeb 

i=l 0 


To simplify the nomenclature, the subscripts 0 are dropped. To avoid the llxll matrix 
entailed by the use of three Lagrange multipliers to find a minimum in g subject to the 
constraints, equations (B7) to (B9) are solved for An^, An^, and AV. 


An^ = -C - An^ + Ang + Ang + An^ 


(BIO) 


An 4 = “^d'^-7 An l““ 


1 3 

An„ — An K - An ft — An„ 

2 u 2 a 2' J 2° b 2‘ 


(BIX) 


AV = — is . ~ s - An-, - An„ - -ZiL (An„ + An„ + An 7 - C - An,-) 
as 1 d 3n 1 1 3n 2 2 3n 3 2 6 7 5 


9V 


as 

an. 


[j (»d - t|) - i An, - i in 2 - i 4n 6 - An,. - 1 an,] 


3s A „ 3s Ari 3s 
3n 5 0 3n g 0 3n 7 ' 1 


(B12) 


Equations (BIO) to (B12) are then substituted into equation (B6). The conditions for a 
minimum in g are then 


JS -=0 (i = 1, 2, 5, 6, 7) (B13) 

aAnj 

For an i of 1, equation (B13) takes the form 

w-^ An^ + Wg An 2 + Wg An g + Wg Ang + w ? An^ - yj = 0 (B14) 

where w^, Wg, Wg, Wg, w^, and are complicated functions independent of An^, Ang, 
An 5 , Ang, and An^. Equation (B 14) is approximately 
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Wjn^ A In + WgHg A In n 2 + w g n 5 A In n g + w g n g A In n g + A In n^ - = 0 (B15) 

Similar equations follow for i = 2, 5, 6, and 7, These five equations can be written in 
matrix notation as 


X9 - Y = 0 

where X is a 5x5 matrix, Y is a five-element column vector, and 


^01 


A In 

^02 


A In XI2 

CO 

o 

a 

= 9 = 

A In 

®04 


A In ng 

®05 

1 

A In ritj 


(B16) 


(B17) 


To find a minimum in g, initial estimates are provided for the seven m and V, 
and i? is found from equation (B17). From equation (BIO) 

-C - nc A In njr + n 0 A In n« + n c A In n c + n„ A in n„ 

A In n 3 = “ — — _ 1 l (B18) 

n 3 


From equation (Bll) 


Tjj - ?? - n, A Inn, - n 9 A In n 
A In n 4 = — — - 1 t .* — . 


n,- A in n, - 2n c A In n„ - 3n„ A in n„ 

5 5 6 6 7 7 

2n 4 


(B19) 


From equation (B12) and the relations 3s/ 5n^ = 3s/ and ds/dn^ = ds/dn^ = 3s/3n g 
= 8s/3n g = 3s/ 3n^ 


A In V = . 


3s 

0V‘ 


s - ~ ii-n, A In n, +(I -2^5. 

2 3n, 1 \2 3n, 3n, 


;)■ 


2 A in n 2 + 


1 _3s^ 

2 3n, 


n 5 A ln n 5 


+ 




A ln n g 


- -2^-^n 7 A ln n„ + C — 

v 2 3 Hj 3n 2 y * ' 0n 2 


1 3s_ 

2 3nj 


(h d - 1?) 


(B20) 
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The logarithms of the seven and V are corrected by means of equations (B17) to 
(B20), and the next iteration may be begun. 

Program Procedure 

The program is in FORTRAN IV, IBM Version 13, for an IBM 709411/7044 Direct - 
Coupled System with 32 000- word storage. The names of the FORTRAN variables and 
their corresponding symbols, equations, or comments are as follows 


FORTRAN 

Mathematical 

Comments or equation 

symbol 

symbol 


A(5,U) 

X" 1 

X" 1 is contained in the square subarray A(l, 6) to 



A(5, 10) 

AKAP 

K 

reciprocal Debye length 

AKP1 

K P1 

equilibrium constant 

AKP2 


s 6^ s l s 2 ^ or *deal gas 

AKP4 


s 2 s 3 /s l tor Weal 

AKP5 

^5 

equilibrium constant 

AKP7H 

K V2 

“p7 

ALGM{7) 


In n^ 

ALPHA (7) 


mass parameter for species i 

ALV 


In V 

ARBI 


arbitrary multiplier for r Q (always 1) 

ARB 2 


arbitrary multiplier for Al (always 1) 

BETA (6, 23) 


Morse parameter (see table II for subscripts) 

C 

C 

excess of negative charges over positive charges 

CN(7) 

VS 

N i 

number density for species i 

•“5 

CNT 


^ X 


X 

1=1 

CN1(17) 

N ln 

number density of H with principal quantum num- 


ber n 
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I 


FORTRAN 

Mathematical 

Comments or equation 

symbol 

symbol 

17 * 

CN1S 


nS” 1 " 

CN4G 

N* 

w 4 

number density of h| 

CN4H(6 , 23) 


number density of a state of (see table n for 

subscripts) 

CN4S 


number density of H 2 (all states) 

CN4V(15) 


number density of h| with vibrational quantum num- 
ber V 

CN6V(20) 


number density of with vibrational quantum num- 

ber V 

CUT1 


(rya// 2 (approximate principal quantum number 
at Be the cutoff) 

CUT2 


1/2 

(Ry/Al) (approximate principal quantum number 

at Debye -Huckel cutoff) 

D(5) 

9 

correction column vector 

DH 


0 for ideal gases; 1 for Debye-Hiickel plasma 
(normally 1) 

DELI 

Al 

absolute value of lowering of ionization potential 

DELIN 


new value of AI 

DELLGM(7) 


A in 

DELLV 


A In V 

ED(6, 23) 


hcD/k (see table II for subscripts) 

EEXP 


e -u/2 

EF(6 f 23) 


hccu’/k (see table II for subscripts) 

EH(I7) 


hcT e /k for H with principal quantum number n 

EH2(6,23) 


hcT/k (see table n for subscripts) 

ERETA 


v d ~v 

ERS 


s d - s 

ETA 

n 

actual gram -atoms of protons 
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FORTRAN Mathematical Comments or equation 

symbol symbol 


ETAD 

^d 

desired gram -atoms of protons (always 3) 

FH2(15, 39) 

f 

rotational term values for Hj 

FH2P(20,42) 

f 

rotational term values for 

FN(7) 


dg/dn i 

FNN(7, 7) 


d^g/da^dn^ 

FNV(7) 


0 2 g/avan. 

FV 


dg/dv 

FVV 


a 2 g/av 2 

GAM(6, 23) 

^h 

contribution of a state of Hg to (no cutoff) (see 

table II for subscripts) 

GAMC(6, 23) 


Fn^h (see table H for subscripts) 

GAMCIM 


with F calculated for Be the cutoff 
n n n 

GH2(15) 

G v 

t 

vibrational term value for 

GH2P{20) 

G v 

vibrational term value for 

GMO(7) 

n i 

number of moles of species i 

IFLAG 


0 for same DH, and IPR for next case; 1 for 

read new DH - - IPR card for next case 

IMAX(23) 


* 

maximum value of first subscript of Hg for value of 
second subscript given in parentheses (see table II 


for subscripts) 

IPR 0 for long output; 1 for short output 

11 number of iterations for cutoff to converge in high- 

temperature iteration {999 if no convergence or if 
medium- or low-temperature composition calcu- 
lated) 

12 number of iterations for n. and V to converge in 

high- temperature iteration (999 if no convergence; 
998 if medium- or low- temperature composition 
calculated) 

19 index for high-temperature iterations 
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FORTRAN Mathematical Comments or equation 

symbol symbol 


LIM 


limit on 19 (always 15) 


OMEGA(20) 

fi v 

contribution of states with vibrational quantum num- 
ber v to q£ 

PD 

Pd 

desired pressure 


PHI(17) 


contribution of states with principal quantum number 
n to qj (no cutoff) 

PHIC(17) 


first member is 2, others are F n <? n 

PHICN 


F n 9 ? n with F n calculated for Bethe cutoff 

PR 


P^V/RT 


PROD(7) 

PRODD 


N o A i 

3 

/ 2 \ 2 

/ e N \ i 7 9 

0 1 T z 2 n 

1/2 


16ttN 0 

\e Q kT f Vi-1 1 



PRODDD 


PO 

Po 

POP 

Po 

Q(7) 

<k 

QUOT 


Q4P 

*4 

Q6P 

q^ 

Q7P 

q 7 



reference pressure for high -temperature iteration 
(always 10 6 N/m 2 ) 

reference pressure for equilibrium constants 
(always 1.01325X10 5 N/m 2 ) 

partition function for species i referred to internal 
energy of atoms and free electrons 

t 3/2 

partition function of referred to rotationless 
ground vibrational state of h| 

partition function of Hg referred to rotationless 
ground vibrational state 

partition function of referred to rotationless 
ground vibrational state 
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FORTRAN 

symbol 


Mathematical 

symbol 


Comments or equation 


Q7PT(67) 

R 

RAT 

RATI 

RAT2 

RE(6, 23) 

RE L{8) 

RES 

RHO 

RO 

RON 

ROOT 
ROOT1 
ROOT 2 

ROOTA(18J 

ROOTA(7 > 23) 

ROOTR 


table of having one-to-one correspondence with 
TT table 

ratio of equivalent concentration to critical equi- 
valent concentration 

vibrational quantum correction 

absolute value of relative change in r Q in one 
iteration 


absolute value of relative change in AI in one 
iteration 

equilibrium internuclear distance (see table n for 
subscripts) 

absolute value of relative change in n- and V in 
one iteration 



density 

half the average distance between neighboring 
molecules 


new value for r r 


10" 16 e 2 N 7 


e Q kTV 


- £ Vi 


negative of average of two adjacent ROOT A in 
Rydberg series 


negative of average of two adjacent ROOT A in 
Rydberg series 

a 1 / 2 for H 
n 

a n^ * or ( see ta ^ le 11 for subscripts) 
(r Q ARBl ) l / 2 
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FORTRAN 

symbol 


Mathematical 

symbol 


Comments or equation 


s 


P/P 0 

SD 


Pd/Po 

SIGMA 1 


£ (2J + l)exp JlSH for h| 

SIGMA2 


2 (2J + l)exp for h| 

even J \ kT / 

SIGMA 3 


2 (2K + l)exp/- for Hj 

odd K V kT/ 

SIGMA4 


Y (2K + l)exp(-M\ f ° r Hj 

e^nK 

SM 


l 2 

E z.n. 

i=l 1 1 

SN(2) 


3s/ 3n. where i is 1 or 2 

SNlSV 


(da/dn^/ids/dV) 

SN2SV 


(3s/3n 2 )/(3s/3V) 

SP(7) 


p/pi 

SPD 


Pd /p i 

SP3 


< K P 4Pi/py 1/2 

STAR1 


average of two adjacent I in Rydberg 

STAR2 


average of two adjacent I in Rydberg 

STARI(7, 23) 

I 

* 

for H 2 (see table II for subscripts) 

STARI(18) 

I 

for H with principal quantum number 

SUM 

“h| 

contribution of h| to 

SUMM 


7 

E n, 

i=l 1 

SUM1 


17 

y F <p 
n=2 


series 

series 


n 
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FORTRAN 

symbol 

SUM4 

SUM9 

SV 

T 

THETA 

THETAT(29) 

TOL1 

TOL2 

TS1 

TS2 

TT(67) 

U 

V 

w 

WEIGHT(6, 23) 
WO 

WOT(29) 

W1 


Mathematical 

symbol 


Comments or equation 


T 

e 


u 

V 


g 


e 


W 


o 


W 


I 


£ Vh 

^2 + ^3 + ^5 + ^6 + ^7 


8s/ 3 V 

temperature 

reduced temperature 

table of $ having one-to-one correspondence with 
WOT, WIT, and W2T tables 

maximum allowable absolute value of relative change 
in r Q and AI in last iteration before freezing 
cutoffs (always 0.001) 

maximum allowable absolute value of relative change 
in n^ and V for last iteration (always 0.00001) 

dividing point of low and medium temperatures 

dividing point of medium and high temperatures 

table of T having one-to-one correspondence with 
Q7PT table 

dimensionless fundamental vibration frequency 


volume of system 


S + c i5- - - — (v d - v) 

8n 2 2 8nj ° 

8S 

8V 


electronic statistical weight (see table II for sub- 
scripts) 

integral called W(0) in ref. 18 

table of W Q having one-to-one correspondence 
with 6 table 


integral called W(l) in ref. 18 
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FORTRAN Mathematical Comments or equation 

symbol symbol 


W1T(29) 

W2 

W2T(29) 

X(5, 5) 
Y{5) 

Z(7) 

ZC 

ZETA{15) 


table of Wj having one-to-one correspondence with 

8 table 

Wjj integral called W(2) in ref. 18 

table of Wjj having one-to-one correspondence with 

9 table 

X matrix for iterative corrections 

Y initial column vector 

»? 

Z c coulomb compressibility 

t 

contribution of v states of to q^ 


(Units are given in appendix A. ) 

A simplified flow diagram of the program is shown in figure 6. The input quantity 
BH determines whether or not coulomb interactions between free charged particles are 
included. Other features of the flow diagram are discussed in the following paragraphs. 

Two sets of criteria were used for convergence because and were assumed 
to be constants when taking partial derivatives of g, Actually, and q^ depend on the 
n. because of the cutoffs, but this dependence must be neglected when taking partial deri- 
vatives if the solutions are to obey Saha's equations (ref. 34) as is customary. This 
results in an inconsistency. Problems due to the inconsistency were averted by using 
two sets of criteria for convergence , as shown in figure 6. The first set was that relative 
changes in one iteration in both r Q and A1 must be less than TOL1 (set at 0. 001), 

After these criteria were met, r Q and AI and, hence, the cutoffs were not changed 
again. The second set of criteria was that relative changes in the seven and V be 
less than TOL2 (set at 0.00001), 

There are two options for output. If the input quantity IPR is 0, the long output (1 
page per case) is used. If IPR is not 0 and the iteration (in the case of high temperature) 
converges, the short output (1 line per case) is used. 

A block diagram of the program and the eight subprograms is given in figure 7. 


Sample Problems 

An input form with data for seven sample cases is given in table VI. Note that 
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IFLAG is 1 on the second card because it is desired to change IPR by means of the next 
card. The output is shown in table VII. Case 1 had long output, whereas cases 2 to 7 had 
short output. The order of the quantities in short output is T and number densities for 
h|, H*, H*, H*, e", H", H+, H+, H + , H 2 , and H. 

Running Characteristics 

A control card is required so that minor underflows are not printed out. With short 
output 168 cases required I. 56 minutes to execute. With long output 45 cases required 
0.71 minute to execute. For all high-temperature problems, II was between 2 and 4, 
whereas 12 was between 3 and 8. 
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Listing of Program 


ilBFTC MAIN DEBUG 

C PROGRAM FOR COMPONENTS OF A HYDROGEN PLASMA INCLUDING MINOR SPECIES 
C 

DIMENSION IMAX ( 23 ) » GAM ( 6,23 ) *GAMC (6,23) ,PHI(I7)»PHIC(l7) t 0t7), 
1PRQD<7> ,ALPHA(7 >,EH( 17 ) * EF I 6, 23 ) , EDI 6, 23 > , EH2 I 6,23 ) , WE 1GHT < 6, 23 I , 
2PE(6,23) , BETA 1 6, 23) , JP1MI15) ,KP1M(20) ,CN<7) ,GH2( 15 > ,GH2P ( 20 > , 

3FH2I 15,39) ,FH2P< 20,42) ,Q7PT < 67 ) , TT < 67 > , THET AT < 29 ) , WOT < 29 ) , WIT < 29 ) , 
4W2T (29) » ZETA ( 1 5 ) ,Z m»CNl(17) ,CN4H(6,23),CN4V< 15) ,CN6V<20) , 

5UMEGA (20),REL(8),SP(7) 

COMMON/ BLOCK 1 /OH ,T ,PO t PROD, ETAD,SD/BL0CK2/PHI »PHIC » SUM! /BLOCK 3/ SUM 
l ,GAM,GAMC t SUM4/BL0CK4/ I MAX/ BLOCK 5/ ARBI , ARB2/BL0CK6/2 
2/ BLOCK 7/ ALPHA ,EH, EF ,E0, EH2 ,WE IGHT, RE, BET A, JPlM,KPlM t GH2,GH2P,FH2/ 
3BL0CK8/FH2P, 07PT ,TT,TH£TAT,W0T,W1T ,W2T 
DOUBLE PRECISION V,GM0(7) ,X(5,5) , Y ( 5 ) , S,ETA,C,SM,R00T ,SUMM,A( 5,11) 
1 , S N 1 2 ) » S V » D ( 5 ) ,DELLGM(7) , D F L L V , A L GM ( 7 ) ,ALV 

2 READ! 5,200) DH,PD,IPR 
WRITE(6,201) OH , PD, IPR 

C ASSIGN VALUES TO CONSTANTS. 

TS1 =1 300. 

TS2=2Q00. 

P0=1 .E+6 
P0P=1.01325E+5 
TOL 1= .001 
T0L2=. 00001 
ARBl=l. 

ARB2=1. 

L I M= 1 5 
ETAD=3. 

K = 0 

1 READ! 5 ,202 ) T , 1 FLAG 

C END OF INPUT. FIND CONTRIBUTION OF GROUND ELECTRONIC STATE OF H2 TO H2 
C PART IT ICN FUNCT ION. 

S0M=0. 

DO 3 KVP1=1,15 
S I GMA1=0 . 

MAX=JP1M(KVP1 > 

DO 4 JP 1 =2, MAX , 2 

4 SIGMA1=SIGMA1+FLUAT<2*JP1-1)*EXPC(— .0143879*FH2<KVPl, JP1)/T) 

S I GMA2=0. 

00 5 JJP1=1,MAX,2 

5 S I GMA2=S l GMA2+FL0AT <2*JJP1-1 )#EXPC (—. 0143879# FH2(KVP1,JJPI) /T) 
ZETAIKVP1 )=.25*EXPC ( - .0143879*GH2 < KVP 1 )/T)*(3.*SI GMA1+S I GMA2 ) 

DEBUG SIGMA l, SI GMA2 »ZETA(KVPl ) 

3 SUM=SUM+Z6TAIKVP1) 

C INITIALIZE VARIABLES. 

R0=1.E-8*(3.2958*T/PD) **.33333 
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DEBUG RO, SUM 

RONaRO 

DEL 1=0. 

DO 46 N I = 1 1 7 
46 DELLGMtNI )=0. 

DELL V=0. 

ERS=0. 

C=0. 

ERETA=0 » 

06P=0 . 

11=999 

12=998 

C TEST IF T IS LOW. 

IF< T-TS1 )7,7,8 

C FIND LOW T COMPOSITION FROM HERE UP TO STATEMENT 8. 

7 0(1) =2. 

CN1 ( 1 )=SORT(PD*SORT(T) /( 1 .0775E-3*SUM ) ) *0 < 1 >*( 1 . £+23*EXP ( -25980.0 
1/T) 1/1.3805 
CN1 S=0 . 

CN( 2 ) =0. 

CN( 3 ) =0 . 

CN ( 4 ) = PD/ ( 1 .3805E-2 3*T ) 

CN4S=0 • 

CN ( 5 ) =0. 

CN ( 6 ) =0 ♦ 

CN ( 7 ) =0. 

CN4G=CN ( 4 ) 

CN ( 1 ) =CN1 ( 1 ) 

Q4P=SUM 

07P=0. 

GO TO LOO 

C FOR MEDIUM AND HIGH T FIND CONTRIBUTIONS OF H2* TO H2 PARTITION 
C FUNCTION ASSUMING NO CUTOFF. 

8 DO 6 M=l,23 
I M A = I MA X ( M) 

DO 6 1 1=2, IMA 
U= E F ( I I , M ) / T 
EEXP=EXP(-U/2. ) 

RAT=U*EEXP/< l.-EEXP**2 ) 

THETA=T/ED( I I » M ) 

V'0 = B I NT ( THETAT , WOT, THETA, 29) 

W1=BINT(THETAT, WIT, THETA, 29) 

W2=BINT (THETAT, W2T , THETA, 29) 

GAM ( I I , M ) =4 «7660*WE I GHT ( I I , M > *EXPC < -EH? ( I I ,M ) /T ) *T** 1 .5* ( l.F+26* 
IRE ( II , M ) **2*W0/BETA ( I I ,M)+2.E+26*RE < I I , M )*W1 /BET A ( II, M )**2+ l #F + 26 
2*W2/BETA( I I ,M>**3)*RAT 
DEBUG U,RAT, THETA, NO, W1 ,W2, GAM ( 1 1 ,M) , I I ,M 
6 CONTINUE 

C FIND CONTRIBUTIONS OF ELECTRONIC STATES OF H TO H PARTITION FUNCTION 
C ASSUMING NO CUTOFF. 

DO 9 N= 1 » 1 7 

9 PF I ( N ) =FLOAT ( 2*N**2 ) *EXPC ( — EH ( N ) /T ) 

DEBUG (PHI ( I ), 1=1, 17) 

C FIND PARTITION FUNCTIONS OF H AND H2 WITH CUTOFF. 

0( 1)=Q1PF(R0,DELI ) 

DEBUG (PHICU) , 1 = 1 ,17) 
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Q4P=G4PF ( RO , DEL I ) 

C FIND PARTITION FUNCTION OF H3+. 

07P=B I NT ( TT,07PT,T,67 ) 

C FIND EQUILIBRIUM CONSTANTS USED FOR MEDIUM AND HIGH T, 

AKP1»1.07756-3*Q4P*P0P*EXP( 5 1959. 9/ T )/(Q( 1 )**2*T**2.5) 
AKP5=15.010*POP*EXP I 8947. 4/ T I / ( 0< 1)*T**2.5 ) 

AKP7H=SQRT ( 1 . 42 1 1 E-3*Q7 P ) . 625*EXP < -65567. / T )/( Q4P**. 75*POP** 
1.25) 

SPD«PD/ POP 

DEBUG 0(1) ,Q4P»Q7P,AKP1,AKP5,AKP7H,SPD,( IGAMC ( I , L ) , 1=2, 6) ,L=1 ,23 > 
C TEST IF T IS MEDIUM OR HIGH. 

IF(T-T$2)10,10,13 

C FIND MEDIUM T COMPOSITION FROM HERE THROUGH STATEMENT LL. 

10 SP(U = (SORT (4.*SPD*AKPL+1. ) -1 . ) / (2.*AKP1 ) 

SP(4)=SP0-SP( 1 ) 

SP( 3) =AKP7H*SGRT ( S P ( 4 ) ** 1 . 5/ ( AK P5*SP ( l ) + !.)) 

SP ( 7) =AKP7H*SP (4 )**.75*SQRT { AKP5*SP< 1 )+i. ) 

SP(5J=SP(7)-SP(3) 

SP ( 2 ) =0 . 

SP( 6 )=Q. 

DO 11 NA = 1 , 7 

11 CN< NA J=SP<NA)=PPOP/< 1 .3805E-23*T ) 

C FIND QUANTITIES NEEDED FOR MEDIUM AND HIGH T, 

12 CN1S=SUM1*CN< L ) / 0 ( 1 ) 

CN1 ( 1 ) *CN( 1 J-CNl S 
CN4S= SUM4#CN ( 4 ) /Q4P 
CN4G=CN(4)-CN4S 

100 DEL IN=.U614E-2*SGRT ( .20998E-3*< CNI2 )+CN( 3)+CN( 5 )+CN( 6)+CN (7 ) )/T ) 
GO TO 40 

C FIND HIGH T COMPOSITION FROM HERE UP TO STATEMENT 40. 

13 S 0= PD/ PO 
12=999 

C FIND H2+ PARTITION FUNCTION. 

DO 14 KKVP1=1 , 20 
S I GMA3=0. 

MA X=K P1M(KKVP1 ) 

00 15 KP1=2,MAX,2 

15 SIGMA3=SIGMA3+FL0AT(2#KP1-1 )*EXPC(-.U143R79*hH2P ( KKVP1,KP1)/T) 

S I GM A4=0 . 

DO 16 KKP 1 = 1 , M AX , 2 

L 6 S IGMA4=S IGMA4+FL0AT (2*KKP1-1 )*EXPC{-.0143H79*FH2P(KKVP1,KKP1)/T) 
OMEGA (KK VP 1 ) =EXPC ( -,0143B79^GH2P(KKVP1 > / T ) *1 1 . 5* S I GMA3+ . 5* S I GMA4 ) 
DEBUG S I GMA3, S I GMA4 , OMEGA ( KKV P 1 ) 

14 06P = Q6 P+OMEGA t KK V P 1 ) 

C FIND PARTITION FUNCTIONS REFERRED TO H AND E-. 

0 (2 )=EXP< -1 57805. /T ) 

0< 3)=2. 

d (4) =G4P*EXP t 51 9 59. 9/ T ) 

Q( 5) =EXP (8947.4/T) 

0(6) =Q6P*EXP ( -1 27044. /T ) 

Q( 7)«Q7P*EXP( -53195./T ) 

C FIND EQUILIBRIUM CONSTANTS USED FOR INITIAL ESTIMATES FOR HIGH T. 

AKP2 = .0010779*06P*P0P*EXP( 30760.5/ T ) / f 0 1 1 )*T**2.5> 
AKP4=.066615*T**2 .5*EXP( -157805. /T ) l (Q(l )#POP ) 

DEBUG SD,06P,AKP2, AKP4, (CM I) , I =1, 7 > 

C MAKE INITIAL ESTIMATES FOR HIGH T. 
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SPm=MDSQRT< l.DO+4.DO*AKPi*{2.DO*!AKP4-DSORT!AKP4)*DSQRT<AKP4+ 
ISPDi H-SPD) ) —X ,D0 > / !2.DO*AKPl) 

SP( 4) =AKP1*$P ( 1 )**2 

SP < 3)=AKP7H*SQRT(SPC4)**1.5/( AKP5*SP(1)+1. )) 

SP3=SQRT{ AKP4*5P( 1 ) ) 

DEBUG SP!1> ,SP!3> ,SP!4> ,SP3 
IF!SP(3)-SP3) 18*18*17 

17 SPm = AKP7H#SP(4)**.75*S0RT t AKP5*SPU ) + l.) 

SP!5)=SP(7)~SPt3) 

SP!2)=SPm*AKP4/SP<3) 

GO TO 119 

18 SP(3)=SP3 
SP( 2)=SP3 

SP! 7 >=AKP7H**2*SP (4)**1 ,5/SP (3) 

SP( 5)=AKP5*SP 1 1 )*SP(3) 

119 SP ( 6)=SP ( 1>#SP( 2 )*AKP2 

V=8.3143*T*ETAD/ ( POP* < SP ( 1) + S P < 2 ) +2 .*SP ( 4 )) > 

PR=POP*V/ < 8.3143*T ) 

DO 19 NB=1 » 7 
GMO ( N8 ) =PR *S P ( N8 1 

19 ALGM(NB)=DLOG(GMO(NBl ) 

A L V=DL OG ( V ) 

SH*0. 

OUOT=T**l .5 

DO 20 JJ=1,7 

PROD ! J J ) =AL PH A { J J ) / QUOT 

20 SM=SM+2(JJ)*GM0(JJ) 

ROOT = DSORT ( 1.2646D4*SM/ ( T*V1 > 

DEBUG V* PR*ALV*SM,ROQT* ( GMO 1 I ) ,ALGM( I ), 1*1*7) 

C ENTER HIGH T ITERATION EXTENDING THROUGH STATEMENT 122. 

DO 122 19=1 jLIM 
SUMM=0 . 

DO 21 NC = 1 * 7 

21 SUMM=$UMM+GMO(NC> 

S= 8.3143*T*SUMM/ ! PO*V 1 S3 10*T*ROOT**3*DH/PO 

ET A = GMO I 1 I+GM0I2 )+2.*GM0( 4)+GM0( 5)+2 .*GMG(6)+3.*GM0 < 7) 

C = GMO! 3 > + GMO( 5 )-GMO! 2 I -GMO I 6 ) -GMO I 7 ) 

DEBUG SUMM,S,6TA,C 

C FIND MATRIX AND INITIAL VECTOR FOR ITERATIVE CORRECTIONS. 

CALL MATR ( V ,GMO* S * ETA*C , 0 , ROOT , SUMM , SM , X , Y , SN, SV ) 

C INVERT MATRIX. 

CALL INVERT I X » A ) 

DEBUG < (A! I ,L) ,L=6, 10) * 1 = 1*5) 

C FIND CORRECTION VECTOR FUR ITERATIVE CORRECTIONS. 

199 DO 23 ND=1»5 
D ( ND ) =0 . 

DU 23 NE = 1 » 5 

23 D { ND ) =D ( MD ) + A ( ND , NE + 5 ) * Y < NE > 

DELLGM! 1 )=D(1) 

DELLGM ( 2 ) = D t 2 ) 

DELLGM ( 5 ) =D{ 3 ) 

DELLGM! 6 ) =D t 4 ) 

DELLGM! 7 ) =D{ 5 ) 

DELLGM(3) = (-C-GM0(5)*DELLGM(5 ) +GMO ( 2 ) *DEL LGM ( 2 ) +GMD (6)*DELLGM(6)+ 
1 GMO { 7 ) *DELLGM < 7 ) ) /GMO! 3) 

DELLGM {4) = ! ETA D— ETA -GMO! 1 ) *DEL LGM! I > -GMO! 2 ) *DELLGM < 2 ) -GMO ( 5 ) * 
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1D£LLGM(5)-2.4GM0(6)*DELLGM(6)-3.*GMQ( 7 )*DELLGM(7) ) /< 2.*GM0( 4 ) ) 
DELL V= ( SLt— S— »5D0*SN ( 1 )*GMO 1 1 ) *DEL LGM ( 1 )+(.5D0*SN( 1 ) -2 ,*SN ( 2 ) ) * 
1GM0(2)*0ELLGM<2)+.5Q04SN( 1 )*GM0<5 )*DFLLGM( 5) + < SN ( 1 ) -2 . *SN ( 2 ) ) * 
2GMG( 6 ) ^OELLGM ( 6 ) + ( 1 . 5DO#SN( 1)-2.*SN(2 ) )*GMO( 7 ) 4DEL LGM < 7 ) +C*SN ( 2 ) 
3«5D0*SN( 1 }*( ET AD-ETA) )/{ V*SV> 

DEBUG DELL V* ( DELLGM { I >,1=1,7) 

DO 24 NF = 1 1 7 

C MAKE ITERATIVE CORRECTIONS. 

ALGMf NF ) =ALGM { NF }+DFLLGM( NF I 

24 GMO(NF)=DEXP(ALGM(NF)) 

ALV=ALV+DELLV 

V=DEXP ( AL V ) 

SM=0. 

DO 25 J J J = 1 » 7 

25 SM=SM+Z I J JJ ) *GMO( JJJ J 

ROOT = OSORT ( 1 .2646044SM/ ( T*V ) > 

DEBUG SM » ROOT * AL V * V* (GMO< I ) » ALGM ( I) , 1 = 1,7 ) 

C TEST IF RO AND DELI HAVE PREVIOUSLY CONVERGED. 

I F ( 11-999)27,26,27 
C FIND NEW VALUES FUR RO AND DELI. 

26 RON=l.E-8*{ .396404V/ (GMO( l)+GMO(2 ) +GMO ( 4 ) +GMO ( 5 ) +GMO t 6 ) +GMO [ 7 ) ) ) 
1#. 33333 

DEL I N= 1 . 1 61 4 E+ 5* ROOT 

RAT1=AKS ( ( RON-RO ) /RON J 

RAT 2 = AHS ( (DEL IN-DEL I ) / DEL I N ) 

C TEST IF RO AND DELI HAVE CONVERGED. 

IF(RAT1»LT.TQLI . AND.RAT2.LT .TDL1 ) 11=19 
C MAKE CORRECTIONS BASED UN NEW VALUES OF RO AND DELI. 

RO=RON 

DEL I =DEL IN 

0(1 )=01PF(R0»DELI ) 

DEBUG II, RATI , RAT2 ,R0 ,DEL I ,0 (1) ,{ PHIC ( I) , 1=1 , 17 ) 

Q4P=04PF ( RO , DEL I ) 

0(4) =04P*EXP ( 51959.9/T) 

OEBUG 04 P ,0(4) , ( < GAMC ( I , L ) t I =2 1 6 ) , L = 1 , 2 3 ) 

GO TO 122 

C RO AND DELI HAVING CONVERGED, FIND RELATIVE CHANGES IN MOLES AND 
C VOLUME. 

27 DO 30 NG=1 , 7 

30 REL ( NG ) =DARS ( DEXP ( DELLGM ( NG ) ) -1 . ) 

REL(8)=DABS(D6XP(DFLLV)-1. ) 

OEBUG (REL ( I ), 1=1,8) 

C TEST IF MOLES AND VOLUME HAVE CONVERGED. 

I F ( REL ( 1 ) .LT.T0L2.AND.RELI2) .LT . T0L2 . AND . R EL ( 3 ) .LT.T0L2.AND 
I .REL (4) .LT.TnL2.AND.REL ( 5 ) . LT. T0L2 . AND. R EL ( 6 ) , LT . T0L2 . AND. 

2REL( 7) .LT .T0L2.AND.REL (8) .LT . T0L2 ) GO TP 39 
122 CONTINUE 
GO TO 97 
39 12=19 
97 ERS=SD— S 

ERETA=ETAD— ETA 

C FIND HIGH T NUMBER DENSITIES FROM MOLES AND VOLUME. 

00 42 NH*1 , 7 

42 CN(NH)=6.0225E+23*(GM0(NH)/V> 

R0!M=l.E-8*( .396404V/ (GMO( X ) +GMO( 2 ) +GMO ( 4 > + GMO ( 5 ) +GMO ( 6 ) + GMO ( 7 ) ) ) 
i*. 33333 



GO TO 12 

C TEST IF SHORT OUTPUT TO BE USED. 

40 I F { IPR.NE.C.AND. 12.NE.999JG0 TO 43 
C TEST QUANTITIES PECULIAR TO LONG OUTPUT. 

CUT 1 « SORT ( RO/ 5 .29 l 7E-1 1 ) 

IF! DEL 1)89. 89.90 

89 CUT2= 1 . E+37 
GO TO 91 

90 CUT2=SQRT( 10969308. /DELI) 

91 SUM9=CN(2)+CN(3)+CN( 5 ) +CN ( 6 ) +CN ( 7 ) 

CNT*SUM9+CN( 1 )+CN(4) 

00 93 KVVP1 = 1 ,15 

93 CN4VIKVVP1 )=CN4G*ZETA(KVVP1 ) /SUM 
DO 94 KKV V P 1=1,20 
IF(T-TS2)102»102,10L 

101 CN6V (KKVVP1 ) =CN ( 6 I IOMEGA (KKVVP1)/Q6P 
GO TO 94 

102 CN6V (KKVVP1 >=0. 

94 CONTINUE 

DO 95 MMM=1,23 

1 W= I MAX ( MMM ) 

DO 95 II 1=2 , IW 
IFIT-TS1 >104,104, 103 

103 CN4H ( I I 1 ,MMM)*CN (4)*GAMC ( I I I , MMM I / G4P 
GO TO 95 

104 CN4H ( I 1 1 , MMM 1 =0 . 

95 CONTINUE 

DO 96 NK=2 ♦ 17 
IF { T-TS1 ) 106,106,105 
10 5 CN1 1 NK >=PHIC (NK >*CN ( 1 ) /Q til 
GO TO 96 
106 CN 1 { NK 1 =0 • 

96 CONTINUE 

RH0=1 .6734E— 27*CN< 1 ) + 1 . 672 5E-27#CN ( 2 ) +9 . 109 1 E-3 1*CN ( 3 ) +3 . 3469E-2 7 
1#CN< 4>+1.6743E-27*CN(5 ) +3 . 34 60E-2 7*CN < 6! +5 . 0 1 94E-2 7*CN ( 7) 

ZC=P 0/ { T*1 .3805E-234CNT ) 

AKAP=S0RT<.2O998E-3*SUM9/T) 

I F ( T-TS1 ) 151.151,150 

151 R=0. 

GO TO 152 

150 R=AKAP**2/ I 4.*(3.1416*SUM9)**.66667 ) 

C WRITE LONG OUTPUT. 

152 WRIT6(6,204)PD,T,RH0,CNT 
WRITE (6,205)ZC,AKAP,R 
WRITE (6,206 ) ( CN ( I > , I =1 , 7 > 

WRITE (6,207 ) ( CN4V (11.1*1,15) 

WRITE (6,208) CN4G, (CM4H< ! , 1 ) , 1=2,4) 

WRITE (6,209) ( CN4H ( I ,2), 1=2,4) 

WRITE! 6,210) ( CN4H ( I ,3) , 1=2 ,4 ) 

WRITE (6,21 1 ) ( CN4H ( I ,4) , 1=2,3) 

WRITE! 6,212 ) ( C N4H ( I ,5 ) , 1=2,3 ) 

WRITE ( 6 , 2 13) ( CN4H ( I ,6) , 1=2,3) 

WRITE (6, 214) (CN4H12, I ) ,1=7, 10) 

WRITE (6,215) ( CN4H ( 1 ,11 ), 1=2,4) 

WRITF(6»216) (CN4H( I ,12) ,1=2,3) 

WRITE (6,2 17) (CN4H( 1,13) ,1=2,4) 
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WRITE(6,218) (CN4H( I, 14), 1*2,3) 

WRITE (4, 2 19) ( CN4H ( 1,15) *1*2,3) 

WRITE (6,220) (CN4H ( I, 16), 1=2, 3) 

WRITE (6 ,221 ) (CW4H(2, I ), 1=17,20 ) 

WR ITE( 6,222) ( ( CN4H ( I , L ) , I =2 , 6 ) ,L=?1 , 22 > 

WRITE (6,223) (CN4H< I ,23) , 1=2,4) 

WR ITE ( 6,224} CN4S 

WRITE (6, 22 5) (CN1U >,1 = 1,17) 

WRITE (6, 226) CNIS 
WRITE(6,227) (CN6V ( I) , 1=1 ,20) 

WRITE(6,228) DEL I , DEL IN, RO, RON 
WRITE(6,229) CUT1,CUT2,Q( 1 ) ,04P 
WR I T E ( 6 ,2 30 ) 06P, 07P , E RET A, C 
WRITE (6,231 ) ERStDELLV, ( DELLGM ( I) , 1 = 1,7) 

WRITE(6,232) ETAD, 11,12 
I F( I FLAG ) 2 , 1 ,2 

C CONTROL LINE SPACING ON SHORT OUTPUT. 

43 K =K+ I 

IF ( K-6 ) 45 , 44, 45 

44 K = 1 

WR I TE ( 6,262 ) 

45 IT=IFIXIT) 

C WRITE SHORT OUTPUT. 

WRI f E (6 ,233 ) IT , CN4G , CN4S ,CN l ( I) , CNIS, CM (3) ,CN(5),CN(7),CN(6),CN(2) 
1 »CN(4> ,CN( 1 ) 

C TEST IF MEW DH, PI), I PR CARD TO H E READ. 

41 I F< I FLAG) 107, 1,107 
107 WR ITE (6 ,261 ) 

GO TO 2 

200 FORMA") ( E10.0 ,E1 5.5 , I 5 > 

201 FORMAT ( 1H0»3hDH , 1 P E 1 1 . 4 , 4X , 3HPU , E 1 1 . 4,4X , 4H I PK ,I?/1H1) 

202 FORMAT (F10.0,I5) 

204 ECRMAT(IH0,9NPRESSI)RF ,IPEU.4,17H N/M**2 TEMP. ,E11.4,14H K 

\ DENSITY ,E11.4,33H KG/M*«3 TOTAL NUMBER DENSITY ,E11.4,7H 1 /M* 

205 FORMAT (1H ,24HC0UL1'MB COMPRESSIBILITY ,FR.5,2RH RECIPROCAL DEBY 

IE LENGTH ,1PEU.4,43H 1/M EOUIV. CflNC. /CRITICAL FOUIV. CONC. , 
2E1I.4) * 

206 FORMAT ( 1H0»65HSPECI ES NUMBER DENSITY, PARTICLES/CUBIC MFTFR/1X,3H 
1H , 1 PE 1 1 .4,4X, 3HH+ , El 1 .4 ,4X ,2HE , E 11 .4 ,4X , 3HH2 , El 1 .4,4X , 3HH- , 
2E11.4,4X, 4HH2+ , E I 1 . 4 , 4 X , 4HH3+ ,E11.4) 

207 FORMAT U HQ , 53HDE TA I LFD NUMBER DENSITY FOR H2, MOLECUL ES/CUft I C METE 

l R / 1 X , 6 7H GROUND ELECTRONIC STATE WITH VIBRATIONAL QUANTUM .NUMBER O 
2E V = (J , 1PE1 1 .4,4X,5HV = 1 ,E1 1 .4 ,4X, 5HV=2 ,E11.4/1X,7H V = 3 , 

3E 1 1 . 4, 4X , 5HV=4 , E) 1 .4,4X , 5HV=5 , El 1 . 4, 4X, 5HV = 6 , E 1 1 .4 , 4X , 5HV= 7 

^ J 1 '2 *4X,5HV-8 ,EH.4/1X,7H V=9 , E 1 1 .4 ,4X , 5HV= 10 ,E11.4,4X,5H 

5V-11 ,E1 1 .4,4X, 5HV = 12 , E 1 1 . 4 , 4X , 5H V= 1 3 , E 1 1 . 4 , 4 X , 5HV= 1 4 ,E11.4) 

208 FORMAT ( 1H ,31H 1SNS SINGLET SIGMA STATES W I TH , 1 1 X , 5HN=1 ,IPE11.4, 

14X, 5HIM=2 ,E1 1 .4,4X,5HN = 3 , E I 1 . 4 , 4X, 5HN=4 ,E11.4) 

209 FURMATflH ,31H 1SNP SINGLET SIGMA STATES W I TH , 1 1 X , 5HN=2 ,IPE11.4, 

14X,5HN=3 , E 1 1 . 4, 4X, 5HN=4 ,E11.4) 

210 FORMAT! 1H ,28H 1SNP SINGLET PI STATES W I TH, 14X , 5HN=2 , 1PE11.4,4X, 

1 5HN= 3 , El 1 .4,4X,5HN=4 ,E11.4) 

211 FORMATIH ,3 1 H 1SN0 SINGLET SIGMA STATES W I TH , 1 1 X , 5HN = 3 .1PE11.4. 

1 4X , 5HN = 4 , E 1 1 .4 ) 

212 FORMAT ( 1 H ,28H 1SMD SINGLET PI STATES WITH, 14X ,5HN=3 ,1PE11.4,4X, 
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15HN=4 ,E11.4) 

213 FORMAT { 1 H ,31H 1$ND SINGLET DELTA STATES W I TH , 1 IX , 5HN=3 ,1PE11.4, 

1 4X , 5HN=4 , Ell. 41 

214 FORMAT I 1 H ,25H 1S4F SINGLET SIGMA STATE ,22X » 1 PEI 1 .4/ IX , 22H 1S4F SI 
1NGLET PI STATE, 2 5X , 1 PE 11 . 4/ 1 X , 25H 1S4F SINGLET DELTA ST AT E , 22 X , 1 P E 
2U.4/1X»23H 1 S4F SINGLET PHI ST ATE , 24X , E 1 1 . 4 ) 

215 FORMAT ( 1H ,31H ISNS TRIPLET SIGMA STATES N I TH, 1 IX, 5HN=2 ,1PE11.4, 

1 4 X , 5HN=3 ,E11.4,4X,5HN=4 ,611.4} 

216 FORMAT ( 1 H ,31H 1SNP TRIPLET SIGMA STATES WITH, 11X, 5HN= 3 ,1PEU.4, 

14X t 5HN=4 , E 1 1 .4 ) 

217 FORMAT ( 1 H ,28H iSNP TRIPLET PI STATES W I TH , 1 4X , 5HN=2 ,1PE11.4,4X, 

1 5HN= 3 ,E11.4,4X,5HN*4 ,611.4} 

218 FORMAT 1 1H ,31H 1SND TRIPLET SIGMA STATES W 1 TH , 1 1 X, 5HN= 3 ,1P£11,4, 

14X , 5HN=4 , E 1 1 .4 ) 

219 FURMATUH ,28H 1SMD TRIPLET PI STATES WITH, 14X,5HN = 3 , 1PE 1 1 .4,4X , 

1 5 HN =4 ,E11.4} 

220 FORMAT I 1 H ,31H 1SND TRIPLET DELTA STATES W I TH, 1 IX, 5HN=3 ,1PE11.4, 

1 4X , 5HN=4 , E 1 1 * 4 } 

221 FORMAT (1H ,25H 1S4F TRIPLET SIGMA ST AT E , 22 X , 1 P E 1 1 . 4/ 1 X , 2 2H 1S4F TR 
1IPLET PI STATE, 2 5X , 1 PE 1 1 . 4/ 1 X , 2 5H 1S4F TRIPLET DELTA S T A TE , 2 2X , 1 PE 
211.4/ IX, 23H 1S4F TRIPLET PHI STATE ,24X , E 1 1 .4 } 

222 FORMAT (1H ,?0H ALL ISM STATES WI TH ,22X , 5HN*5 , 1 PEI 1.4,4X»5HN=6 , 

1E1 1.4, 4X, 5HN=7 , E 1 1 .4 ,4X, 5HN=8 ,E11.4/1X,7H N=9 , E 1 1 . 4 , 4X , 5HN= 

210 ,EH.4,4X,5HN = 11 , El 1.4, 4X ,5HN*12 , E 1 1 .4 ,4X , 5HN= 1 3 ,611. 4, 4X, 
35HN= 1 4 , E 1 1 . 4 ) 

223 FORMAT ( 1 H ,7H N»15 , 1 P E 1 1 . 4 , 4X , 5HN= 1 6 , E 1 1 .4 ,4X , 5HN= 1 7 ,E11.4} 

224 FORMAT ( 1 H ,35H TOTAL IN EXCITED ELECTRONIC S T A T ES , 1 2 X , 1 PE 1 1 .4 ) 

225 FORMAT ( 1 HO ,4RH0E T A I L FC NUMBER DENSITY FOR H, ATOMS/CUBIC METER/ IX, 

1 47 H STATES WITH PRINCIPAL QUANTUM NUMBER OF N=1 , 1 PE 1 1 .4, 4X , 5HN* 

22 ,E11 .4,4X,5HN=3 , E 1 1 .4 , 4X , 5HN=4 ,E11.4/1X,7H N=5 ,E11.4,4X, 

35HN=6 ,E11.4,4X,5HN=7 , E 1 1 . 4 , 4X , 5HN = 8 , E 1 1 .4 , 4X , 5HN = 9 ,E11.4, 

44X , 5HN= 10 ,E11.4/1X,7H N=ll , E 1 1 .4 , 4X , 5HM= 12 , E 1 1 .4 , 4X , 5HN= 1 3 , 
5E11 .4,4X,5HN=14 , E 1 1 .4 , 4X , 5HN= 1 5 , E 1 1 . 4 , 4X , 5HN= 1 6 ,E11.4/1X,7H N* 
617 , F 1 1 .4 ) 

226 FORMAT ( 1 H ,35H TOTAL IN EXCITED ELECTRONIC ST A T E S , 1 2 X , 1 PE 1 1 . 4 > 

227 FORMAT ( 1 HO ,49HDE T A I L ED NUMBER DENSITY FOR H2+, IONS/CURIC METER/1X 
1 , 67H GROUND ELECTRONIC STATE WITH VIBRATIONAL QUANTUM NUMBER OF 

2 V=0 » 1 PE 1 1 .4,4X,5HV = 1 , E 1 1 . 4 , 4X , 5HV = 2 ,EU.4/1X,7H V = 3 ,E11.4 

3,4X, 5HV = 4 ,E11.4,4X,5HV = 5 » E 1 1 . 4 , 4 X , 5H V = 6 , F 1 1 . 4 , 4X , 5HV = 7 , 

4E11.4,4X,5HV=fi , F_ 1 1 . 4/ 1 X , 7H V = 9 , E 11 .4 , 4X , 5HV = 1 0 , E 11 .4 , 4X , 5H V = 
511 ,E11.4,4X,5HV=12 , E 11 . 4 , 4 X , 5HV = 1 3 , E U . 4 , 4X , 5HV = 14 ,E11.4/1X, 

67H V= 1 5 , R 1 1 ,4 » 4X , 5HV= 1 6 , £ 11 . 4 , 4X , 5H V= 1 7 , E 1 1 .4 , 4X , 5HV = 1 8 , 

7 E 1 1 * 4, 4X » 5HV = 19 ,611.4} 

228 FORMAT ( 1H0,24H MISCELLANEOUS QU ANT I T I E S/ 1 X , 10H DELI ,1PEU.4, 

1 8 X , 9HDEL I N ,EI1.4,8X,9HR0 , 6 1 1 .4 , 8 X , 9HR0N ,E11.4) 

229 FORMAT ( 1H ,10H CUT1 , 1 PE 1 1 .4 , 8 X , 9HCUT 2 , E 1 1 .4 , 8X , 9HQ t 1 } 

1 , Ell .4,8X, 9HQ4P ,E11.4) 

230 FORMAT 1 1 H ,10H Q6P , 1 PE 1 1 .4 , 8X , 9HQ7P , E 1 1 . 4 , 8 X ,9HERET A 

1 ,611 .4,8X,9HC ,011.4} 

231 FORMAT f 1 H ,10H ERS , 1 PE 1 1 . 4 , 8X , 9H0EL L V , D 11 . 4 , 8 X , 9HD6L LGM < 

1 1 } ,D1 1 .4,8X,9HDELLGM(2 ) , D 1 1 . 4/ 1 X , 1 OH DEL LGM I 3) , D1 1 . 4 , 8 X , 9HDEL LGM ( 4 

2 > ,D11.4,8X, 9HDEL LGM ( 5 ) ,D11.4,8X, 9HDELLGM (6) , Oil. 4/ IX, 1 OH DEL LGM (7} 
3,011.4} 

232 FORMAT UH+,29X,9HETAD , IP El 1 .4 ,fiX ,9HI 1 ,I3,16X,9HI2 

1 ,I3/1H1) 

233 FORMAT ( 1 H , I 5 , 2X , 1 P 1 1 E 1 1 .3 } 
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261 FORMAT (lHlf 

262 FORMAT <1H > 
END 



$ 1 BFTC LQ1PF 

C SUBPROGRAM FOR PARTITION FUNCTION OF H ATOMS 
C 

C THIS SUBPROGRAM FINDS THE CONTRIBUTION OF EACH EXCITED STATE TO THE 
C PARTITION FUNCTION USING 2 DIFFERENT CUTOFFS AND ADOPTS THE SMALLER 
C VALUE. THE CONTRIBUTIONS OF THE EXCITED STATES ARE THEN ADDED TD THE 
C GROUND STATE CONTRIBUTION TO GIVE THE PARTITION FUNCTION. 

FUNCTION QL PF ( RO»DEL I ) 

COMMON/BLOCK2/PHI , PHICtSUMl/BLOCKS/ARBl , ARB2 
DIMENSION PHI 1 17 J ,PH1C 1171 , STARH 18 ) , ROOT A I 18 ) 

DATA ST AR I / 10967 380. ,2741970., 1218650. ,685490. » 438720. , 304670. , 
1223840. , 171380. , 135410. , 109680. » 90650. , 76170. T 64900. t 55960 . , 
248750. » 42 8 50. .37960. , 33 360./ 

DATA R00TA/7.2744E-6, 1 .4549E-5 , 2 . 182 3E-5 , 2 .90986-5, 3 .6372E-5 , 
14.3646E-5,5.0921E~5, 5. 81956-5, 6. 54706-5,7.27446-5,8.00186-5, 
28.72936-5,9,45676-5,1 .0 1 84E-4 , 1 .091 2E-4 , 1 . 1 639E-4 , 1 .23666-4, 

31. 30946—4/ 

PHIC(l)=PHim 

SUM1=0. 

R OflT R= SORT ( RO^ ARB 1 ) 

DO 1 N = 2 ,17 

2 STAR1 = (STARKNUSTARI IN+1U/2. 

1 F( ST AR1 -DELI *ARB2 >5,4,4 

4 PH I C ( N ) = PH I { N I 
GO TO 9 

5 IF(STARIIN) —DEL I*ARB2 >7,6,6 

6 PHIC (N)=PHI (N)*( (STARI (N)-DELI : < ! ARR2)/ (STARI (N)-STAR I (M+l U + .5) 

GO TO 9 

7 STAR2=(STAP,I (N-l J+STARI ( N ) ) /2 . 

1F(DELI*ARB2-STAR2)8,3»3 

8 PHIC(N)=PHI (N>*U STARI { N 1 -DEL I* ARB2 )/ (STAR! I N-l ) -STAR I ( N ) ) + , 5 ) 

GO TO 9 

3 PH1C<N)=*0, 

9 ROOT 1 = {-R0GTA(MI-RQCTA{N+1 ) )/2. 

I F ( ROOT 1+R DOT R ) 11,10,10 

10 PH 1CN= PH I { N ) 

GO TO 16 

11 I F (—ROOT A ( N > + ROOT R ) 13,12,12 

12 PHICN= PHI (N)=M 1 — ROOT A ( N ) +ROOTR ) / ( — ROOT A ( N ) +ROOT A ( N+ 1 ) ) + .5) 

GO TO 16 

13 RQQT2= I —ROOT A ( N— 1 ) — ROD TA ( N ) > / 2 , 

I F { -ROOTR— R00T2 114, 15,15 

14 PHI CN= PHI (N)M ( -ROOT A ( M ) + ROOTR ) / t -ROOT A(N— l)+ROOTA(MJ ) + .5) 

GO TO 16 

15 PH1CN- 0. 

16 PHIC<N)= AMIN1(PHIC(N>,PHICN} 

SUM1= SUM 1+PH I C ( N ) 

1 CONTINUE 

OlPFsPH ECU l+SUMl 

RETURN 

END 
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SIBFTC L04PF 

C SUBPROGRAM FOR PARTITION FUNCTION OF H2 MOLECULES 
C 

c this subprogram finds the contribution of each excited state to the 
C PARTITION function using 2 different cutoffs and adopts the smaller 
C VALUE. THE CONTRIBUTIONS OF THE EXCITED STATES ARE THEN ADDED TO THE 
C GROUND STATE CONTRIBUTION TO GIVE THE PARTITION FUNCTION. 

FUNCTION Q4PF ( RO » DE L I I 

CUMMON/BLOCK3/SUM,GAM t GAMC, SUM4/BLOCK4/ I MAX/BLQCK5/ARB 1 , ARB2 
DIMENSION GAM (6, 23) ,GAMC ( 6,23) , 1MAX (23 > , STAR I I 7,23) , RUQTA I 7 , 23 I 
DATA! (STARK I , J ) , I * 1 , 5 ) , J* 1 , 3 ) / 1244 1 300 . , 252492 0. , 11 45640 . , 593 1 80 . 

1.458380.. 12441300. .3420970.. 1393480. .752770. .458380.. 12441300. , 

22531970. . 1154100. . 65 7930 .. 45 8380. / 

DATA! ( STAR I (I , J ) , 1 = 1 , 4 ) , J =4 ,6 ) /2 892650 . , 1 260140 . , 7 1 1 820 . ,4 583 80 . , 

12892650. . 1240260. .685140. .458380. .2892650. . 1 204780 677040 . , 
2458380. / 

DAT A ( (STARK I , J > , I = 1 , 3 ) , J = 7 , 1 0 ) / 1 288480 . , 7 33350 . ,458380. , 12 88480., 

1652200. . 458380. . 12 88480 724700 . . 458380 . . 12 88480 . . 70 1 600 . . 458 380 . / 
OAT A (STAR I (1,11), I = L , 5 )/ 1 244 1 300 ., 2942520 1 262970. , 706510 . , 

1458380./ 

DATA ( STAR K I , 12 ) , 1 = 1,4)/ 2892650., 1781960. , 889930. , 458380. / 

DATA ( STAR Id, 13), 1 = 1, 5) / 1 244 1 300 . , 296 1 190 . , 1 2 80330 . , 7305 70. , 
1458380./ 

DAT At ( STAR! (I , J) ,1 = 1 ,4) , J = 1 4 , 1 6 ) / 2892 650 . , 1 2 70300 ., 703460 . ,4583 80 . 

1.2892650. . 1248340. .693310. .458380. .2892650. . 1 2 1 1 650 . . 692 720 . , 
2458380./ 

DA T A ( (STARK I , J ) , 1=1 ,3) , J=1 7 ,20) / 1288480. ,696660. , 458380. , 12 88480. 

1 .682150. . 458380. . 1288480. .724760. .458380. . 1 2 88480 . . 70 1 660 . , 
2458380./ 

DA T A < ( STAR I ( I , J ) , 1 = 1 ,7) ,J=21,22 ) / 708 7 30 . , 45 8 380 . , 3 1 7 1 30 . , 2 32 380 . , 

1177660.. 140210. . 1 13520. . 140210. . 1 L 3520. ,93760. , 78800 ., 67 200 . , 

257980. . 50540./ 

DATA ( STAR 1(1,23), 1=1 , 5 > / 5 7980. , 50540 . , 44450 394 10 . , 35190./ 

DATA! (ROOT At I ,J) ,1 = 1,5) , J =1 , 3 > / . 6803 E-5 , 1 , 454E-5 , 2 , 1 93 E-5 , 

1 3.020E-5, 3.575E-5, . 68 03E -5 , 1 . 3 1 2 F-5 , 1 . 91 9 E-5 , 2 . 656 E-5 , 3 . 575 E-5 , 
?.6803E-6,1.452E— 5,2.1 82E-6 , 2 .88? F-5 , 3. 575E-5/ 

DAT A ( { ROOT A ( I , J ) , I = 1 , 4 ) , J =4 , 6 ) / 1 . 4 1 5 F-5 , 2 . 1 35 E-5 , 2 . 86 1 E-5 , 

13.575E-5, 1 .415E-5,2. 1 50E-6 , 2 .899F- 5 , 3 . 575 E-5 , l .4 1 5 E -5 , 2 . 1 94E-5 , 

? 2.916 F— 5,3.575 E-5/ 

DATA t (ROUTA (I, J), 1=1, 3), J =7, 10)/?. 1 20 E-5 , 2 . 862 E-5 , 3 . 575 F-5 , 

12 . 1 20 E-5,2.934E -5,3.57 5E-5,? . 120E-5,2. 842 E-5, 3 , 5 78F-5 , 2 . 1 20E-5 , 
2?.«63E-6,3.575E-5/ 

DA T A ( ROOTA ( 1,11), I = 1 , 5 ) / . 6803E-5 , 1 . 4 54E- 5 , ? . 1 93E-5 , 3 . 020E- 5 , 
13.575E— 5/ 

DA TAt ROOTA ( 1,12), I * 1 , 4 ) / 1 . 41 5E-5 , 1 . 9 1 9F-5 , 2 . 65 6E-5 , 3 . 576E-5 / 

DAT A { ROOTA { I , 13) , I = 1 , 5 ) / . 6803 F-5 , 1 . 452E-5 , 2 . 1 82. E-5 , 2 . 882 E-5 , 
13.575E-5/ 

DAT A ( { ROOTA ( 1,0) ,1 = 1,4) , J = 1 4 , 1 6 ) / 1 .41 5E-5 , 2 . 1 3 5E-5 , 2 . 86 1 E-5 , 
13.576E-5, 1.415E-5,2. 1 50E-5 , 2 , 899 E-5 , 3 . 57 5 E-5 , 1 .41 5 E-5 , 2. L 94E-5, 

2 2.9L5E-5,3.575E-5/ 

DA TAt { ROOTA ( I , J > , I =1 , 3 ) , J = 1 7 , 20 > /2 . 1 20E-5 , 2 . 862 E-5 , 3 . 57 5E-5 , 
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12.120E-5,2*934E-5,3.575E-5,2.120E-5,2.842E-5,3.575E-5,2..120E-5, 

22* 863E-5,3,575E— 5/ 

DATA l (ROOTAt I, J), 1 = 1, 7), J =2 1,22 ) /2 . 8 6 1 E-5 , 3 . 5 7 5 fi~5 , 4 . 302 E-5 , 
15.030E-5,5.758E~5,6.48SE-5,7.212E-5,6.485E-5, 7 .2 1 2 E-5 , 7, 940E-5 , 

28 *668E-5,9,395E-5,10. 12E-5, 10.8 5E-5 7 
OA T A ( ROOT A { I , 2 3 ) , I = 3 . , 5 ) / 1 0 . 1 2E-5 * 1 0 . 85E-5 , 1 1 . 5 8E-5 , 12 . 30E-5 , 
113.03E-5/ 

ROOT R= SORT ( RG*ARB 1 > 

SUM4=0. 

DO 1 M=l,23 
I MA= I MAX ( M ) 

DO 1 1=2, IMA 

2 STAR1=( STARK 1 ,M)+STAR I( I + 1,M) ) / 2. 

IF(STAR1-0ELKARB2)5,4,4 

4 GAMC ( I , M ) =G AM ( I ,M > 

GO TO 9 

5 I P( START ( I ,M 1-DELI *ARB2 ) 7,6,6 

6 GAMCt I ,M)=GAM( I ,M)*< ( STAR I( I ,M l-DEL I$ARB2 ) / ( STAR I < I,M)-STARKK1, 
1M 1 ) + * 5 1 

GO TO 9 

7 ST AR2=( STAR I < 1-1 ,MH-$TARK I ,M> )/2. 

I F ( DEL I*ARB2-STAR2 >8,3,3 

8 G AMC ( I , M ) =GAM ( I , M ) * ( ( ST AR I ( I , M ) -DEL I * AR82 > / (STARHI-1,M>-STARI(I,M 
1 > > + * 5 > 

GO TG 9 

3 GAMC ( 1 , M > =0 * 

9 ROOT 1 = 1 —ROOT A ( I , M 1 -ROOT A ( I + 1,M) > /2. 

IE ( R0DT1+R00TR ) 11,10,10 

10 GAMC I M= G AM ( I ,M> 

GO TO 16 

11 I F ( — ROO TA ( I ,M ) -t-ROOTR > 13,12,12 

12 GAMC IM= GAM ( I , M > - ( (-ROOTA< I , Ml+ROOTR > / ( -ROOTAt I ,M> +ROOTA ( I+1,M>| 

1 + . 5 > 

GO TO 16 

13 R00T2= ( —ROOT A ( I — 1 , M ) -ROO TA( 1 ,M> )/2. 

IF (-R00TR-R00T2 > 14,15,15 

14 GAMC 1 M= G AM ( I ,M)*( (-R00TA ( I ,M ) +R00TR1 / (-ROOTA ( 1-1 ,M ) +ROOTA ( I ,M) > 

1 + .5) 

GO TO 16 
L 5 GAMC I M= 0. 

16 GAMC ( I , M ) = AM INK GAMC (I , M ) , GAMC I M > 

SUM4= SUM4+GAMC ( I ,H) 

1 CONTINUE 

04PF=SUM+ SUM4 

RETURN 

END 
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£ ! FSFTC LRLO 

SUBPROGRAM GIVING CONSTANTS FOR H, H+, E~, H2» H- t H2+, AND H3 + 

BLOCK IT A T A 

COMMON/ BLOCKA7 I MAX/BL0CK6/Z 

1/BL0CK7/ALPHA,EH,EF,ED,EH2,'JEIGHT,RE,BETA, JP 1M,KP1 M, GH2 ,GH2P, FH2 
DIMENSION I MAX (23 I ,7. (7) .ALPHA 17) , EH < 17>,EF(6,23),ED(6,23),EH2(6, 
123 > .WEIGHT (6,23) ,RE (6,231 T BETA! A, 23 1, JPlfM 15 1 ,I<P1M( 20) ,GH2( 15} , 
2GH2P(20I ,FH2( 15,39) 

DATA I MAX/4, A, A, 3, 3, 3, 2, 2, 2, 2, 4, 3, A, 3, 3, 3, 2 ,2, 2, 2 ,6, 6, 4/ 

DATA Z/Q„, l.,l.,0. ,1. ,1 

DATA ALPHA/3. 1673E-3,3*1 700E— 3 » 2.49*-0E+2 , 1. 1 1 986-3 , 3 . 1 649E- 3 , 

11 .1203E-3.6.0974E-4/ 

DATA EH/ 0., 118354. , 140271. , 147942 i 51493. ,153421., 154534, ,155339. 

1 .155856.. 156227. .156500. . 1 56709 .. 1 5 6571 . , 157000 157103., 1571^8. , 
2157259./ 

DATA EE/O. ,3249. , 32 34 . , 2 7 38 , , 3*G . , 19 18. ,2332., 2935. ,3*0. ,3475. , 

13236. . 3164. .3*0. ,3342 . ,3269. ,4*0. , 32 B6 . , 3064. , 4*0. , 3 1 2 0. , 3 1 40 . , 4*0 

2. . 2 987. .5*0. ,3359. ,5*0., 29 79. ,5*0., 3340. , 5*0 ., 3672 ., 3264 ., 2823 . , 
33*0. ,3144. ,3085. ,4*0. , 3352., 3224., 3394. , 3*0. , 30 86. , 30 87 . ,4*0. , 

43271 . . 31*8. ,4*0. ,3100. ,3 109., 4*0. ,2953. , 5*0 . , 33 2 5 . , 5*0 . , 2979 . , 5*0 . 

5.3340. . 5*0., 304 5. ,3067. ,3085. ,3097. ,3108. ,0., 311 6. ,312 3. ,3127., 

631 30. . 31 33. .0. .31 36. .31 38.. 3140. . 2 *0 . / 

DATA ED/O., 29 339., 3 1396., 24349., 3*0. ,41506. ,12655. ,26748. ,3*0., 

1297 15. . 31 520. .32012. .3*0., 14779. ,26344. ,4*0. , 14*64 . ,25847 . ,4*0 . , 

232183. . 32 2 73.. 4*0. ,2649 7., 5*0., 25*46. , 5*0 . , 26 369 . , 5*0. , 32734. , 5*0. 

3. 35561.. 3309 7.. 2 602 2. .3*0., 18642., 10362. ,4*0. ,35661. ,33324., 

426683.. 3*0., 16026., 26123. ,4*0., 15826, , 26083. , 4*0. , 32272. , 32*80. , 
54*0. ,2 5951 . ,5*0. ,32**6. ,5*0., 2 6 369. ,5*0. ,32 735. ,5*0. ,29363., 

630334.. 30908.. 31236.. 3 1*78. ,0. ,31661 . , 3178*. , 31882. ,31961. ,32025. 

7. . 0.32072.. 32 112. .32 147.. 2*0./ 

DATA EH2/0. , 14097 5. , 1608 35. , 1690 35.1 3*0. ,1288 07,,). 57 658. ,1666 36., 
13*0., 140751., 16071 1 . , 167890., 3*0., 159016. , 16704 1 . ,4*0. , 15933 1 . , 

2 167537. . 4*0. , 160048 1 6 1629 . ,4*0 . , 166888. , 5*0 . , 167843 . , 5*0. , 

31 670 16.. 5*0., 167 168., 5*0., 134752.t 159133., 167363, ,3*0., 151672., 

4164401. . 4*0., 134652. ,158907., 166701 . ,3*0. , 159045. , 1672 61 . , 4*0. , 

5159245. . 167301 ..4*0., 159959., 1674? 2, ,4*0. ,167433. ,5*0. ,167466., 
65*0. , 1670 L 5. , 5*0. , 1 67 1 67. , 5*0. , 170820. , 172843. , 1 7*053 . , 1 74835. , 

7175369. . 0 . . 175749. . 176030. . 176243. .176408.. 176540. .0. .176645., 

8 176732. . 1 76803. .2*0./ 

DATA WE IGHT/O. , 3*1 . , 3*0. ,3*1 . ,3*0. , 3*2 . , 3*0 . , 2* 1 . ,4*0 . , 2*2 . , 4*0 . , 
12*2., 4*0., 1., 5*0., 2., 5*0 .,2., 5*0., 2., 6*0. ,3*3. ,3*0. ,2*3., 4*0., 

23*6 . ,3*0. ,2*3. ,4*0 . ,2*6 . ,4*0 .,2*6. ,**0. ,3. ,5*0. ,6. ,5*0. ,6. ,5*0., 

36. . 5*0. , 100. ,144. , 196. ,256. ,324. ,0. ,400. ,484. , 576. , 676. , 784 . ,0. , 

4900. . 1024. . 1 156. .2*0./ 

DATA RE/0.,1.0334E-10, 1. 1181E-10, 1.129<lE-l(), 3*0 . , 1 . 2 9076- 1 0 , 
11.1270E — 10,1, 111E-10, 3* D..1.0316E — 10,1. 0440E- 10,1. 0555F-10, 3*0., 
21.0740E— 10, 1.42 05 £—10,4*0. ,1.1639E-1 0, 1.3972 8-1 0 , **0 . , , 99 1 6E- 1 0 , 
3.89598-10,4*0. , 1 . 0809 E- 10 , 5*0 . , 1 .058 3 E- 1 0 , 5*0 . , I . 08? IE- 1 0 , 5*0 . , 
41.02538-10,5*0.,. 99556-10, 1 .0613E-10, 1 . 1 132E-10, 3*0. , 1 .0869E-10, 

51 .06 54E — 10,4*0. , 1.06996-10,1 .06 67F-10, ] .0603E-1 0 , 3*0 . , 1 .4461 E-10 , 
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6.92986-10,4*0. ,1 . 1466E-10, 1. 1675 E-l 0,4*0. ,.9937E-10 t . 0348 E- 10, 

74*0. ,1. 08668-10, 5*0., 1.027A8-10, 5*0., 1.0R21 8-10,5*0., 1.D253C-10, 
85*0. » 1.0732E-10,! .07038-10, 1.0676E-10, 1.06578-10,1.06418-10,0. , 

91 .Q628E-10,1.0617E-in, 1 .06108-1 0 , 1 . 0605E-10, 1.0600E-10 , 0 . , 

A 1 .05958—1 0,1 .05928-10, 1.05898-10,2*0./ 

DATA BE7A/0. ,1.4529 8+10,1.39078+10,1 .34488+10 , ?*0. , .694 8 E + 10 , 
l 1. 70708+10, 1.375OE+10, 3*0., 1.54938+10, 1.33956+10,1.34518+10,3*9. , 
22. 2 7 70E+ 10, 1. 5551 E+l 0,4*0., 2.26588+10,1 .46628+10,4*0. ,1 .32048+10 , 
3 1.328 18+10, 4*0., 1.40688+ 10, 5*0., L . 6348 6+ Hi , 5*0 . ,1.40688+10,5*0., 
41. 40688+10, 5*0., 1.48448+10, 1.364 18+10, 1.33858+10,3*0. ,1 .83098+10, 
52. 66948+10, 4*0., 1.34978+10,1. 34076+1 0, 1 .6084E+10, 3*0., 1.9680E+10, 
61 .46938+10,4*0. ,2. 1220E+10, 1 . 50228+ 10 ,4*0 1 .3 1 03 E+ 1 0 , 1 . 30868+10 , 
74*0. ,1. 4068E+1 0,5*0., 1 .4068E+ 10 , 5*0. , 1 . 40686+1 0 , 5*0 ., 1 . 4 068 E+ 1 0, 
85*0., 1.35568+ 10, 1,34068+ 10, 1.33 52 8 +10,1 .33308+10,1.33198+10,0., 
91.33146+10,1.33118+10, 1.3309E+10, 1.33088+10, 1.33078+10,0. , 

A3*l .33078+10,2*0./ 

DATA JPlM/39,36,34,32 ,30,28,26, 24 ,22,19,17, 15, 12,9,5/ , K81 M/42 , 8 0 » 
138,37,35,33,31 ,30,28,26,24,22,20,18,16, 14,11 ,8,5,2/ 

DATA GH2/0. ,41611 4. ,808574., L 178092. , 1 524910 . , 1 849093 . , 2 150499. , 

22 42 8 77 7. . 26831 39. .2912441. .3115129. .32 88751. .3430366. .3535284., 
23597480. /,GH2P/0.,2 191 34., 425545. ,619656. ,801821. ,972336, , 

31 131427.. 1279260.. 14 15937. .1541490. .165 6ft 79., 175R991 .,1850623., 

41930484. . 1998171. .2053168. .2094829. .2122424.. 2 135548.. 2137849./ 



DATA ( FH2 

( I, 1),I= 

It 15)/ 





X 

0., 

0., 

-0., 

— 0 * , 

-0., 

0. 


X 

0 . , 

“0., 

o. t 

“0 . ? 

"0 . , 

0 . 


X 

-0., 

-0., 

0./ 






DATA ( FH2 

( r , 2 ) , i = 

1,15)/ 





X 

11852. , 

11257., 

10679* , 

10111., 

9565. , 

8988 . 


X 

8413., 

7821., 

7211*, 

6 6 67 . , 

5849. , 

5068. 


X 

4152., 

3079., 

1727./ 






DATA ( FH2 

(It 3) , 1= 

It 15)/ 





X 

35453. , 

33672., 

31932., 

30233., 

28582., 

26865. 


X 

25136., 

23363., 

21527., 

19679 . , 

17463., 

15123. 


X 

12365. , 

9114., 

4994./ 






DATA t FH2 

f I , 4) , 1 = 

1,15)/ 





X 

70587., 

67029., 

63564. , 

60171., 

56881. , 

53456. 


X 

50009., 

46467. , 

42801*, 

38910. , 

34663. , 

29976. 


X 

24489. , 

17879., 

9439./ 






DATA t FH2 

< I , 5 > » I = 

1,15)/ 





X 

116932., 

111037., 

105286., 

99643. , 

94 1 62 . , 

88500. 


X 

82760., 

76875. , 

70737., 

642 79. , 

57196., 

49365. 


X 

40218. , 

29090, , 

14500./ 






DATA t FH2 

( I , 6 > , I = 

1,14)/ 





X 

174094., 

165299., 

156715., 

148294. , 

140077* , 

131651, 


X 

123062., 

114259., 

105089. , 

95367. , 

86730., 

72917. 


X 

59127. , 

42240./ 







DATA < FH2 

(1, 7 > * I — 

1,141/ 





X 

241600. , 

229358., 

217404. , 

205686., 

194196. , 

182501* 

, 

X 

170530., 

158230. , 

145398 . , 

131803. , 

116864., 

100218, 

f 

X 

80722. , 

56691 ./ 







DATA { FH2 

( I , 8) , I = 

1,14)/ 





X 

318852., 

302680., 

286854. , 

271338. , 

256050. , 

240547. 

, 

X 

224697., 

208344. , 

191254. , 

173066. , 

153109. , 

130746. 


X 

104418, , 

71609,/ 







DATA ( FH2 

(1, 9) ,1 = 

1,14)/ 
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X 

405346* * 

384685. t 

364505., 

344702* , 

325119., 

305303. , 

X 

2 8 3058 * , 

264112* f 

242156., 

218740*, 

192981*, 

16394-7 . T 

X 

129518*, 

85558./ 





DATA ( FH2 

( 1,10) » I 

= 1,13)/ 




X 

5QG392 * , 

474767 * , 

449781., 

425219., 

400851*, 

376206. , 

X 

351079. , 

325015* f 

297609., 

268282., 

235916. , 

199246. , 

X 

155153./ 






DATA ( F H2 

U , 1 1) f I 

= 1*13)/ 




X 

603335. , 

572282. t 

542057*, 

512290*, 

482679. , 

452715. , 

X 

4-22177., 

390486 * , 

357056., 

321116* , 

281300., 

235841 . , 

X 

180225./ 
DATA ( FH2 

f I, 12) , 1 

= 1,13)/ 




X 

713500. , 

676570. , 

640698. , 

605301., 

569997*, 

534273. , 

X 

497810, , 

459917* , 

419966., 

376580* , 

328393. , 

272721 . , 

X 

20304-2./ 
DATA { FH2 

U f 1 3 ) , I 

= 1,12)/ 




X 

830194. , 

786987., 

745062. , 

703666* , 

662221. , 

620288. t 

X 

577390., 

532733. , 

485385., 

433965., 

376318., 

308806./ 


DATA ! FH2 

u , 1 4 ) , I 

= 1,12)/ 




X 

952734., 

902894 * , 

854558*, 

806738*, 

758744. , 

710158. , 

X 

660331 . , 

608326 • , 

553002. , 

492575., 

424316. , 

342614./ 


DATA ( FH2 

t F , 15) f I 

= 1, 12)/ 




X 

1 08044 L , , 

1023684*, 

968516. , 

913920., 

853987. , 

803306. T 

X 

746059 . , 

686133., 

622111 

551738*, 

471342., 

371564./ 


DATA < FH2 

U , 1 6 ) , T 

= 1*11)/ 




X 

1212683. , 

1 148721., 

1086399. , 

1024604- , 

962369. , 

899130. , 

X 

833977. , 

765632*, 

691994* , 

610507. , 

515770./ 



DATA ( FH2 

U ,17 > , I 

= 1,11)/ 




X 

1348894. , 

1277402*, 

1207582. , 

1138229., 

1068309*, 

997097. , 

X 

923479. , 

845805. , 

761758. , 

667760* , 

555461 ./ 



DATA ( FH2 

(I , 18), 1 

= 1 , 10)/ 




X 

1488465. , 

1 4091 42. , 

1331524*, 

1284241., 

1176241., 

1096586. * 

X 

1013929. , 

926305., 

330773., 

722407./ 




DATA { FH2 

H , 1 9 } , I 

=1,10)/ 




X 

1630829. , 

1543379., 

1457681., 

1372094*, 

1285621. , 

1197038., 

X 

1104753,, 

100632*.* 

897973., 

772062./ 




DATA ( FH2 

< 1 ,20) , I 

= 1, 9)/ 




X 

1775451 ., 

1679592. , 

1585510., 

1491280. , 

1395881., 

1297904. , 

X 

1195296. , 

1084950* , 

96192 1 ,/ 





DATA { FH2 

(I , 2 1 > , I 

=1, 9)/ 




X 

1921791 . , 

1817272., 

1714504*, 

1611313*, 

1506573*, 

1398545. , 

X 

1284780. , 

1161330. , 

1020958./ 





DATA { FH2 

( f ,22 ) , I 

= 1* 9)/ 




X 

2069372 . , 

1955961 . , 

1H44169., 

1731686* , 

1617105., 

1498380. , 

X 

1372526. , 

1234226*, 

1070223*/ 





DATA ( FK2 

II ,23) , I 

= 1* 8)/ 




X 

2217699. , 

2095209*, 

1974039. , 

1851871., 

1726957., 

1596827. , 

X 

1457534. , 

1301644./ 






DATA ( FH2 

( I ,24) , I 

= 1, 8)/ 




X 

2366299. , 

2234566. , 

2103658., 

1971389. , 

1835555., 

1693041 . , 

X 

1538665. , 

1359827./ 






DATA ( FH2 

t I *2 5) * I 

=1* 7)/ 




X 

2514773., 

2373579., 

2232613., 

2089735, , 

1942227*, 

1786214. , 

X 

1614187./ 
DATA < FH2 

( I ,26), 1 

= 1, 7)/ 




X 

2662745. , 

251 1326., 

2360462., 

2206391 *, 

2046374., 

1875100. , 


X 1680320./ 

DATA ( FH2 <1,271,1=1, 6)/ 

X 2809896., 264-8894., 2486726., 2320839., 2147107., 1958130./ 
DATA < FH2 (I, 28), 1=1, 6)/ 

X 2955879., 2784347., 2610949., 2432370., 2243405., 2031554./ 
DATA ( FH2 (I, 29), 1=1, 5)/ 

X 3100329., 2917813., 2732578., 2540291., 2333464./ 

DATA ( FH2 (1,301,1=1, 5)/ 

X 3242913., 3048916., 2851009., 2643523., 2414117,/ 

DATA ( FH2 (1,311,1=1, 41/ 

X 3383239., 3177186., 2965595., 2740537./ 

DATA ( FH2 (1,321,1=1, 4)/ 

X 3520920., 3302143., 3075429., 2828060./ 


DATA i FH2 

( I ,33) *1 = 1, 

3)/ 

X 3655593* * 

3423200., 3179162./ 

OAT A ( FH2 

(1,345,1=1, 

3)/ 

X 3786839. t 

3539599., 3273966./ 

DATA ( FH2 

( I ,35) ,1-1, 

2 ) / 

X 3914204., 

36500S 1 * / 


DATA ( FH2 

( 1,36) f 1 = 1, 

2)/ 

X 4037023., 

3752564./ 


DATA { FH2 

( I ,37) ,1 = 1, 

1 )/ 

X 4154504./ 

DATA { FH2 

{ f ,38) ,1=1, 

1)/ 

X 4265188./ 

DATA ( FH2 

( 1,39) ,1=1, 

1 )/ 

X 4364780,/ 


END 
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SI8FTC LBIP 

C SUBPROGRAM GIVING CONSTANTS FOR H2+, H3+, AND DIATOMIC PARTITION 
C FUNCTION APPROXIMATION USED FOR H2* 

C 

BLOCK DATA 

COMMON/ 8LOCK8/FH2P ,07 PT ,TT,THETAT,W0T,W1T,W2T 

DIMENSION FH2P ( 20,42 ) ,Q7PT ( 67 ) , TT ( 67 ) , THETAT ( 29 ) , WOT ( 29 > , WIT ( 29 ) , 
IW2TI29) 


DATA { FH2 H 

U, 11*1=1*201/ 




X 

-0., 

— 0 . * ■" 0 » * 

0* , 

0., 

-0*, 

X 

0.* 

— 0 . * —0 * » 

0., 

0., 

-o.. 

X 

0 » T 

—0 * * 0 . * 

o.. 

0., 

O.t 

X 

-o.. 

0./ 




DATA I FH2P 

U* 21*1=1,20)/ 




X 

5826. t 

5519., 5221., 

4936 # t 

4656., 

4384., 

X 

4110., 

3845 * , 3579., 

3313* , 

3044. , 

2770* , 

X 

2484*, 

2193., 1884., 

15154. , 

1192., 

790. , 

X 

342 

63./ 




DATA < FH2P 

II , 3), 1 = 1 ,191/ 




X 

17418, , 

16511., 15626., 

14773* , 

13951 . , 

13101*, 

X 

12309, , 

11501., 10702., 

9896., 

9106., 

8268. , 

X 

7432.* 

6548., 5624, , 

4627. t 

3542., 

2328., 

X 

950./ 





DATA { FH2P 

U , 41,1=1,19)/ 




X 

34714., 

32891., 31125., 

29420. , 

27761., 

26098., 

X 

24499., 

22896., 21305., 

19702., 

1 H 1 10 . , 

1 6462 . , 

X 

14769. , 

13006., 11153., 

9159*, 

6981. , 

4537* , 

X 

1771./ 






DATA ( FH2P 

(I , 5) ,1=1,19) / 




X 

57562., 

54529., 51897., 

48763. , 

45996. , 

43255. , 

X 

40588., 

37931., 35288., 

32627. , 

29971., 

27233., 

X 

24414., 

21477., 18367., 

15093* , 

11405. , 

7290*, 

X 

2567./ 






DATA t PH2P 

(I , 61,1 = 1,18 1/ 




X 

85804. , 

81259., 76880., 

72653*, 

68504., 

64428 . , 

X 

60438., 

56481., 52536., 

*8568. , 

44 569. , 

40470. , 

X 

36248. , 

31846., 27178., 

22158., 

16674. , 

10415./ 


DATA { FH2P 

U » 7) , 1=1 , 18)/ 




X 

119199. , 

112880., 106783., 

100893* , 

95115., 

89449. , 

X 

83891 

78374., 72874., 

67340., 

61754., 

56032., 

X 

50126. , 

43966., 37415., 

30368* , 

22612. , 

1366 7,/ 


DATA ( FH2P 

( I , 8) , 1=1, 13)/ 




X 

157500., 

1^9161., 141087., 

133252., 

125623. , 

118125., 

X 

110752., 

103434., 96132., 

88791 * , 

81362., 

73750* , 

X 

6 5887 * , 

57670., 48936., 

39471., 

28995. , 

16648./ 


DATA { FH2P 

t I , 91,1=1,171/ 




X 

200508., 

189847., 179545., 

169572*, 

159802., 

150242* , 

X 

140817., 

131465., 122124., 

11271b*, 

103195. , 

93443. , 

X 

83341. t 

72772., 61500., 

49244. , 

35545./ 



DATA ( FH2 P 

( I , 10) ,1=1*17)/ 




X 

247871., 

234658., 221891., 

209501., 

197419., 

185550* , 
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X 173848., 

162231. , 

150619., 

138934., 

127075., 

114906., 

X 102289.* 

89046., 

74899., 

59454., 

41872./ 


DATA ( FH2P 

U, 111,1 = 

1,17)/ 




X 299287.* 

283321., 

267857., 

252860., 

238213., 

223821. , 

X 209629.* 

195526., 

181418., 

167178., 

152762., 

137909., 

X 122509., 

106286., 

88881., 

69717., 

47436./ 


DATA ( FH2P 

(1,121,1= 

1,16)/ 




X 354507., 

335523., 

317158., 

299338., 

281915., 

264794* , 

X 247890., 

231081., 

214267., 

197292., 

180031., 

162258. * 

X 143766.* 

124214., 

103156., 

79723./ 



DATA { FH2P 

(1*13) ,I*» 

1,16)/ 




X 413131., 

390969., 

369504., 

348652., 

328260. , 

308204*, 

X 288397., 

268687. , 

248932., 

228981 . , 

208650. , 

187679 * * 

X 165801.* 

142591., 

117391., 

88908./ 



DATA ( PH2P 

(1,141,1= 

1,16)/ 




X 474906.* 

449325., 

424587. , 

400526., 

376973. , 

353795* * 

X 330899., 

308056., 

285182., 

261999. , 

238367. , 

213926* * 

X 188320.* 

161065., 

131180., 

96138./ 



DATA ( FH2P 

(I ,15) , 1 = 

1,15)/ 




X 539480.* 

510340. , 

482122., 

454676. , 

427770. , 

40 1 2 95 * f 

X 375092., 

348967., 

322728., 

296142. , 

268931 . , 

24073b * T 

X 211086.* 

179298., 

143982./ 




DATA ( FH2P 

( I ,16) , 1 = 

1,15)/ 




X 606555., 

573682., 

541824., 

510823., 

480405., 

450453* * 

X 420801.* 

391179. , 

361384. , 

331113. , 

300078. , 

267822., 

X 233725., 

196869., 

154868./ 




DATA 1 FH2P 

( 1 , 1 7 ) , I = 

1,14)/ 




X 675813.* 

639063., 

603420. , 

568673., 

534604. , 

501009* , 

X 467 684. * 

434374., 

400829., 

366666. , 

331539., 

294877*, 

X 255913., 

213220./ 





DATA ( FH2P 

(I ,18) *1 = 

1,14)/ 




X 746961.* 

706187., 

666609. , 

628013., 

590104. , 

552705*, 

X 515559., 

478369. , 

440836. , 

402539. , 

363021 

321573. , 

X 277166., 

227433./ 





DATA ( FH2P 

( I , 19 ) , I = 

1,13)/ 




X 819714., 

774786., 

731121. , 

688521 . , 

646665. , 

605298. , 

X 564153.* 

522908., 

481149., 

438458., 

394216. , 

347533*, 

X 296924./ 






DATA ( FH2P 

(1,201,1= 

1,13)/ 




X 893796., 

844572. , 

796731., 

750004, , 

704026. , 

658540, , 

X 613226., 

567696. , 

521531., 

*74 135., 

426 794. , 

372301* , 

X 314226./ 






DATA ( FH2P 

( 1 , 2 1 ) , 1 = 

1,12)/ 




X 968947., 

915328., 

863177., 

812177., 

761957., 

712195* , 

X 662532., 

612526. , 

561662. , 

509262, , 

454356. , 

395284./ 

DATA ( FH2P 

(I ,22) , 1 = 

1,12)/ 




X 1044880., 

986793. , 

930222., 

874836. , 

820241 ., 

766038. , 

X 711857., 

657127., 

601295 * f 

543481., 

482431 . , 

41 5404*/ 

DATA ( PH2P 

(1,231,1= 

It 11 >/ 




X 1121411., 

1058739. , 

997639* , 

937744., 

878605. , 

8198 l 3* , 

X 760888., 

701235., 

640093 • , 

576436. , 

508364./ 


DATA < FH2P 

(1,24),!= 

1,11)/ 




X 1198293., 

1130945. , 

1065209*, 

1000683. , 

936873. , 

87 3 310. , 

X 809445., 

744572., 

677761., 

607587., 

531009./ 


DATA l FH2P 

(1,251,1= 

1,10)/ 




X 1275312., 

1203199., 

1132724*, 

1063444., 

994828. , 

926300* * 



X 

857250., 

786830. 

f 

713H79.* 

636234./ 




DATA { FH2P 

( I ,26), 

I 

= 1 T 1 0 J / 




X 

1352270. , 

1275300. 

t 

1199980 » t 

1125842., 

1052205. , 

978528. , 

X 

904028. , 

827675. 

* 

747830., 

661026./ 




DATA ( FH2P 

(1,27), 

I 

*1* 9)/ 




X 

1428972. , 

1347059. 

t 

1266786., 

1187602., 

1108834., 

1029763,, 

X 

949478. , 

866633. 

* 

778916./ 





DATA ( FN2P 

( I ,28) , 

I 

= 1, 9 }/ 




X 

1505235. , 

1418292. 

, 

1332947., 

1248595. , 

1164471 . , 

1079732 . , 

X 

993247. , 

903208. 

f 

805429./ 





DATA ( FH2P 

( 1,29) , 

I 

= lt 8)/ 




X 

1580909., 

1488816. 

* 

1398277., 

1308589., 

12L8877. , 

1 128133., 

X 

1034872. , 

936456. 

/ 






OATA { FH2P 

(1,30), 

I 

= 1, &)/ 




X 

1655766. , 

1558459. 

* 

1462588., 

1367370. , 

1271798., 

1 1 74615. , 

X 

1073797. , 

964079. 

/ 






DATA < FH2P 

(1,31), 

I 

= 1* 73/ 




X 

1729725,, 

1627072. 

* 

1525687., 

] 424714. , 

1322934. , 

1218714. , 

X 

1108939./ 
DATA ( FH2P 

(1,321, 

l 

= 1, 63/ 




X 

1802543., 

1694420. 

r 

1887375* * 

1480369. , 

1371923., 

1259755./ 


DATA ( F M2 P 

(1,33), 

I 

= 1 , 6}/ 




X 

1874138. , 

1 760389. 

T 

1647441 

1534041 . , 

1418290. , 

1296477./ 


DATA ( 6H2P 

(1,34), 

I 

= 1 i 5 } / 




X 

1944330. , 

1824775. 

t 

1705651 . , 

1585373. , 

1461285./ 



DATA ( FH2P 

(1,35), 

I 

- 1 , 5 3/ 




X 

2012959. , 

1887378. 

t 

1761698. , 

1633861 . , 

1499343./ 



DATA { FH2P (I,3A!,1=1, A i / 

X 2079863., 1947977., 1813243., 1678877./ 

DATA ( FH2P (1,371,1=1, 4! / 

X 2144863., 200629V., 1865786., 1717821./ 

DATA { FH2P (1,381,1=1, 3) / 

X 2207751., 2062003., 1912420./ 

DATA t FH2P (1,391,1=1, 21/ 

X 2268294., 2114594./ 

DATA ( FH2P (1,401,1=1, 21/ 

X 2326179., 2163145./ 

DATA ( FH2P (1,411,1=1, 11/ 

X 2380966./ 

DATA < FH2P ( I ,42 1,1=1, 11/ 

X 2431910./ 

DATA 07 PT/ 4. 8 506, 7.4266, 1 .0359F+1 , 1.362 3 P+1 , 1 . 72206+1,2. 11675+1, 
12.55038+1,3.02 756 + 1 , 3 . 5 5446 + 1 , . 1 372 6 + 1 , 4 . 7830 6 + 1 , 5 , 499 1 6 + 1 , 
26.2936F+1 , 7. 1 7466+1 , 8 . 1 5096+1 , 9. 231 96 + 1 , 1 .04276+2 , 1 .17488+2, 

31 .48096+2, 1 .851 16 + 2 ,2 .29656+2 ,2.82946+2 ,3.46376+2,4.21496+2 , 

45. 1 005F+2 ,6. 1 40 1 6 +2, 7, 3555 6+2, 8. 77128+2,1.1 33 16+3,1.45096+3, 
51.842 86+ 3, 2. 32356+3, 3. 132 06+3, 4.1 768 P + 3 , 5 . 9064F+3 , 8 . 2440P+3 , 

61.1 374E+4, 1 .55206+4,2 .0951 E+4, 2.79888+4,3. 701 1 £+4,4.84606+4, 
78.0735E + 4, 1 . 2969 6+5 , 2 . 0 1 29 6+5 , 3 . 0243 E+ 5 , 4. 4063 £+5 , 6 . 236^F + 5 , 
88.58816+5,1 .15256+6, 1.50916+6,1.93106+6,2.96356+6,4.20186+6, 
95.54326+6,6.83876+6,8,21736+6,9.66106+6 , 1.11 54E+7,1 .26826+7, 

A1 .42336+7 ,1 .57976+7, 1.81476+7,2.04856+7,2.27906+7,7.50526+7, 

B2 .79856+7/ 

DATA TT/ 300. ,400. ,500. , 600 . , 700 . , 8 00 . , 900 . , 1000. , 1100. , 1200. , 

1 1300. . 1400 . . 1500. .1600. . 1700. . 1800. . 1 900 . . 2 000 . . 2 200 . . 2400 . . 2600. 

2.2800.. 3000. . 3200. . 3400. .3600. .3800. .4000. .4300. .4600. .4900. , 


51 



35200. . 5600.. 6000. .6500. .7000. .7500. , 8000 8 500 9000 9500 . » 

410000. . 11000. . 12000.. 13000. . 14000. .15000. . 16000. .17000. . 18000. , 
519000. ,20000. ,22000. ,24000. ,2 6000. ,28000., 30000. , 32000 . , 34 000 . , 

636000.. 38000. .40000.. 43000. .46000 49000 .. 52000 .. 5 6000 . / 

DATA THETAT/ .030 , . 035 , . 040 , . 050 , .060 , .070, .080, . 100 , . 1 30 , . 1 6 , . 20 , 

1.25. . 30. .35, *4, *5,. 6,. 7, *8, 1.0, 1.3, 1.6, 2. 0,2.5, 3. 0,3. 5, 4., 5., 6./, 
2W0T/2.7634E-1, 2. 9930E-1, 3. 2085 E-l, 3. 6077E-1, 3. 9754E- l,4.3202E-l , 
34. 64 77 E- 1,5.2 644 E- 1,6.11L4E-1,6*8566E-1»7.6545£-1,8.315BE-1, 

48 . 6490E— 1 , 8 . 7356 E-l , 8 . 65 30E— 1 , 8 . 2038E-1 , 7 . 5978 F-l , 6 . 9A97E— 1 , 

56. 3748E-1, 5. 346 5E- 1,4.19078-1,3. 3762 E-l, 2. 6264E-1, 2. 01 05 F-l , 

61 . 6003E-1 , 1 .31 19E-1 , 1 . 1004E-1 ,8 . I486 E-2 ,6.34286-2/ , NIT/ 6 .54235-3, 
78.3437E-3»1 .0320E-2, 1 .4798E-2, 1 .9991 E-2 ,2 . 5941E-2 , 3.2704F-2 , 
84.8882E-2 ,7.9674E-2, 1 .1599E-1 , 1 . 66 1 8 E-l , 2 . 2 1 69 5-1 , 2 . 6369E- l , 
92.9218E-1, 3. 0945E-1, 3.20836-1 ,3. 1439E- 1 , 2 . 9978 E~ 1 , 2 . 82 OOE- 1 , 
A2.4567E-1 ,1.9914E— 1 , 1 .6 37 3E - 1 , 1 . 29 58 E-l , i . 00 54E-1 , 8 . 0 733H-2 , 

B6 • 6599E— 2 , 5 • 6 1 1 9E— 2 , 4 . 1 828 E— 2 , 3* 2 698E— 2 / , W2T /4. 68 88 E— 3 , 6 . 05 77 E— 3 , 
C7.5925E— 3,1.1192E-2,1.5575E-2,2.087flE— 2,2.7282F— 2,4.4229E-2, 

D8, 1832 E-2, 1. 32 54E- 1,2. 1022 E-l ,3. 04046-1,3.8063 E-l, 4. 3657F-1, 

64. 7 39 3E-1 ,5. 0763 E-l , 5 .0774E-1 , 4 . 9096E-1 , 4. 6657E-1 ,4.12086-1 , 

F 3. 38 156-1, 2. 8008E- 1,2. 2305 E-l, 1.7391 E-l, 1 .401 1F-1 , 1 . 1586E-1 , 
G9.7787E-2 , 7 .3060E-2 , 5 . 7204E-2 / 

END 
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AIBFTC LMATR DEBUG 

C SUBPROGRAM TO FIND MATRIX AND INITIAL VECTOR FOR ITERATIVE CURRECTIUNS 
C 

c this subprogram finds ist derivatives of pressure with respect to 

C MOLES AND VOLUME AND 1ST AND 2ND DERIVATIVES UF FREE ENERGY WITH 
C RESPECT TO MOLES AND VOLUME. FROM THESE IT FINDS THE 5X5 MATRIX X AMD 
C INITIAL VECTOR Y FOR ITERATIVE CORRECTIONS. 

SUBROUTINE MA TR { V , GMO, S * E T A , C , Q» ROOT * SUMM t SM , X , Y , SN , S V > 
COMMON/BLOCK1/DH,T,PO,PROD,ETAD, SD/BLQCK6/Z 

DOUBLE PRECISION V,GMO( 7) , SN ( 2 > , S V , FN( 7 > ,FV, FVV, FNV ( 7 I , FNN (7, 7 ) , 
1X15,5), YI 5 ) , SM, ROOT, SUMM, PRODO, S, ETA, C, PROODD, RES, W,SN IS V,SN2SV 
DIMENSION PRGO(7 ) ,0(7) ,Z(7) 

SN { 1 ) = 8.3143*T/(P0*V> 

SN( 2>» SN( 1 )-3.4733D3+UH*K00T/ (PO#V) 

SV* -8.3143*T*SUMM/ { P0*V#*2 ) + . 2 7465*DH*T*rtU0 T#*3 / I PQ*V) 

PRODD= 9.91000-2*DSQRT ( ( 1 .264604/T )**3#SM/V ) 

DO 3 L = 1 , 7 

3 FN(L)* -DLOGI V*Q( L ) / ( PROD! L )*GMO< L ))) + l.-Z < L )*PP.ODD*DH 
DEBUG SN( 1 ) ,SN(2) ,SV, ( FN{ I) , 1 = 1 ,7) 

F V = -SUMM/V + .O33D33#RO0T##3#DH 

FVV=5UMM/V**2-.049550*RUnT**3*DH/V 

PRODDD® .049350*1 1 .264604/ ( f*V ) >**1 . 5D0+0S0RT( SM) 

DO 4 M=l,7 

4 FN V ( M ) = -l./V+Z (M)#PRODOD*DH 

R ES= . 049 550*DSQR T ( ( 1 . 2 64604/ 7 ) **3/ ( V*SM) ) 

DEBUG FV, FVV, (FNV(I), 1=1,7) 

DO 5 1=1,7 
DO 5 N= 1,7 
IF ( I -N ) 7,6,7 

A FNN ( I , N ) = l./GMO( I )-Z ( I)*RES*DH 
GO TO 5 

7 FNN ( I , N ) = -Z ( I )*Z [NJ#RES*OH 

5 CONTINUE 

W= ( SD-S+C*SN<2 )-.500*SN( 1 ) *< E TAD-ETA ) ) / S V 
SN1 SV= SN ( 1 ) /SV 
SN2 S V= SNI21/SV 

DEBUG W,SN1SV, SN2SV , ( ( FNN ( I , N ) , N=I,7),I = 1,7) 

X( 1 , 1 ) = ( FNN ( 1 , 1 )+. 2 500* FNN (4, 4 )+SNlSV#( . 2 5D0*SN1 S V*F V V-FNV ( 1 >+.500 
1*FNV (4 ) ) —FNN ( 1,4)) *GMO< 1 ) 

X( 1 ,2 ) = ( ,2aD 0*FNN ( 4,4)+SNlSV* I ( — . 2 5D0*SN 1 S V + SN2 S V ) *F V V- . 5D0* F N V ( 2 i 

1 — • 50O*FN V ( 3 ) + .5D0*FNV( 1 ) ) -2 . *SN2 S V* ( FNV I 1 ) - . 5 DO*FNV ( 4 ) )+FNN< 1 ,2 > + 

2 FNN ( 1 , 3 ) - . 5D0*FNN f 1 ,4 )-. 5lH.*FMN< 2 ,4 ) 5D0*FNN( 3,4) )*GM0(2) 

XI 1 ,3) = { ,25D0*FNN (4,4 )+SNl S V* t - . 25D0*SN1 S V*F V V+ . 5D04 FNV I 3 J-.5D0# 
1FNVI 5) + .5D0*FNV( 1 ) » — FNN ( 1 , 3 ) - . 5D0*FMN t 1,4) + FNN ( 1 , 5 ) + . 5 DO* FNN ( 3 , 4 ) - 
2.5 DO*FNN 14,5) )*GM0(5) 

X I 1 , 4 ) = ( , 5 DO* FNN ( 4, 4 ) +SN 1 S V* I (SN2SV-. 5D0*SN 1 S V ) *F V V- . 5D0* FNV (3)- 
1 .5D0*FNV(6)+FNV (1> )-2.*SN2SV* ( FNVI 1 )-.500*FNV 14) )+FNN( 1 ,3)-FNN( 1 ,4 
2 ) + FNN( 1 »6}-.5D0*FNN( 3,4)-.500*FNN(<i,6) )*GMil( 6) 

XI 1 ,5) = < .7 500* FNN (4»4)+SMlSV*( I - . 75 00*SN 1 S V+ SN2 S V > * FVV- . 5U0*FNV (3) 
1+1 .5D0*FNV{ 1 )-. 5D0*FNV { 7 )) +SN2SV* ( -2.*FNV ( 1 )+FNV ( *, I ) + FNNI 1 , 3 )-l. 3D 
20*FNN<1 ,4)+FNNI 1 ,7 ) - . 5D0*FNN ( 3 , 4 ) -. 5 DO*FNN I 4 , 7 > )#GMI’J( 7) 
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Y ( 1 ) =- ( F N ( 1 >-.5D0*FN<4 >-. 5DQ*SN 1 S V*F V+ ( ETAD-ETA ) * ( -,25D0*FniN ( 4 , 4 )- 

1 ,25DO*SN1SV*FNV(4)+.5DO*FNN( 1,4) ) -W * .5 DO*SN 1 S V* FVV+U* < FNV ( 1 )-.50o* 
2FNV<4) ) +C*( ,5D0*SN1SV*FNV( 3)-FNM{ 1 , 3 ) +. 5D0*FNN ( 3 , 4 ) ) i 

X ( 2 , 1 ) = ( • 2 5 DO* FNIM ( 4 , 4 ) +S N 1 S V * ( ( 25D0*SN1S V+SN2S V ) * F V V — • 5D0*FNV (2 ) 
l-.5DO*FNV( 3) + .5DO*FNVU ) ) +SN2 SV* < -2 . *FNV ( 1 )+FNV<4) )+FNN< 1,2)+ 

2FNN ( 1 ,3)-.5D0*FNNU ,4)-.5T)0*FNN ( 2 ,4) 5DO*FNN ( 3 ,4) )*GMO< 1 ) 

X(2,2)=( FNM (2,2 ) +FNN ( 3 , 3 ) + . 2 5DO*FNN ( 4, 4 ) + ( . 5D0* SN1 SV-2 . *SN2 S V ) **2* 
1FVV+(SN1SV-4.*SN2SV )*( FNV( 2J+FNV ( 3 ) 5DO*F NV ( 4 ) > +2 . *FNN ( 2 , 3 ) - 

2 FNN (2,4 ) — FNN (3,4) >*GMO( 2 ) 

X (2 ,3)=l-FNNt3,3)+.25D0*FNN (4,4 )+SNlSV*( l .2 5DO*SN 1 SV-SN2S V }*FVV- 
1 .5D0*FNV(4)+.5D0*FNV<2 )+.5D0*FNV ( 5 ) ) +SN2 SV* ( 2 .* FNV (3>--2.*FNV(5) + 

2 FNV { 4 > ) —FNN (2,3)— , 5 00* FNN (2,4) + FNN' (2,5) + PNN! 3 ,5 ) -,5C0*FNN ( 4,5 ) >* 
3GM0 1 5 ) 

X ( 2 , 4 ) = ( F NN ( 3 , 3 ) + • 5 DO* F N N ( 4 , 4 ) + ( . 5DQ*SN 1 SV-2 ,*$N2 SV > * ( SN1SV-2.* 
1SN2SV)*FVV+SN1SV*( 1.5D0*FNV (3)-FMV(4)+FNV( 2 ) + . 6DG* FNV ( 4) ) + SM2SV* ! - 
24.*FNV ( 3>+3.*FNV (4)-2 ,*FNV( 2 )-2 ,*FNV(6) ) + FNN (2,3) —FNN (2,4) + F NN (2,6 

3) -l. 5D0*FNN( 3,4 > + FNN( 3,6 >-.5D0*FNN 14,6 ) )*GMH( 6) 

X{2,5)=( FNN I 3, 3> + .75D0*FNN(4,4) + < , 5D0*SNlSV-2 , *5 N2 S V ) * ( 1 ,5D0*SN1SV 
1— 2.*SN2SV )*FVV+SN1SV*( 2 . *FNV (3 ) — 1 .500 *F MV (41+1,5 DO* FNV ! 2 )+.5D0* 

2 FNV ( 7 ) )+SN2SV*(-4.*FNV( 3 ) +4 ,*F NV l 4 ) -2 . *FNV ( 2 ) -2 . *FNV ( 7) )+FNN(?,3)- 
31, 5D0*FNN(2,4)+FNN( 2,7 )-2 . *FNN ( 3 , 4 ) ~ . 5D0*FNN ( 4, 7 ) + FNFK 3,7) )*GMO( 7 i 

Y (2 ) =-( FN(2 ) + FN( 3 )-. 5D0*FM< 4 ) +F V* ( ,5D0*SN] S V-2 . *SN2 SV ) +C* ( -FNN (3,3 
1 )-( . 5D0*$N1 SV-2 ,*SN?SV )*FNV ( 3) -FNN (2, 3 ) + . 5 DO* FNN (3,4) ) + ( FTAD-ETA1* 
2 ( —.2 5D0*FNN ( 4 ,4 } + ( . 2500*SM1 SV -SN2SV ) *FNV ( 4 > + . 5D0*FNN ( 2 , 4 } + .500* 

3 F NN (3,4) ) + W * ( ,500*SN1SV-2.*SN2SV) *FVV+W* ( FNV ( 2 )+FNV{ 3)-. 500* + NV (4) 

4) ) 

X( 3,1 > = ( . 2 5D0*F NN ( 4 , 4 ) +SN 1 S V* ( - . 2 5 00*5N 1 S V* FV V+ . 5D0*FN V ( 3 ) -.500* 

1 FNV ( 5 ) + . 5D0*FNV ( 1 ! ) - 1- i N l 1 , 3 ) - . 5D0*FNN ( 1 ,4) + FNN( 1 , 5 ) + . 5D0*FNN ( 3 , 4 ) - 
2,5 DO* FNN ( 4,5} ) *GMl)( 1 ) 

X ( 3,2 ) = (-FMN1 3,3)+.25Q0*FNN ( 4 , 4 ) +SN 1 S''* ( ( , 2 5D0*SN 1 S V-SN2S V )*PVV- 
1 . 5D0*FN V ( 4) + ,5D0*FNV (2 )+.5D0*FNV ( 5 } )+SN2SV*( 2.*FNV ( 3 >-2.*FNV ( 6 ) + 
2FNVI4 } ) — F NN (2,3) -. 5 DO*FNN ( 2 , 4 ) + FNN ( 2 , 5) + FNN( 3 , 5 )-,5D0*FNN (4,5))* 
3GH0 ( 2) 

X (3, 3) = t FNN ( 3,3) + , 2 5 DO* FNN (4,4) +FNN ( 5,5 )+SNlSV* ( . 2 5Q0*SN 1 S V*F V V- 

1 F NV ( 3 ) +F NV ( 5 ) — . 5D0*FNV 14) ) +F NN ( 3 , 4 ) — 2 » *F NN (3,5) —FNN (4,5) )*f,ND< 5 > 

X< 3 , 4 ) = ( — F NN ( 3 , 3 ) + . 5 00* FNN (4,4 ) +SN1SV* ( .500* ( SN IS V-2 ,*SN2 SV )*FVV 

1-. 5D0*FNV ( 3)+.5DO*FNV ( 6 )-FNV (4 J + FNV ( 5 ) ) +SN2 SV* ( 2 . * FNV { 3) -2.* 

2 FNV ( 51+FNV (4) ) + .5D0*FNN ( 3 ,4 ) -F NN ( 4 , 5 ) -. 5D0*FNN ( 4 , 6 ) +FM N (5,6) + 

3FNN1 3,5 )-FNN( 3,6 ) ) *GMU(6 ) 

X ( 3 , 5 ) = ( - FNN (3,3) + , 7 6 1)0* FNN (4,4)+SNlSV*( ! . 75D0*SN 1 S V-SN2 S V ) *1" V V - 
1 F NV ( 3 ) — 1 , 5D0*F NV ( 4 ) + 1 . 5D0*FNV ( 5) +.5 00*FNV ( 7 > ) +S N2S V* ( 2 .* FNV ( 3)-2. 
2*FNV( 5 ) + FNV(4) } + FNN ( 3,4)+FNN< 3 , 5 ) -1 . 500* FNN ( 4 , 5 ) - , 5D0*F NN ( 4 , 7 ) + 
3FNN( 5, 7 )-FNNl 3,7 ) ) *GMD ( 7 ) 

Y ( 3 ) =- ( — F N ( 3 ) — . 5 DO* F N ( 4 ) + F N ( 5 ) + . 600* S N 1 S V * F V+ C * ( FNN ( 3, 3 }-.5D0* 

1 SN 1 S V*FNV ( 3 ) + , 5 DO* FNM ( 3 ,4)-FNN<3, 5 ) ) + ( ETAD-ETA ) * < - . 2 500* FNN ( * , 4 ) + 

2 ,25D0*SN1SV*FNV(4)-.5D0*FNN i 3 , 4 ] + . 500* FNN < 4 , 5 ) ) + N*. 500* SN 1 SV*FW+ 
3W*( FN V ( 5 ) — FN V (3) — .500* FNV (4 ) ) ) 

X (4, 1 } = ( . 50 0 * F NN (4,4 ) +S N 1 S V* ( ( -. 500* SN1 SV+SN2 SV ) *F VV- . 500* FNV ( 3 )- 
1 ,5D0*FNV ! 6 ) +FN V ! 1 ) ) -2 - *$N2 S V* ( FNV ( 1 )-.5D0*FNV (4 ) ) + FNN( 1 ,3 }-FMN( 1,4 

2 ) +FNN ( 1,6)— . 5 DO* FNN ( 3 , 4 ) — . 5 DO*FNN ( 4 , 6 ) )*G/K) ( 1 ) 

X( 4,2 )= ( FNN( 3,3 )+. 500*FNN ( 4 , 4 ) + ( . 5 DO*SN 1 SV-2 . *SM2 S V ) * ( SMI SV-2.* 
1SN2SV ) *FVV+SN1SV* ( 1 ,5D0*FNV ( 3 ) - FNV ( 4 ) + FNV ( 2 ) + . 500* FNV ( 6 ) 1 + SN2SV* 

2 ( -4 , * FN V ( 3 ) +3 . * F N V ( 4 ) - 2 . * F N V ( 2 )-2 ,*FHV (6) )+FNN( 2,3) — FNN (2,4) + 

3 FNN (2,6) — 1. 500*FNN ( 3 , 4 ) + F NN ( 3 , 6 ) - . i>00*FNN ( 4 , 6 ) )*(,MIJ ( 2 ) 

X ( 4 » 3 ) = ( — F NN { 3 , 3 ) +. 5 DO* FNN (4,4)+SNlSV*( { . 5D0*SN1 S V-S M2 S V ) *F V V + 

1 FNV ( 5 ) + . 5 DO* FNV (3 ) — . 6D0*FNV (4 ) + . 500*F N V ( 6 ) — FNV ( 3)-. 5D0*FNV ( 4 ) ) + 
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2SN2SV*(2.*FNV< 3 )-2.*FNV ( 5 ) +FNV < 4 > ) +. 5D0*FNN ( 3, 4 l-FNN ( 4 i 5 5D0* 
3FNN< 4,6 l+FNN < 5,6 )+FNN( 3, 5) -FNN <3, 6} )*GM0(5) 

X(4,4) = <FNM3,3 ) + FNN (4,4) +FNN (6, 6 ) + ( SIUS V-2 . *SN2 S V > **2*F VV+2 .* ( 
ISN1 SV-2 .*SN2 SV > * ( FNV ( 3 ) +FNV l 6 )-FNV l 4 > ) -2 . *FNN < 3 , 4 ) + 2 . *FNM I 3 , 6 > -?. . 
2*FNN<4,6) >*GM0<6) 

X(4,5) = (FNN(3,3) + 1.500*FNN(4,4} + <SN1SV-2.*SI"2SV)*< 1 . 5DQ*SN 1 SV-2 .* 
1SN2SV)*FVV+SN1SV*( 1 . 5D0*FNV ( 6 )-3 . *FN V < 4 ) +2 . 51)0*1= FNV ( 3)+FNV( 7) ) + 
2SN2SV* (—4.* FNV ( 3 ) + 5.*FNV (4 ) -2 .*FNV ( 6 ) -2 . *F NV { 7 > > -2 . 50O*FNN ( 3 , 4 > + 
3FNNI 3,6)-1.5D0*FNN(4,6)-FNM4,7)+FNI'l< 6, 7 )+FNN( 3,7 > )*GM0( 7) 
Y(4)=-(FN<3)-FN<4)+PN(6)+FV*< SN 1 S V-2 .*SN2 S V ) +C* ( -FNN ( 3 , 3 ) -FNV ( 3 ) * 

1 ( SNISV— 2 »*SN2SV )+FNN( 3,4 )— PNN (3,6) ) + ( ET AD-ET A ) * ( 5D0* FNN { 4 , 4 ) + 
2.5DO*FNV(4>*(SN1SV-?.*SN2SV) + .5DO*FNM3,4) + .5DQ*FNM<4,6 ) )+W*{ 

3SNL SV— 2 »*SN2S V > *FVV+W* ( FNV ( 3 ) — FNV { 4 > + FNV ( 6 ) ) > 

X< 5,1) = < .75(J0*FNN (4,4)+SNlSV*( < 75D0*SN1 S V+SN2 S V ) *F V V- . 5D0* 

IFNV(3)-.5D0*FNV<7)+1.500*FNV(1 ) ) + SN2 SV* ( -2 . *FNV ( 1 )+FNV(4 ) } + 

2FNN< 1 , 3 ) -1 . 500* FNN (1,4) + FNN ( 1 , 7 )-.5D0*FMN ( 3 r £, ) 5D0 *Fp'IN ( 4 , 7 )> * 
3GMO( 1 ) 

X (5,2)<*(FNN (3,3)+. 7500* FNN (4, 4) + ( . 500*SN L S V-2 .*SN2 SV ) * < 1.500* 
1SN1SV-2.*SN2SV)*FVV+SN1SV*(2.*FNV< 3)-l .500* FNV (4)+1.5D0*FNV(2)+ 

2 » 5D0*FNV ( 7 ) >+SN2SV*< -4.*FNV( 3)+4.*l-NV (4)-2 ,*FNV < 7 ) -2 . *F NV ( 2 ) ) + 
3FNN( 2,3>-1.5D0*FNN(2,4)+PNN< 2,7 )-2.*FNN ( 3,4 ) - . 5U0* FNN < 4 , 7) + 

4FNN( 3,7) ) *GMO< 2 ) 

X < 5,3 ) = <-FNN( 3,3 ) +. 75D0*FHN ( 4,4 ) +SN1 SV * ( ( .75t)0*SNl SV-SN? SV ) *F VV + 

1 l .5DO*FNV ( 5 ) -FNV < 3 ) — 1 »5I)0*FNV (4)+.5D0*FNV ( 7 ) ) +SN2 SV* ( 2 .*FNV ( 3 ) -2 . 
?*FNV( 5)+FNV<4) ) +FNN ( 3,4 ) +FNN ( 3 »5 ) —1 • 500* FNN ( 4, 5 ) - . 500* FNM 4 , 7 ) + 

3 FNN (5,7) -FNN (3,7)) *GM() ( 5 ) 

X( 5,4)»(FNN(3,3 )+1.5D0*FNN(4,4)+( 5N1SV-2.*SN2SV)*( 1 . 5D0*SNlSV-2.* 
1SN2SV )*FVV+SN1SV*< 1 .500* FNV < 6 ) -3 .*F NV ( 4 ) +2 . 500* FNV ( 3 ) + FNV ( 7 ) ) + 
2SN2 SV* ( -4. *FNV ( 3 ) + 5.*FNV (4 )-2.*TNV ( h ) -2 . *F!m V { 7 ) )-2.5D0*FNN( 3 ,4 ) + 
3FNN( 3 , 6 ) — 1 . 500* FNN (4,6) -FNN < 4, 7 ) +FMN ( 6, 7 ) + FNl'l ( 3 , 7 ) ) *GM0 ( 6 ) 

X ( 5, 5) = ( FNN ( 3, 3 ) + 2.2*>D0*FNN< 4,4) + FNN< 7, 7 ) + ( 1 . 50O*SN 1 S V-2 . *$ N2 S V > * 
1*2*FVV+ ( 3.*SN1SV-4.*SN2SV )*< FNV ( 3 ) + FNV ( 7 ) - 1 . 5D0*FNV ( 4 ) > -3 . * 

2FNN( 3 ,4)— 3.* FNN (4 , 7 ) +2 »*FNN (3,7) ) *GMC) ( 7) 

Y ( 5 ) =-( FN ( 3 ) -l . 5 DO * F N ( 4 ) + FN ( 7 ) + FV*( 1 . 50C)*SN 1 S V-2 . *SN?S V ) +C*(- 
1 FNN (3, 3 )-( 1 .500*SN1SV-2.*SM?SV ) *FNV ( 3 ) -FNN( 3 , 7 ) +1 . 500* f NN< 3,4))+ 

2 ( E T A 0-E T A ) * ( — . 7 5 00* F NN ( 4 » 4 ) + ( .75*>0*SN1 SV-SN2SV ) *FNV ( 4 ) ■(■ . 51)0* 

3FNN< 3,4) + . 5l)0*FNN (4,7) )+W*( l .SDO-SHl S V-2 . *SN?S V ) *F V V+W* ( FN V ( 7 ) 
4+FNV ( 3 ) —1 . 5t)0*FNV ( 4 ) ) ) 

DEBUG ( (X( I , J) , J = 1 ,5) , 1 = 1,5 ) 

DEBUG ( Y ( I ) , I =1 , 5 ) 

RETURN 

END 
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$ I8FTC LBINT 

C SUBPROGRAM FOR 4 POINT INTERPOLATION 
C 

C THIS SUBPROGRAM PERFORMS 4-POINT INTERPOLATION BY THE METHOD OF 
C LAGRANGE. X AND Y ARE THE INDEPENDENT AND DEPENDENT VARIABLES, 

C RESPECTIVELY, IN A TABLE OF N ENTRIES. Z IS THE VALUE OF X DESIRED. 

C BINT IS THE INTERPOLATED VALUE OF Y. 

FUNCTION B INTI X,Y,Z ,N3 
DIMENSION X(67),Y(673 
Ml = I 
M2-N 

5 I F ( M2— Ml — 1 >7,7,6 

6 MT£M=M 1 + I NT ( .5*1 FLOAT t M2-M1 ) + .0001 3 ) 

IF( Z-X(MTEM) 18,1,9 

8 M2 = MT EM 
GO TO S 

9 Ml =MT EM 
GO TO 5 

1 B I NT = Y ( M T E M ) 

RETURN 

7 IFIMl— 1 32,2,3 

2 I -2 

GO TO 4 

3 IFIM2-N til, 10*10 

10 I*N~2 
GO TO 4 

11 1=M1 

4 01=Z-X{ 1-1 ) 

02 = Z— X ( I > 

03=Z-XI 1 + 1 ) 

04=Z— X ( 1+2 I 
012=01*02 
034=03*04 

P21 = X( I 3 — X ( 1-1 ) 

P32=Xf 1+1 )-X( I ) 

P42 = X ( I + 2 > —X ( I ) 

P31 =X ( I + I ) — X ( 1-1 ) 

P41 = X { I + 2 > —X ( 1-1) 

P43 = X I 1+2 > — X < 1 + 1 } 

BINT=034*(01*Y( I 3 / ( P32*P42 > -02* Y ( i -1 > / f P3 1 * P4 1 3 3 / P2 1+012* < Q3*Y ( 1+2 
1 3 / ( P41*P42 )-04*YU + l 3 / ( P31*P32 3 ) / P43 
RETURN 
END 
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SIBFTC LEXPC 

C SUBPROGRAM FDR ANTICIPATING AND PREVENTING POSSIBLE UNDERFLOWS 
C 

C THIS SUBPROGRAM SETS EXPONENTIALS EQUAL TO 0 IF THEIR ARGUMENT IS 
C SMALLER THAN -72. 

FUNCTION E X PC ( A ) 

2 IF< A+72.)3,4,4 

3 EX PC=0 • 

GO TO 5 

A EX PC = £XP (A) 

5 RETURN 
END 
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$ I BFTC LINV 

C SUBPROGRAM FOR INVERTING A 5X5 MATRIX 
C 

C THIS SUBPROGRAM INVERTS A MATRIX BY THE METHOD OF CRUUT. THE ELEMENTS 
C ARE NOT REARRANGED. THE DIAGONAL MAY HAVE NO ZERO ELEMENTS. THE MATRIX 
C IS ASSUMED TO BE STORED IN REGION AA. THE INVERSE IS STORED IN THE 
C SQUARE SUBARRAY At 1,6) THROUGH A(5,10). 

SUBROUTINE INVERT ( A A , A ) 

DOUBLE PRECISION A < 5 , 11 ) , A A { 5 , 5 ) ♦ SUM 

M AG=5 

BIG = 1.1 

SMAL = 0.9 

L 1 =MAG+ 1 

L2=2*MAG 

L3 =L 2 + 1 

LA=MAG-1 

C SAVE MATRIX AND AUGMENT MATRIX. 

101 DO 1003 E1=1,MAG 

102 DU 1001 J 1 = L 1 »L3 

1001 At 1 1 , Jl> = 0.0 

103 OU 1002 J 2 = 1 , MAG 

A < IL,J2)=AA< II, J2) 

1002 A( II ,L3 I=AII1 ,L3)+A( I 1 , J2 ) 

A ( II ,L3 > = A ( 1 1 ,L3 > + l .0 

J 3- 1 1 + MAG 
A ( I ]., J 3 > = 1 . 0 

1003 CONTINUE 

C BEGIN FORWARD SULJTIIJN. 

201 DO 2001 J4=2,L3 

2001 At 1 ,J4)=A(1,J4)/At 1, 1) 

LORA =2 

2002 LORI = LOR A- 1 
L 0R2=L OR A+ 1 

202 DO 2003 I 2 = LOR A , M AG 

203 DO 2003 K1=1,L0R1 

2 00 3 At I 2 , L OR A )=A( I 2 T L UR A ) - A ( I 2 , K 1 ) * A ( K 1 , L UR A ) 

204 DU 2005 J5=L0R2,L3 
SOM = 0.0 

205 DO 2004 K2=l,L0Rl 

2 004 SUM=SOM+A ( L UR A , K2 ) #A ( K 2 » J 5 ) 

2 00 5 A t LORA, J5>= l A t LOR A , JS ) -SUM) / A t LORA , LIJRA ) 

I F { MAG— LOR A ) 999,901, 999 
999 LURA = LURA+1 
GO TO 2002 

C BEGIN BACK SOLUTION. 

901 DU 9001 N2=1,L4 
1 4=MAG— N2 

402 DO 4001 J 7 = L 1 ,L3 

403 DO 4001 KSUM= 1 ,N2 
K 3 = L 1— K SUM 

At I 4 , J 7 ) = A ( I 4 , J 7 ) - A ( I4,K3)*AtK3,J7) 
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4001 CONTINUE 
10000 RETURN 
END 
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APPENDIX C 

CUTOFF OF H AND H 2 PARTITION FUNCTIONS 

If the energies of the excited electronic states of H and Hg were calculated considering 
all neighboring particles, there would probably be a finite number of bound states; there- 
fore, no difficulty would arise in evaluating the partition functions. Unfortunately, no 
rational experimentally verified method has been developed for calculating these perturbed 
energy levels for conditions in this report, and in most cases the perturbations caused by 
the neighboring free particles are small (ref. 10). Consequently, it is customary to use 
the unperturbed energy levels (or term values) and to cutoff the series in the partition 
function at some reasonable energy level to prevent the series from diverging to infinity, 

Cox (ref. 35) has pointed out that for high degrees of ionization the perturbation of the 
energy levels is due principally to coulomb forces , so that one method of cutoff should be 
used, whereas for low degrees of ionization the perturbation of the energy levels is due 
principally to neutral particles, so that another method of cutoff should be used. In this 
report the two methods used were the Debye -Hiickel cutoff and a modification of the Bethe 
cutoff (ref. 36). The versions of these two cutoffs used herein are discussed in the fol- 
lowing paragraphs. 

In the Debye- Hiickel cutoff (ref. 10) states are not counted in the partition function if 
their unperturbed ionization energy is less than the absolute value of the lowering of the 
ionization potential A I, which is given by (ref. 34) 

2/2 \l /2 

AI = — — f— — E z 2 N. ] (Cl) 

4 I rhce o \€ o kT i V 

Here the summation is over all species. 

For hydrogen atoms, strict adherence to this cutoff results in large discontinuities 
in the partition function as the charged- particle density is varied at high temperatures* 

This can cause serious difficulties in numerical methods for finding equilibrium and is 
physically unrealistic because states near the cutoff are greatly broadened* A more 
realistic approach is to cutoff each state gradually. In this report this gradual cutoff is 
accomplished by essentially the same method used by Krascella (ref. 1). Let F r be the 
fraction of the contribution of the unperturbed state to the divergent partition function to 
be counted in the cutoff partition function. Then 
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F n = 1 


F = 
n 


Vl - *n 


^ + ^ n+1 - - AI ^ 0 j 


o\ 


In- AI i 

F n = — + — 

V*n + 1 2 


2 

*n + *n+l 
2 

and 


- Al< 


A - 4151 


^ - 441 i 


+ — 



(C2) 


F n = ° 


AI - 


Vl + I n 


°) 


where I n _jt I n , and I n+1 are all for the same Rydberg series. (If the n - 1 state did 
not exist, a weighted average of the ^ j for Rydberg series that did exist was used.) 

To summarize equation (C2) (1) if AI was less than (1^ + I n+1 )/2, there was no cutoff; 

(2) if AI equaled ^ the state was half cut off; (3) if AI was greater than (I n _j + 1^/2, 
the state was completely cutoff; and (4) for intermediate values of AI, the cutoff was 
linear in AI. 

In the original Bethe cutoff (ref. 36), all states were cutoff whose highly excited 
electron had a radius a n greater than half the average distance between neighboring 
molecules. This half distance r Q was taken as the radius of a sphere whose volume was 
equal to the gas volume per molecule. The radius a n was given by 


a n = a o n 


(C3) 


This cutoff was used by Olsen (ref. 37). A modified Bethe cutoff with the cutoff accom- 
plished gradually in a manner linear in n was used by Bond (ref. 38) and Pike (ref. 39). 

Equation (C3) is realistic for H atoms. However, for H 2 , a more realistic approxi- 
mation was made. For Rydberg atomic orbitals of an atom 
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The equations for F for the modified Bethe cutoff were obtained by replacing I ,, 
n 1/2 1/2 1/2 1/2 n- l 

1^, and AI in equation (C2) by ^ , -a^ 7 , -a n ^j > and -r Q 7 , respectively. 

Thus, for H atoms F was the same as in the Bond modification (ref. 38) to the Bethe 


cutoff . 


In computing and q^, F n was calculated for each state by both the Debye -Huckel 
cutoff and the modified Bethe cutoff, and the smaller value was used. Hence, for low 
degrees of ionization the Bethe cutoff was always effective. At high degrees of ionization 
the method of cutoff depended on the conditions, and no general statement can be made. 
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APPENDIX D 


METHODS OF OBTAINING CONSTANTS OF EXCITED ELECTRONIC STATES H ? 

For all excited Rydberg states Hg, Hund's coupling case (b) (refs. 17 and 20) was 
assumed to hold. For each Rydberg state, a Morse potential (ref. 17) was used as a 
model to calculate the partition function. The formula for a Morse potential is 

„ r -/3(r-r)T 2 

U(r) = T + D 1 - e e (Dl) 

The vibrational term values T + + G (see fig. 8) of Morse oscillator are (ref. 17, p. 101) 


T + + G y = T + jS' Dh 


k l/2 


27T 2 C/M/ 




(D2) 


measured from the rotationless ground vibrational state of Hg. If 

Dh N 1 / 2 


3 fit 


(,2tt 2 c/u, 


(D3) 


equation (D2) becomes 


T + G y = T + uj 


H)-^K ) 2 


(D4) 


For electronic states with no hump in their potentials (see lower Hg potential in 
fig. 8) and known vibrational term values, the appropriate method of fitting the Morse 
potential to the actual potential has been given in reference 18. It includes making W 
for the Morse potential the same as the actual W, where W for these states is the 
energy difference between the dissociation asymptote and the rotationless ground vibra- 
tional level of the state (see fig. 8). From equation (D4) 


D = 


W +~+ Vw(W + O)) 

2 


<D5) 


The constant tu was determined by fitting the excited vibrational energy levels (ref. 18). 
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In general, D for the Morse potential was not the same as the actual well depth, and to 

determine it a knowledge of the actual dissociation asymptote was, of course, required. 

* 

As shown in figure 8, the potential curves of some excited electronic states B 2 

possess a hump, which presented a problem because the Morse potential has no hump. 

However, the principal effect of the decrease in the potential energy to the right of the 

hump is to broaden those bound energy levels lying above the dissociation asymptote but 

not shift them much. Hence, a good approximation was to set W for the Morse potential 

equal to the energy difference between the top of the hump and the rotationless ground 

* 

vibrational state (see upper Hg potential in fig. 8), The Morse well depth is then given by 
equation (D5)* A knowledge of the actual dissociation asymptote and the hump height H 
(if any) are required to find W and D. 

For all excited electronic states of Hg In this report, the dissociation products are 
the atoms H(l) and H(n), where the quantity in parentheses is the principal quantum 
number* The energy of the dissociation asymptote is thus determined by the value of n. 
The state H(n) was determined for each electronic state of H 2 by the correlation rules 
(ref. 17) also known as noncrossing rules* (Note that fig. 157 of ref, 17 does not apply 
to Hg despite the label Hg on it. ) The values of n in H(n) are tabulated in table n 
(pp. 76 and 77), 

* 

Some of the hump heights of H 9 were obtained from the literature, whereas others 
had to be predicted. Six electronic states H 2 have been reported (refs. 15, 41 , and 42) 
to possess humps, and the heights have been given. Of these, five have promoted molec- 
ular orbitals (MO). The sixth has a much smaller hump than the other five. On the 
other hand, two states with promoted molecular orbitals but no humps have been reported 
(ref. 41). Unfortunately, there are many electronic states H 2 for which there are no 
reports on the presence or absence of humps. For these it was assumed that there was 
no hump if there was no promoted MO* For states with promoted molecular orbitals, 
an approximation was made* The ratios H/6 (fig. 9) were calculated for the seven elec- 
tronic states with promoted molecular orbitals for which values of H were available 
(including the two states with H = 0). An average of the ratios was then calculated by 
using the electronic statistical weights g^ as weighting factors. This average had the 
value 0. 15036. It was assumed that all other states with promoted molecular orbitals 
had H/6 - 0* 15036. Values of H are tabulated in table H. 

It was necessary to have Morse potentials for every Rydberg state H 2 with n of the 
excited electron of the molecule having values from 2 to 17. A different Morse potential 
was used for each electronic state with n from 2 to 4, and experimental term values 
were used when available in determining the Morse potential. For large n all states 

+ 1 

Hg have potentials nearly the same as H 2 , and for n greater than 4 there is little exper- 
imental information. Hence, for electronic states with n from 5 to 17, the same Morse 
potential was used for all states with a given value of n. The constants for Morse 
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potentials with n from 5 to 17 were estimated from the constants for Hg and the averages 
of the constants for the n = 3 states. (Constants for many of the n = 4 states were 
unreliable or had to be estimated.) The most used methods of determining or estimating 
Morse constants for n = 2 to 17 are given below. Less used methods are described in 
footnotes in table 33. 

Morse potentials were most easily fit to states for which the rotationless ground 
vibrational term value, one or more excited pure (rotationless) vibrational term values, 
and several rotational term values were available. Setting v and G v equal to 0 in 
equation (D4) gives 


T + - T = - « - (D6) 

2 16D 

Substituting equation (D6) into equation (D4) 

G„ = wv - ~ v(v + 1) (D7) 

v 4D 

Equation (D5) may be substituted into equation (D7) to get G y in terms of w, v, and W. 

In accordance with reference 18, a> was then varied to minimize the sum of the squares 
of the deviations of the calculated quantities G y from the experimental quantities G v - 
The well depth D was then found from equation (D5), and ft was found from equation (D3). 
The last step was to find r . To do this, the Morse potential was expanded in powers of 
(r - r a )/r Q so the rotational term values f(v,k) for each vibrational level could be 
expressed by Dunham* s formulae (ref, 43), The equilibrium intemuciear distance 
was then varied to minimize the sum of the squares of the deviations of the calculated 
quantities f from the experimental quantities f, 

A more difficult problem was presented by electronic states for which no rotationless 
ground vibrational term values were known , but for which a number o! vibration* rotation 
term values T + + + f were known for more than one value of v. An expression was 

derived from equation (D2) and Dunham’s formulae which gave vibration* rotation term 
values in terms of the energy of the dissociation asymptote , the hump height , D, r 0j 
v, and K, The quantities D, 0, and r e were then simultaneously varied to minimize 
the sum of the squares of the deviations of the calculated term values from the experi- 
mental term values , using weighting factors. To be as nearly consistent with refer- 
ence 18 as possible 5 each vibrational level for which any term values were available was 
given the same weight. 

For some n = 4 states there was only one vibrational term value or none at all 
available. For such cases there is no good way to estimate ^ or w. However, fair 
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estimates were obtained by interpolation of f} along a Rydberg series with as its limit. 

The interpolation was accomplished linearly between n = 3 and n = <» (i.e. , Hit) in the 

ft - 1/n plane. When no n = 3 state of the Rydberg series existed , fi(n = 3) was taken 

as the average for all n - 3 states, using g as a weighting factor. In most cases, 

e 

was then found from equation (D3). 

For all states with n = 5 to 17, all Morse constants were estimated. First, the 
average I (fig. 8) for all n = 3 states was calculated, using the g as weighting 
factors. From this value the average quantum defect (ref. 40) was calculated. It was 
then assumed that the average quantum defect was independent of n (which is not strictly 
correct because the quantum defect would be expected to be independent of n only for 
large n along a Rydberg series). Thus, average values of I and, hence, T were cal- 
culated for n = 5 to 17. An average value of p was found for each value of n by inter- 
polation, as described in the preceding paragraph. For each n, an average value of D 
was found from the average T, the various dissociation asymptotes, and the various 
hump heights, using the g g of the H 2 states as weighting factors. 

For some n = 4 states and all states with n = 5 to 17, it was necessary to estimate 
r . Herzberg (ref. 17, p. 456) gives the empirical relation 

r^u' = Constant (D8) 

for electronic states of the same molecule. The extension of equation (D8) to the molec- 
ular ion resulting from removing an electron from that molecule is obvious. Many of the 
values of r c for Hg are unreliable, so for this report equation (D8) was used in the form 



to estimate r g for various electronic states Hg, where r e and w for Hg were obtained 
from the results of Wind (refs. 21 and 44). 

All values of D, p, r , and T are tabulated in table n (p. 76). In addition to D, 

Psl ■ 

p, r , and T, equations (17) and (18) also require values of the fundamental frequency of 
vibration o>' (fig. 8). From equation (D4) 


* 

This completes the determination of the constants of Rydberg states of Mg* 


(DIO) 
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APPENDIX E 


SOURCES OF EXPERIMENTAL ENERGY LEVELS OF EXCITED ELECTRONIC 
STATES OF H 2 AND CORRECTIONS TO THEM 

$ 

The energy levels (or term values) of states of have been determined from spec- 
troscopic observations by a large number of investigators , but much of the early work is 
inaccurate. Spindle r (ref, 14) took statistically weighted averages of all experimentally 
determined spectroscopic constants for the ls<?2pcr IscSpo *£, lsa2prr *n, and 
lsaSpTT 1 n states, and his values were used for vibrational term value differences (no 
rotation) and for equilibrium inter nuclear distance r e for these four states (see table I 
(p. 75)). For all states it is necessary to know the absolute value of at least one term 
value in addition to the term value differences. The absolute values are usually mea- 
sured from the rotation less ground vibrational state of An absolute term value for 
each of the four previously mentioned states was taken from the latest measurements 
(refs, 22, 45, and 46) (see table I), For other states all experimental data were taken 
only from the latest measurements (refs, 33, 47, and 48). 

Although for many states Dieke (ref. 33) had the most reliable term value differ- 
ences, some corrections to his absolute term values were necessary because he relied 
on old, far ultraviolet measurements for the absolute term values of excited singlet 
states and found the relative position of excited singlet states and triplet states by extrap- 
olation to the limits of Rydberg series. However, visible, infrared, or near ultraviolet 
line positions were accurately measured by Dieke to determine the term values of all 
excited singlet states relative to each other and of all triplet states relative to each other. 
Consequently, Dieke *s term values are accurate to within an additive constant for excited 
singlets and to within a different additive constant for triplets. 

Field, Somerville, and Dressier (ref. 49) have examined the latest far ultraviolet 
measurements on singlet states (refs_ 22, 45, 48, 50, and 51) and concluded that 700 
meters" * should be added to Dieke's term values for excited singlets, which was done in 
this report (see table I). 

Determination of the correction to Dieke’ s term values for triplet states was more 
complicated because there are no spectral lines connecting singlet and triplet states. 
Kolos (ref. 52) states that Professor R. S. Mulliken believed a correction was necessary, 
and Kolos determined the correction by calculations from first principles. In this 
report reliance was put more on experimental observations, although Wind’s very 
reliable term values and dissociation energy (ref. 21) based on calculations from 
first principles had to be used also. The derivation follows that of Beutler and Junger 
(ref, 53) with two exceptions (I) the energies of all H and states are assumed to be 
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known exactly relative to the rotation less ground vibrationa l state of h| by means of the 
dissociation energy of h| (table in, p. 78) , the term values and ionization potential of H 
(ref. 54), and the term values and dissociation energy of (ref. 21); and (2) the relative 
positions of all triplet levels are known as a result of infrared measurements (refs. 33 
and 55). 

Observations (refs. 56 and 57) of the lines of three electronic transitions were used 
to draw the following conclusions: 

(1) From the transition lscr3pw 3 II — lsa2scr 3 S it was concluded that lsa3p?r 3 n 

(v = 3, K = 3, para) does not predissociate into H(l) + H(2), where the numbers in paren- 
theses following each H are principal quantum numbers. 

O 

(2) From the same transition it was concluded that lscr3pir n (v = 3, K = 4, ortho) 
does predissociate into H(l) + H(2). 

(3) From the transition lscr3 pw 3 n - lso2sa 3 S it was concluded that lscr3{«r 3 H 
(v = 7, K = 2, para) preionizes to H £ (v = 0, K = 0) + e“. 

(4) From the same transition it was concluded that ls<x3pfr 3 n (v = 6, K = 3, ortho) 
does not preionize to Ht (v = 0, K = 1) + e~. 

^ o o 

(5) From the transition lsa3scr S — lsa2prr H it was concluded that lscr3scr 
(v = 3, K = 1, ortho) does not predissociate to H(l) + H(2). 

(6) From the same transition it was concluded that lso3so 3 S (v = 3, K = 2, ortho) 
does predissociate to H(l) + H(2). 

o o 

(7) From the transition Is<x4d6 A — lso2pff n it was concluded that 
lso4d6 3 a (v = 3, K = 4, para) does not preionize to Hg (v = 0, K = 2) + e”. 

(8) From the same transition it was concluded that lsa4do 3 A (v = 3, K = 5, ortho) 
does preionize to Hg (v = 0, K = 3) + e". 

Let d be the correction to be added to Dieke's triplet term values (ref. 33). Then, 
from Dieke’s term values, supplemented where necessary by term differences from 
reference 56 or estimates of term differences from reference 53, and from 
conclusions (1) to (8) it follows that 


11 850 900 + d < 11 837 270 m" 1 

(El) 

11 870 940 + d > 11 837 270 m" 1 

(E2) 

12 514 840 + d > 12 441 300 m~ 1 

(E3) 

12 374 420 + d < 12 447 120 m* 1 

(E4) 

11 851 770 + d < 11 837 270 m" 1 

(E5) 

11 862 270 + d > 11 837 270 m" 1 

(E6) 
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12 481 360 + d < 12 458 720 m 


(E7) 

12 516 970 + d > 12 476 010 m' 1 (E8) 

respectively. From equations (El) to (E8) 

d = -23 82Q±1180 m" 1 (E9) 

Hence. 23 820 meter -1 was subtracted from each of Dieke’s triplet term values (see 
table I, p. 75). This puts lscr2pji n within 760 meter of the term values relative to 
Hg calculated by Beutler and Jiinger (ref. 53) and within 1460 meter -1 of the term values 
relative to Hg calculated by Richardson (ref. 56). Beutler and Jiinger and also Richardson 
had available less complete, less accurate data and calculations from first principles 
than were available for the present study. 
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TABLE X. - REFERENCES FOR EXPERIMENTAL TERM VALUES OF EXCITED 


ELECTRONIC STATES Hg AND CORRECTIONS TO THEM 


Excited electronic 

Reference for 

Reference for 

Reference for 

Amount added 


state 

vibrational 

vibration - 
rotation 

absolute 

term 

to absolute 

Term 

Small- 

term value 

term values 

symbol 

internuclear~ 

distance 

designation 

differences 
(no rotation) 

term value 
differences 
or r e 

values 

obtained 
from reference 
in preceding 

column,^ 

1 

m 

l Z 

Iso2scr 

33 

33 

33 

700 


I so 3 so 

33 

33 

33 

700 

'a 

lso2po 

14 

14 

22 

0 


lso3po 

14 

14 

45 

0 


lso4po 

47 

47 

47 

0 

X 0 

lso2p^ 

C 14 

14 

46 

0 


lso3p? 

C 14 

14 

45 

0 


lso4pir 

C 4S 

47 

48 

0 


Isolde 

33 

33 

33 

700 


lsa4do 

<d> 

33 

33 

700 

l u 

IsoSdf 

(C) 

33 

33 

700 


lso-4dflr 

(C) 

33 

33 

700 

l A 

IsoSdO 

(c) 

33 

33 

700 


!so4d8 

(c) 

33 

33 

700 

h 

lso2so 

33 

33 

33 

-23 820 


lsa3so 

33 

33 

33 

-23 820 

h 

IsoSpo 

33 

33 

33 

-23 820 


lsa4po 

33 

33 

33 

-23 820 

h 

iscr2pjr 

(c) 

33 

33 

-23 820 


lso3pf 

(e) 

33 

33 

■23 820 


lso4p^r 

(c) 

33 

33 

23 820 

3 S 

Isolde 

33 

33 

33 

“23 820 


lso4do 

33 

33 

33 

-23 820 

3 n 

lso3dff 

fc) 

33 

33 

-23 820 


lso4d^ 

fc) 

33 

33 

-23 820 

3 a 

iso3d6 

(c) 

33 

33 

-23 820 


Isa4d6 

fc) 

33 

33 

-23 820 


rt 1 

^Measured from rotationless ground vibrational state (leolso Z) of 
fe See appendix E lor explanation. 


c Virtual levels. 
^Not given. 



TABLE E, - CONSTANTS AND SUBSCRIPTS OF EXCITED ELECTRONIC STATES B* 


Excited electronic 
state 

Elec- 

tronic 

Principal 
quantum 
number 
n of 
excited 
dissoc- 
iated 
atom a 

Bump 

height* 

*1 

m 

Well 

depth, 

Morse 

parameter, 

ft 

m- 1 

Equilibrium 
interim clear 
distance, 

V 

m 

Energy 
relative to 
electronic 
ground 
state, 

T. 

m -l 

m 

Subscripts 

tn 

computer 

program 

Term 

symbol 

Small- 

internuclear^ 

distance 

designation 

degen- 

eracy* 

&e 

h 

leaSstr 

1 

2 

b o 

2.0391x10 s 

l,4529x!0 10 

1*0334x10' 10 

9, 798l5xlO S 

2, 1 


lso3se 

1 

3 

b 0 

2. 1821 

1.3907 

1, 1181 

11. 17846 

3* 1 


lso4sa 

1 

3 

C 8G 170 

1.6923 

d 1.3448 

*1, 1290 

f 11. 74843 

M 

»r 

lso2pu 

1 

2 

*0 

2 r 8848xlO S 

0, 6948x1 0 i0 

1.2907x1 0 -10 

8. 95240x10 s 

2,2 


lso3po 

1 

2 

b o 

. 8796 

1.7070 

1, 1270 

10. 95767 

3,2 


lscr4pe 

1 

3 

C 0O no 

1,8591 

1.3750 

1. Ill 

11. 58169 

4, 2 

l u 

lso2pjj 

2 

2 

h 10 550 

2, 0653x1 0 S ! 

1.5493xl0 1S 

1. 0316x10" 10 

9. 78257X10 6 

2, 3 


lscr3pir 

2 

3 

b 0 

2. 1907 

1.3895 

1, 0440 

11. 16907 

3, 3 


lecr4pir 

2 

4 

l G 

2.2249 

1.3451 

1. 0555 

11.66084 

4,3 

*£ 

lSo3do 

1 

2 

b 241 970 

3. 0272x1 0 6 

2. 2770X1O 10 

h. 0740xl0' 10 

11. 05203X10 6 

2,4 


lse4do 

1 

3 

c 80 no 

1.8310 

d l. 5551 

1. 4205 

11.60980 

3,4 

l a 

lso3dv 

2 

2 

b 241 970 

1. 0053x10 s 

2,265axl0 i0 

1. 1639X10^ 10 

11,07394x10® 

2, 5 


latr4dtr 

2 

3 

C 0G 170 

1. 7965 

1.4662 

1.3972 

11. 64430 

3,5 

>A 

IspSdfi 

2 

3 

b 0 

2.2368xl0 S 

1.3204X10 10 

0. 99 16x10' 10 

11. 12376x10 s 

M 


Isolds 

2 

4 

k o 

2.2431 

d 1.3281 

.8959 

11,65069 

3, 6 

l r 

lso4fc 

1 

3 

c 90 170 1 

1. 8416x10 s 

*1, 4068X10 10 

e -10 
1. 0809x10 u 

m 11. 59920xlO S 

2, 7 

*11 

lso4fjt 

2 

3 

‘73 600 

H. 76a6xlO S 

i 1.6348xl0 10 

*1. 0583x10' 10 

*11. 66560xlO S j 

2 t fl 

J A 

Isc4f5 

2 

3 

c 80 170 

1. 8327x10 s 

i l. 4068xl0 10 

e l, 0821x10" 10 

m ll. 60810x10 s 

2, 9 

1 

<p 

lso4f ^ 

2 

4 

k o 

2.2751xlO S 

1 1.4068xl0 10 

e 1.0253x10' 10 

m il. 61865x10 s 

2, 10 

h 

IbpZsu 

3 

2 

b 0 

2. 4716x10 s 

1. 4844X10 10 

0. 9955x10" 10 

9. 36566x10 s 

2, 11 


lso3sa 

3 

3 

b o 

2. 3004 

1.3641 

1* 0613 

11. 06023 

3*11 


lscr4so 

3 

3 

c so no 

1. 8088 

d 1.3385 

e l* 1132 

^ll. 63219 

4, 11 

h 

lsd3pa 

3 

2 

b G 

l,2957xlO S 

1.8309x10*° 

1. 0869x10 ' 10 

10. 54161xlO S 

2. 12 


lSff4po 

3 

2 

c 309 210 

.7202 

2.6694 

1.0654 

11,42631 

3, 12 

3 n 

lso2p?j 

6 

2 

b 0 

2.4785x10 s 

1.3497*10 10 

1. 0499xl0‘ 10 

9. 35873x10 s 

2, 13 


IscSpr 

6 

3 

b 0 

2,3161 

1.3407 

1.0667 

11. 04446 

3, 13 


IscHpr 

6 

3 

c 80 no 

1.8545 

1. 6084 

1. 0603 

11. 58621 

4 t 13 


^Excited electronic states of K^, this dissociate to one ground state atom and one excited H atom. 
b ReL 41* 

c CaIculated assuming H/6 - 0. 15036. 

^Interpolated between n = 3 and along Rydberg series, 

e Assumed value of same as for Kg. 

f Ref. 58. 

e Ref. 59. 

h Ref. 15. 

Energy level s t hump height t and from ref, 42. 

^Irregular term values; thus, r fi determined by method on p. 112 of ref. 56. 

Assumed value of 0 because of unprom oted MO. 

^Interpolated between n = 3 (average for all states) and H^ + 
m Estimated. 
n Bef . 60. 
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TABLE U. - Concluded. CONSTANTS AND SUBSCRIPTS OF EXCITED ELECTRONIC STATES Eg 


Ex cited electronic 

Elec- 

Principal 

Hump 

Welt 

Morse 

Equilibrium 

Energy 

Subscripts 


state 

tronic 

quantum 

height, 

depth, 

parameter, 

inter nuclear 

relative to 

in 

Term 

Small- 

degen- 

number n 


D, 

a, 

distance. 

electronic 

computer i 

symbol 

intern u clear- 

eracy t 

of excited 

m 

m“* 

m " 

r e> 

ground 

program 


distance 


dissoc- 




m 

state. 



designation 


iated 





T, 





atom a 





m' 1 


3 S 

lso3do 

3 

2 

b 330 700 

1. 1139x10® 

1.9680X10 10 

1.446 1X10“ 10 

11.05411x10® 

2, 14 


lso4d(7 

3 

3 

c 80 170 

1, 8156 

!, 4693 

j.9298 

11.62514 

3,14 

3 n 

Iscr3d?r 

6 

2 

b 330 700 

1. 1000x10® 

2. 1220x1 0 10 

1. 1466x10* 10 

11. 06801x10® 

2, 15 


lso4d3T 

6 

3 

c 80 170 

1.8129 

1.5022 

1. 1676 

11.62789 

3, 15 

3 a 

Isa3d6 

6 

3 

b 0 

2.2430x10® 

1. 3103x10*° 

0. 9937x10' 10 

11. 11759x10® 

2, 16 


lso4d5 

6 

4 

k o 

2.2575 

1.3086 

.8348 

11.63629 

3,16 

3 r 

lso4fo 

3 

3 

C 8G 170 

1. 8037x10® 

*1. 4068X10 10 

e L O366X10 -10 

n ll. 63704x10® 

2, 17 

3 n 

lso4f?r 

6 

4 

k o 

2.2551x10® 

1 1. 4068X10 10 

e l. 0276x10" 10 

n ll, 63865<1G® 

2, 18 

3 a 

Isc4f5 

6 

3 

C 80 170 

1.8327x10® 

'i.4068xio 10 

e l. 0821* 10" lb 

m 11. 60804x10® 

2, 19 

V 

lsa4£<p 

6 

4 

k o 

2. 2752 <10® 

1 1. 4068xl0 10 

e l. 0253x10" 10 

m ll. 61859x10® 

2,20 

Various 

lso5 

10G 

3 to 5 

Various 

2. 0408x10® 

1 1. 3556X10 10 

e l. 0732 -10" 10 

11,87246x10® 

2,21 


1s<t6 

144 

4 to 6 


2 r 1083 

' 1. 3409 

e l. 0703 

12, 01307 

3, 21 


1sct7 

196 

4 to 7 


2. 1432 

^ 1. 3352 

e l. 0676 

12,09721 

4,21 


lsa8 

256 

5 to 8 


2. 1710 

1 1,3330 

e l. 0657 

12. 15153 

5, 21 


1so9 

324 

6 to 9 


2,1373 

*1.3319 

e l. 0641 

12. 18862 

6,21 

Various 

IsolO 

400 

6 to 10 

Various 

2,2005x10® 

1 1 3314<10 10 

e l. 0628' 10" 10 

12,21506M0 6 

2, 22 


Iso 11 

484 

7 to U 


2.2091 

1 1,3311 

e l. 0617 

12.23457 

3, 22 


lsol2 

576 

8 to 12 


2.2159 

' 1.3309 

e l. 0610 

12,24938 

4, 22 


lsal3 

676 

9 to 13 


2.2214 

1 1.3308 

e l. 0605 

12, 26083 

5. 22 


IsoH 

784 

9 to 14 


2. 2253 

* 1.3307 

e 1.0600 

12.27000 

6,22 

Various 

lsol5 

900 

10 to 15 

Various 

2.2291-10 6 

1 1, 3307xl0 10 

e l. 0595x10" 10 

12. 27734x10® 

2, 23 


Iso 16 

1024 

11 to 16 


2.2319 

1 1, 3307 

e l. 0592 

12.28335 

3.23 


lso!7 

1156 i 

11 to 17 


2.2343 

1 1. 3307 

*1. 0586 

12,28832 

4, 23 


b Ref . 41. 

Calculated assuming H/5 = 0, 15036, 


e Assumed value of wr 2 same as for Hj. 

e l 

J Ir regular term values t thus, determined by method on p. 112 of ref. 56. 
k 

Assumed value of 0 because of unpromoted MO. 

1 Interpolated between n = 3 (average for all states) and K^. 
m Estimated. 
n Ref, 60. 
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TABLE HL - DISSOCIATION, PHOTODETACHMENT, 


AND IONIZATION ENERGIES 


Reaction a 

Energy absorbed 

L 

in reaction, 
m- 1 

Reference 

h| (lSCTlsff *2) - 2H* 

w| = 3. 61136X10 6 

46,61 

Hg (Iso 2 S) - H* + H + 

Wg = 2. I3794xl0 6 

21 

H 3 < 1 a i> - 2H * + H+ 

W 7 = 7. 2707x10® 

13 

H' ~ H* + e' 

l g = 6. 2187x10® 

62 

— H + + e' 

l| = 1. 096788X10 7 

54 


a Molecules and molecular ions are assumed to be in 
their rotationless ground vibrational states. 

^All values given are differences in term values 
(energy /he) of reactants and products; hence, 
units are m” 1 
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TABLE V. - NUMBER DENSITIES OF SPECIES IN HYDROGEN PLASMA a,b 
(a) Pressure, 1. 0l325xl0 5 N/ro 2 (1 atm) 


Temper^ 

ature. 


H 2 

++ 

X 


* 

H 


e“ 




| 


| H 2 


H + 


















T, 

K 







c Number density, N. , ro ~ 







300 

2,447£ 

25 

0* 


1 * 03 OE 

-10 

0. 


0- 


0, 


0, 


0- 


0- 

400 

1- 035E 

25 

0- 


2, 449E 

-01 

0* 


0* 


3- 


3 * 


0, 


0- 

500 

1.468E 

25 

0- 


1 -021E 

05 

0, 


0* 


3, 


3* 


0- 


0- 

600 

1.2 2JE 

25 

0, 


5* 64 f£ 

08 

0* 


0- 


3- 


0* 


0, 


0- 

700 

L*049E 

25 

0. 


2 , 646 E 

11 

0- 


0. 


0, 
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a See fig. 2(a). 

b With coulomb interactions between tree charged particles for T > 2000 K r 
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TABLE V, - Continued. NUMBER DENSITIES OF SPECIES IN HYDROGEN PLASMA* 3 
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23 

15000 

3. 34 3 E 

20 

8.744E 

la 

3. 50UE 

24 

3.U09E 

22 

6.994E 

23 

2 -5U6E 

20 

' 4,2 96b 

16 

2* 30 7E 

20 

6 * 994E 

23 

16000 

U736E 

20 

7.7S9E 

is 

2. 775E 

24 

3.955E 

22 

9. 106t 

23 ; 

2-2636 

2 0 

2E 

16 

2.080E 

20 

9. 10 6 E 

23 

isooo ; 

3.S34E 

1 9 

4.495E 

IS 

1. 534t 

24 

5. 355E 

21 

L.279E 

24 ; 

I * 3S4E 

20 

7 . 148 E 

| 

15 

L.272E 

20 

L, 279= 

24 

20000 

6.731E 

la 

l.ai 2 fc 

16 

I 7.352b 

23 

6.664E 

22 

1 . 475E 

24 

5 .2 14E 

19 

l . 293b 

15 

5. 72 BE 

19 , 

L.4 75E 

24 

22000 

1.090b 

IS 

6 . 1B5E 

17 

' 3* 31 4t 

23 

5. 1 7fi£ 

22 

1*5096 

24 | 

2 * 3S4E 

19 

1.925E 

14 

2, 2076 

19 

1.509E 

24 

24000 

L,082E 

1 7 

2 , 06 Ob 

17 

1. 52Gfc 

23 

4.4S6E 

22 

1 .4 58E 

24 

8 + 9b2 1 

18 

2.B75E 

13 

B, 32BE 

IS 

1 . 4 5a E 

24 

26000 

3. t 22 b 

16 

. 7.2 S4E 

16 

7. 372 E 

22 1 

3.7S3E 

22 

1 .378E 

24 1 

3 .54 V t 

ia 

4.745b 

12 

3 * 305E 

IB ! 

1. 37 BE 

24 

2 BOO 0 

8.54LE 

15 

2 . Slot 

16 

3,616b 

22 ; 

3.241b 

22 

1 *2 93E 

24 | 

1,5326 

LBj 

a *B46E 

1 l 

1 *406E 

LS 

1 *2936 

24 

.100 00 

2,254b 

1 5 

1 . 17 7E 

16 

2.103E 

22 

2*7S2E 

12 

1.214b 

24 

6.857E 

L 7 

1 .9 19fc 

11 

6.452b 

17 

1,2 L 4b 

24 

12000 

6 . 7 2 BE 

14 

5.269E 

15 

1.225E 

22 

2.4 L4E 

22 

1 - 141E 

24 

3.345b 

£7, 

4 . 7 6 B E 

10 

3.16GE 

17 

L* 1416 

24 

34000 

2.236E 

14 

2 * 547E 

1 5 

7. 4 sat 

21 

2. 1 02E 

22 

1 1.075E 

24 

1 . 7 3 L a 

17 

L. 33 4 E 

10 

1 .641 £ 

17 

1.0 75E 

24 

36000 

a. 1 4 h E 

13 

1.296E 

*5 

4. 773 t 

2 1 

1 . S63E 

22 

1.016E 

24 

9.433E 

16 

4. 136E 

09 

S.973E 

16 

1.0166 

24 

3600Q 

3.2176 

13 

6 .96 at 

14 

3- 155E 

21 

L . 6 5 7 E 

22 

9. 631E 

23 

5 -379E 

1 6 

1. 403E 

09 

5,1 33E 

16 

9.6 3 LE 

23 

40000 

1. 162E 

13 

3.0776 

14 

2. 153b 

21 

1.483E 

22 ' 

9.1526 

23 

3 * 191 E 

1 6 

j 5.145b 

OS 

3.054E 

16 

9. 1526 

23 


^With coulomb interactions between free charged particles for T > 2000 K. 
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TABLE V, - Continued. NUMBER DENSITIES OF SPECIES IN HYDROGEN PLASMA 5 


(c) Pressure, 1, 01325*H> 7 N/m 2 (1O0 atm) 


Temper ~ 

4 


H 2 


H* 


H* 


e~ 


IT 


4 


h 2 


H + 


ature 





| 














T, 

K 







Number density. 

K 

m 








400 

1 ,8 356 

27 

0, 


2 „ 44 9 E 

00 

0, 


0- 


0, 


0. 


0, 


0. 


600 

1 * 46 at 

27 

0, 


1-02 It 

Oo 

0, 


Op 


0, 


! 0 . 


0, 


0, 


6 00 

1*2 2 if 

27 

0- 


5,64 ft 

09 

0, 


0* 


0* 


Op 


0, 


0- 


700 

1 * 0 49 E 

27 

0* 


2, 646b 

12 

0, 


0. 


0, 


0- 


0- 


0- 


800 

9- 1 75E 

26 

0, 


2*6526 

14 

0* 


0. 


0 - 


0 . 


0. 


0* 


900 

8* 15 3fc 

26 

£>, 


y, 5006 

1 5 

0, 


Op 


0, 


0. 


0. 


0. 


1000 

7* 340c 

26 

0- 


t * 66 1 1 

1 1 

0- 


0- 


0- 


0- 


0. 


0. 


UOQ 

6,6 72 E 

26 

0. 


1 , 72 Ot 

18 

0, 


0. 


0 . 


0 . 


0. 


0 * 


1200 

6*1 1 q if 

26 

0. 


1.2038 

1 9 

0 , 


0. 


0 - 


0. 


0- 


0- 


i 300 

5,64 6E 

26 

! 0, 


6-223E 

19 

0, 


0, 


0 - 


0. 


0- 


0- 


1400 

5.243E 

26 

0, 


2,5ME 

20 

0- 


2, t 73E 

05 

5*5356 

0 1 

2,1716 

05 

0* 


0. 


1500 

4,6431 

26 

0. 


B-5B3t 

20 

0. 


4,8 55E 

06 

2.B956 

03 

4,B58E 

Ofa 

0. 


0- 


L6 00 

4, 587 £ 

26 

0, 


2 * 4 85 l 

21 

o. 


7.3696 

0 r 

7*9556 

04 

7, 3776 

07 

Op 


0- 


1000 

4 - U 7 7 1 

26 

9 . L 7ab 

-05 

1*45^6 

22 

1 * 6 l SE 

-Q6 

6 * B 4 6 E 

0^ 

1,3536 

0 T 

6-8866 

09 

0. 


0- 


200 0 

3.6696 

26 

1 * 359E 

- 01 

5, 947fc 

22 

4 ■ 7 4 46 

-03 

2 .5886 

1 1 

1,5636 

09 

2.6036 

U 

0 . 


0. 


22JO 

3-3346 

26 

4,95ut 

01 

l - 8 72 1 

23 

3-2406 

00 

5 . 04 7t 

12 

3*5376 

LO 

5- 102E 

12 

7.979E 

05 

6.5L lb 

05 

Z4uG 

3-U53E 

26 

6 , 7 7 tit 

03 

4, 84 3 1 

23 

7 , 42 Ot 

02 

6 .005 E 

13 

1 -366E 

12 

6- 1126 

13 

6,3596 

07 

6- 3636 

07 

2600 

2,8126 

26 

4, 3tB6 

05 

1- 07 7t 

24 

7, 332E 

04 

4 . 8026 

14 

1.2846 

13 

5,0106 

14 

2,5906 

09 

3 . 0 8 BE 

09 

zeou 

2,600E 

26 

1 -5546 

07 

2, 12 8t 

24 

3, 742E 

35 

2.9406 

15 

1.0686 

14 

3,0466 

15 

6, 2 73E 

10 

0.6 406 

10 

100 0 

2-4086 

26 

3 . 430E 

OB 

3, 0 1 7 fc 

24 

1 * 125E 

OB 

1,3 91 6 

L 6 

6,6076 

14 

L.457E 

L6 

9.9396 

11 

1.55TE 

12 

32 ou 

2 ,2 30E 

26 

5- 12 76 

09 

6, 324b 

24 

2* i 97t 

09 

5.402E 

1 9 

3-2086 

L 5 

5 ,7 L3E 

16 

I . LITE 

13 

1 + 960b 

13 

3400 

2 ,U60E 

26 

5-54'it 

1 0 

9, 821 1 

24 

3 * 00 BE 

LO 

L . 7816 

t 7 

L ,2 706 

18 

l *9056 

17 

9,4546 

13 

1-8396 

14 

3600 

l,d95£ 

26 

4 ,56 it 

U 

1.440t 

25 

3 , 05 7 1 

L L 

5- 1 I4fc 

17 

4,2426 

16 

5,5186 

17 

6. 305E 

14 

1.3496 

15 

36 00 

1.7306 

26 

2* 9736 

12 

2,01 IE 

25 

2.4156 

12 

1,306 E 

i d 

1.223b 

17 

1 . 4 1 fa E 

L 8 

3-4356 

15 

8.0456 

15 

40 00 

L - 56 6t 

26 

1,5846 

13 

2 , 690 1 

2 fa 

1 ,53 78 

L 3 

3,01 IE 

18 

3- 1066 

1 7 

3-2666 

13 

1.5716 

16 

4-OZOt 

16 

4300 

1 .3206 

26 

l - 40 9E 

14 

3- B71 1 

25 

1 * 7 546 

1 4 

0,9976 

1 3 

1,0256 

18 

9. 5686 

IS 

| 1 * 1 6 L E 

17 

3.389E 

17 

4600 

1 . 07 4 E 

lo 

0.9906 

14 

5, 1 6 7 1 

25 

1 -42 IE 

15 

2 * 3Q6E 

L9 

2.7686 

18 

2 . 3056 

19 

6, 3 90 1 

17 

2- 1 42E 

18 

4900 

0.5256 

25 

4-345t 

15 

' 6 4 45 3 b 

25 

0. 66 1 E 

15 

5 *2866 

1 7 

6- 397 1 

18 

4.6056 

19 

2,7056 

10 

t , 0 5 OE 

19 

5200 

6.5046 

25 

1 - 65 4t 

16 

7- 60tit 

25 

4, 1638 

16 

1.1286 

20 

1 -3256 

19 

7.6396 

L9 

8, 9526 

ia 

4. 075E 

19 

5600 

4- 330E 

25 

7*2186 

l 6 

8- 7 76 1 

25 

2 . 4 9 OE 

17 

2-9206 

20 , 

3 , 1 3 1 E 

1 9 

1 - 141E 

20 

3. 1 L2E 

19 

1.7026 

20 

6000 

2, 7 72 £ 

25 

Z * 3 7 5t 

l 7 

9-461E 

25 

L - 12 56 

IB 

7, 1 50E 

20 ! 

6 , 7 0 9 fc 

13 

1 - 3286 

20 

7.B28E 

19 

5, 72 OE 

20 

6500 

1.5546 

25 

7 * 828t 

1 7 

9* 73 71 

25 

5 , 4 7 OE 

18 

1.9886 

21 ; 

1.5166 

20 

1,3066 

20 

1 ,B30E 

20 

L-828E 

21 

7000 

B-8L2E 

24 

2 , 04 2t 

IB 

9.6036 

2 5 

2 . 04 4 b 

l 9 

4.8726 

21 

2-9716 

20 

1 - 1596 

20 

3.4926 

20 

4. 705E 

21 

7500 

5- I69t 

24 

4*5zt>6 

IB 

9,26 7 b 

2 5 

6,439b 

19 

1 - 0 5 6 E 

2 2 

5. L93E 

20 

9 . 9fa9E 

L9 

5-8376 

20 

L.Q49E 

2 2 

a oou 

3-1926 

24 

B-9L3E 

1 8 

B * 854E 

25 

L, 735E 

20 

2,1176 

22 

8 ,306E 

20 

B. 573E 

19 

B . 9786 

20 

2- L02E 

22 

90 00 

1,34 It 

24 

2.692d 

19 

8- 00 7t 

25 

9, 072E 

20 

6 ,5746 

22 

1.7526 

2 1 

6-5856 

19 

1- 7046 

21 

6- 66 4E 

22 

1000 0 

6- z>3Bt 

23 

6 +4 1 Bt 

19 

7-2406 

25 

3 - 4 3 9t 

21 

1 -679E 

23 

3,0836 

21 

5,3576 

19 

3-0176 

21 

1.680E 

23 

1 luuo 

3,5426 

23 

1 . 29 Ot 

20 

6, 5656 

25 

L, 0326 

22 

3.5856 

23 

4.76BE 

21 

4*5496 

19 

4,5506 

21 

3-507E 

23 

12000 

2,0 69 6 

23 

2- 27 5t 

20 

5. 962E 

25 

2,5866 

22 

6,7506 

23 

5*6856 

2 1 

3,9566 

19 

t* 279E 

21 

6- 754E 

23 

l 30 GO 

1.2/26 

23 

3.602E 

20 

5, 405t 

25 

5.635t 

22 

1.1516 

24 

8.6566 

2 L 

3 .46 26 

19 

8.0476 

21 

1, 1 52E 

24 

14000 

6 - J 66c 

22 

5 * 1B0E 

20 

4* 07&t 

25 

1. 0946 

23 

1 .8 L06 

24 

1 *0466 

22 

3 .0006 

19 

9 , 66 7 E 

21 

1 . a L 16 

24 

15000 

5.143E 

22 

6 , 05at 

20 

4, 361 t 

25 

L- 924t 

23 

2 , 659E 

24 

1 - L08E 

22 

2 ,5376 

19 

1-0946 

22 

2.660E 

24 

160 00 

3 * 3 5 OE 

22 

7. 86 76 

20 

3, B56t 

25 

2-B9fa£ 

23 

3 , 6 B 4 £ 

24 

1.2726 

22 

2,0726 

19 

1-1696 

22 

3-685E 

24 

ISO 0 0 

1, 350E 

22 

7- 30 7t 

20 

Z, B80t 

25 | 

4.7196 

23 

6,09 L6 

24 

1 -2376 

22 

l. 199E 

19 

l- 137E 

22 

6.092E 

24 

20UO U 

4. 472 t 

21 

5 * 3B3E 

20 

i, 99BE 

25 

6-4606 

23 

8.5266 

24 

9,7616 

21 

5,520t 

IB 

8,9996 

2\ 

8- 52 76 

24 

22000 

1*647 £ 

21 

3-2426 

20 

1-2846 

25 

7 - 40 IE 

23 

1.045E 

25 

3 * 4 2 2 £ 

2 L 

2 ,0 1 7£ 

18 

5,9436 

21 

1.0456 

25 

24000 

5,0 21 E 

20 

1 +6 7 4t 

20 

7, 850t 

24 

7.5156 

2 3 

l , L59E 

25 

3,8816 

2 1 

5 , 101 E 

17 

3-4216 

21 

1. 1596 

25 

26000 

L,476fc 

20 

7.916b 

19 

4, 642 1 

24 

7, 0&7t 

23 

l . 200E 

25 

1*9436 

2 i 

1.6396 

17 

I . 8136 

21 

1.200b 

25! 

2BQG0 

4, 195t 

19 

3 ■ fa 2 36 

19 

2* 73dt 

24 

fa* 346E 

23 

1 - 1 92E 

23 

9,9346 

20 

4* 196£ 

It 

9.3106 

20 

1 - 1 9 Zb 

25 

30000 

1 *3736 

19 

1 - 666fc 

19 

1 , 642 h 

24 

5, 554E 

23 

1,1576 

25 

5.10 16 

20 

1 , 1 L5E 

16 

4.800t 

20 

l* 157E 

25 

32000 

4-5856 

LB 

B.022E 

IB 

UOUt 

24 

4. 79dt 

23 

l . 1096 

25 

2,685£ 

20 

3- 1596 

15 

2-536E 

20 

1 . 1 Q9E 

25 

34000 

1 . 6 4 1 £ 

1ft 

4* 0326 

IB 

6 , 4 1 6 E 

2i 

4, 24 1 1 

23 

1.057E 

25 

1,4596 

20 

9 * fa 30 E 

14 

1*3836 

20 

1.05 76 

25 

36000 

6* 304t 

17 

2,04 7b 

IB 

4, L99t 

23 

3, 754t 

23 

1,0066 

25 

3*2 17E 

19 

3* 1696 

14 

7-8 1t£ 

19 

1- OObb 

25 

36000 

2, 58 Be 

l 7 

1 - 1426 

IB 

2* 830t 

25 

3.3136 

23 

9.57B6 

24 

4,7986 

19 

1 , I22£ 

14 

4.578E 

19 

9 . 5 7 0E 

24 

4UGU0 

1.1286 

1 7 

6*479t 

17 

1 , 959 t 

23 

2,960E 

23 

9.125E 

24 

2 * B96 E 

19 

4 ,2486 

13 

2-7716 

19 1 

9- l 25E 

24 


b With coulomb interactions between free charged particles for T > 2000 K. 
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TABLE V ♦ - Concluded. NUMBER DENSITIES OF SPECIES IN HYDROGEN PLASMA b * d 


(d) Pressure, 1. 01325x3 0^ N/rn 2 (1000 atm) 


Temper- 



* 

H 2 


H* 


* 

H 


e“ 


H" 


H 3 


H 2 


H + 


ature, 















T t 

K 







Numbe 

:r density, N. 

, m- 3 








1600 

4-60/f 

27 

0* 


7, 859E 

21 

0, 


4,1 J9E 

08 

1 *413c 

06 

4. L53£ 

08 

o. 


' 0* 


iaoo 

4# 0786 

27 

9,1766- 

04 

4.6016 

22 

5. 1 LQE 

“06 

3. 8496 

10 

3 -4b 36 

08 

3,8846 

10 

0. 


0, 


2000 

3*67Q£ 

2 7 

U36UE 

00 

u aatt 

23 

1 * 5Q0£ 

-02 

1*4466 

12 

2* 7626 

10 

1 ,4746 

12 

0. 


0, 


220 0 

1.136c 

2 7 

4,9526 

02 

5. 92 ot 

23 

1,0256 

01 

2.B06E 

13 

J. 733£ 

u 

2.9036 

13 

1*4356 

06 

3* 704E 

05 

24 00 

3*057£ 

2 7 

6- 7066 

04 

1 , 5326 

24 

2*34 76 

03 | 

3.3 lae 

14 

i,8b36 

13 

3.5046 

14 

1* 152E 

oa 

3* 64 4E 

07 

2600 

2,3l4£ 

2 7 

4,3796 

06 

3*4 1 it 

24 

2. 321 £ 

05 

2*6776 

15 

2.2296 

14 

2,9006 

L5 

4,7496 

09 

U 7836 

09 

2auo 

2-b L5E 

27 

1* 5626 

08 

6,74/fc 

24 

1, id6£ 

07 

1,6006 

15 

1*3446 

15 

L -7S5E 

16 

1,1596 

11 

S.034E 

10 

IGOG 

2*4 34E 

27 

3, 4666 

09 

1 , 214t 

25 

3.5T2E 

OB 

7 ,5226 

16 

1 .1366 

15 

8*6586 

16 

1*8586 

12 

9.1 53£ 

11 

3200 

2,27 IE 

27 

5*2246 

LO 

2,0216 

25 

7.004c 

09 

2,9G9E 

1 7 

5.510E 

Lb 

3*4596 

17 

2-1166 

13 

l. 1636 

1 3 

340 0 

2, 12 76 

27 

5,71 at 

11 

3, 15&E 

25 

9. 641t 

10 

9,5756 

17 

2 ,1946 

1 7 

L , L 77 6 

ia 

I .8 t7E 

14 

1- LG 06 

14 

1600 

U992E 

27 

4 , 79 Ot 

12 

4 , 6 7 UC 

25 

9*0 79t 

11 

2,7576 

13 

7,4 146 

17 

3.4966 

Lb 

1 .2 32E 

15 

8- L29E 

14 

3Bl>u 

1 * ubb£ 

27 

1, 149t 

13 

6, o04t 

25 

7 , ti 9 1 £ 

12 

7,oa?t 

H 

2.13 16 

18 

9*2566 

IS 

6. 8396 

15 

4.878c 

15 

4000 

1- 746E 

27 

L * 7 ft it 

14 

6- 97dt 

25 

5 » 1 G 06 

1 3 

1 ,6556 

19 

5,6936 

18 

2,2 19E 

19 

3,1996 

16 

2.4516 

16 

4100 

l ,b 736 

27 

1 ,6 72c 

15 

1,33 7c 

2b 

b, 0016 

14 

5 , 075 1 

19 

1.997b 

l 9 

7, 0266 

19 

2,4686 

17 

2.088E 

1 7 

46 00 

i * 4 oot 

27 

L , 1 atib 

lb 

L, 6676 

26 

5,0676 

15 

1 . 3 39E 

20 

5,3076 

19 

1 ,d92£ 

20 

1 , 45 IE 

L8 

1*3466 

18 

4900 

1 . 26 Lc 

2 t 

6,3146 

lo 

2,4726 

26 

3.258E 

16 

3. I25E 

20 

1.449c 

20 

4,4376 

20 

6 * 8 05 E 

18 

6 . 8 9 8t 

ia 

62 00 

L.G99E 

27 

2,755c 

17 

3, 125E 

2b 

i . 670E 

1 7 

b . 5 9 96 

20 

3,1 846 

20 

9 , 2 28 E 

20 

2*6326 

19 

2,9166 

19 

6600 

9.0866 

26 

1 .4846 

ia 

4,020fc 

2o 

1 * 10 Jfc 

ta 

1,5806 

21 

7,7606 

20 

2, 078 t 

21 

1.2 38 6 

20 

1 , 5 4 7 £ 

20 

6000 

7 , i5ae 

26 

6- 125E 

ia 

4, 874E 

2b 

5,5 3 7c 

i a 

3.4306 

21 ; 

1,54 it 

21 

3.949E 

21 

4,5176 

20 

6, 406i 

20 

6600 

6*4 99E 

26 

2, 654b 

1 9 

5,7926 

2b 

3,05 46 

t 9 

7. 9896 

21 

3.624 t 

2 l 

7- 122£ 

21 

Ub7 it 

21 

2- 81 36 

21 

7000 

4,00d£ 

2b 

6* 7*56 

19 

6* 476 1 

26 

1 . 2 d Oc 

20 

1 ,7366 

22 1 

7,1426 

21 

1,0496 
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b VVith coulomb interactions between free charged particles for T > 2000 K. 


d See fig. 2(b)- 
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TABLE VI. - INPUT FORM WITH DATA FOR SAMPLE 


PROBLEMS FOR COMPUTER PROGRAM FOR 
COMPONENTS OF HYDROGEN PLASMA 


f Top row of numbers are card column numbers. ] 


1 10 

11 25 

26 30 

DH 

PD 

1 . 01325E+8 

IPR 1 

- oJ 

1 . 


1 10 

11 15 

T 

I FLAG 

30000. 

1 


1 10 

11 25 

26 30 1 

DH 

PD 

IPR ) 
1] 

1 . 

1 . 01 325E+8 - 


1 10 

11 15 

T 

I FLAG 

500. .. ...... 

0 

1500. 

0 

10000. 

0 

20000 r 

0 

30000. 

a. 

40000. 

a 


a Case I, 
^Cases 2 to 7. 
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TABLE TO, - OUTPUT FOR SAMPLE PROBLEMS FOR COMPUTER PROGRAM FOR 


COMPONENTS OF HYDROGEN PLASMA 


a {ON |»0OOOE tit) P(j 1,01326 OS J. PH 0 


PRESSURE I*0l32£ OS N/M**2 


3*00006 04 K 


oensity 2,8727e-3L lg/m**3 


TOTAL NUMBER DENSITY 2*64386 lb 17M**3 


CQJLOH& COMPRESS! HI ul 17 0,92540 

R6U PHOCAL 

DEBYE LEMS IN 

I * 1 4 1 5 £ ua i/M 


63JIV, C39C*/CFU TICAL E^DIV, 

CONS* 4*6571£-0I 

SPECIES NUMBER DENSITY* P AR T 1 C L £ S /C U6 1 G 

ME Ttft 










H ?*8l65E 25 H* 9. 2¥24fc 25 

6 

9*29146 

25 H2 

3.5S7&E 22 

rl- 

1*75256 23 

H2* 1,64946 

2 3 

Nit L*117e£ 20 

DETAILED NUMBER DENSITY FUR H2, HOLE t ULfc S/C UB L C 

METER 









GfttlUNJ LL 6CTRUN IL STATE Hi ?H VlPRAl 

I UNA L VUAMTOH NUMBER DP 

V>0 

4,92276 

2 1 

V-l 

4.01486 21 

V- 2 

3*3414= 21 

v-a 2.7fr5re 21 v-4 2 , 2 a* it 

21 

Y-5 
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21 

V-5 

1*53836 

21 

V- 7 

1,2499* 21 

V«8 
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V-9 T.02S9E 20 V- 10 5*19016 

2D 

V“l L 

4. 03826 

20 
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20 
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N=l 
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22 
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9-4 

0, 
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N-2 
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2f 

N-l 

4*15936 

19 
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0, 
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M-3 

2,20796 

19 

N-4 

0* 






ISND SINGLET Pi STATES WITH 


N=3 

A * 91 96 E 
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0. 


19= L 4 

0* 
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DETAILED NUMBER DENSITY FDR H2+* 1 ON 5/C Util C METER 
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V= 0 
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22 

V= 1 
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22 
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DELI L. 32586 06 DEL 3 N 


1*32586 

06 

ftO 


1*115/6- 

09 

^ 09 1 * 1 1 6 7E 

-09 


CUJ l 4*5 93/6 00 CLJT2 


2 * 3 7 fr 4 fc. 
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41 

2*51 76E' 

01 

U£ L L 

S4I93 - 

2*1206£- 

07 

oeu 

. Gi ( 63-2*12066 

-07 
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1*00006 

00 

U 


4 


12 
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22 
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25 
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25 

20000 
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24 
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22 
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26 
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24 
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25 
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23 
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21 

6.139E 

23 
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30000 
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7* 02 OE 

25 

7 . 96 36 

24 
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25 
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23 
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23 
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25 

>0000 
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20 
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25 
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24 
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25 
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it 

2*0506 

18 
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22 

8, J95£ 

25 

1*4326 

21 
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25 


First page of ease l. 
^Second page of esse l. 
°First page of cases 2 to 1, 
^Second page of cases 2 to 1. 
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fb) At oressure of L 0B2MG S newtons per square meter 
UOOO atm). 

Figure 2. - Number dens? ties of components of hydrogen 
plasma with coulomb Interactions between free charged 
particles for T > 2000 1C 
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Figure 3. - Calculated number densities of e" and H“ 

In hydrogen plasma at pressure of 1. 01325x10 s newtons 
per square meter (1000 atm} with and without Hjl 
Coulomb interactions between free charged particles 
were included. 



Figure 4, - Ratio of calculated number densities of ft e“, 
and H” with and without coulomb interactions between 
f ree c ha rged par t i c I es. P ress u re of L 01 325x 10® newton s 
per square meter (1OO0 atml. 
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Figure 6. - Simplified flow diagram for computer program for components of hydrogen plasma. 
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— Actual potential energy 

— Actual dissociation asymptote 

— - Dissociation asymptote if orincipal Quantum numbers were 

Independent of r 

Potential energy If principal Quantum numbers were Inde- 
pendent of r 



I nternuclear distance, r 

Figure 9 . - Quantities in approximation for hump height for excited elec- 
tronic state with promoted molecular orbital, I Mot to scale , } 


NASA-X^gtfly, Its® 25 £-4346 
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